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SUMMARY 

The objective of this work was to assemble and analyze fragmentation 
data for exploding liquid propellant vessels. These data were to be retrieved 
from reports of tests and accidents, including measurements or estimate 
blast effects, fragment velocities, masses, shapes, and ranges, ^orrea- 
tions were to be made, if possible, of fragmentation effects with type of acci- 
dent, type and quantity of propellant, blast yield, etc. - signi ican amoun 
of data was retrieved from a series of tests conducted for measuremer 
blast and fireball effects of liquid propellant explosions (Project PY ), a 
few well -documented accident reports, and a series of tests to determine 
autoignition properties of mixing liquid propellants. The data were re nee 
and fitted to various statistical functions. Comparisons were ma e wi 
methods of prediction for blast yield, initial fragment velocities, an rag- 
ment range.' Reasonably good correlation was achieved. Methods presente 
in the report allow prediction of fragment patterns, given type an quan 1 > 
of propellant, type of accident, and time of propellant mixing. However, 
more work must be done before the results of this study can be easily applied 
to estimation of damaging effects of fragments from exploding l'quid pro- 
pellant vessels. 


INTRODUCTION 


Background 

The primary hazard relating to large-scale explosions has « m the past 
been assumed to be the blast wave generated by the explosion. 1 erma 6 ' 
fects have been considered next, and effects of damaging fragments last. 

This study attempts to partially rectify this situation by provi mg a compre 
hensive analysis of fragmentation effects of bursting liquid propellant vessels. 

In storage or in a launch configuration wivhin tankage in a rocket motor, 
liquid propellants are initially contained within vessels of various sizes, ge- 
ometries, and strengths. Various modes of failure of these vessels are^ 
possible, from either internal or external stimuli. If the vesse is presou 
ized with static internal pressure, one possible mode of failure is simp y 
fracture, instituted at a critical size flaw and propagated throughout t e 
vessel. A similar kind of failure can occur if the vessel is accidenta y i 
mersed in a fire, and pressure increases internally because of vaporization 
of the internal propellant. Some launch vehicles have the liquid fuel and oxi- 
dizer separated by a common bulkhead. Accidental over-pressurization of 
one of these chambers can cause rupture of this bulkhead, and subsequen 
mixing and explosion of the propellant. External stimuli which can cause 
vessel vailure include high-speed impact by foreign objects, accidental deto- 
nation of the warhead of a missile, dropping of a tank to the ground (as m 
toppling of a missile on the launch pad), as well as many other external 
sources Vessel failure can result in an immediate release of energy or it 
can cause subsequent energy release because of mixing of propellant and 
oxidizer and the subsequent ignition. Other modes of failure which have re- 
sulted or could result in violent explosions are fall-back immediate y a er 
launch due to loss of thrust, or low-level failure of the guidance system after 
launch with subsequent impact into the ground at several hundred feet per 

second. 

Failure of a vessel containing liquid propellants can result in various 
levels of energy release, ranging from negligible to the full heat value of tne 
combined propellant and oxidizer. Toward the lower end of the scale of 
energy release might be the failure of a pressurized vessel due to crack pro- 
pagation He r “ the stored pressure energy within the compressed propellant 
or gas in an ullage volume above the propellant could accelerate fragments of 
the vessel or generate a weak blast wave. In the intermediate range of energy 
releases could lie vessel failure uy external stimulus and ignition, either very 
rapidly or at very late times, so that only small proportions of mixed pro- 
pellant and oxidizer contribute to the energy release. At the upper end of the 
scale could be the explosion of a mixed propellant in a vessel wherein a pre - 
mixed propellant and oxidizer detonate in much the same fashion as a hig i 
explosive, and explosions resulting after violent impact with the ground. In 
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past studies of possible blast and fragmentation effects from vessel rupture, 
a critical problem has been to accurately assess the energy release as a re- 
sult of the accident or incident. A common method of assessment of possible 
energy release or a correlation of the results of experiments has been to 
assess the energy release on the basis of equivalent pounds of TNT. Ibis 
method is used because a large body of experimental data and theoretical 
analyses exist for blast waves generated by TNT or other solid explosives 
(refs. 1 and 2). Although the comparison with TNT is convenient, the corre- 
lation is far lrom exact. Specific energies which can be released, i. e. , 
energy per unit volume or mass of reacting material, differ quite widely be- 
tween TNT and various liquid propellants or mixtures of liquid propellants 
and oxidizers (ref. 3). 

Dependent on the total energy release and the rate of this energy re- 
lease, the sizes and shapes of fragments generated by liquid bursting pro- 
pellant vessels and their appurtenances cover a very wide spectrum. At one 
extreme is the case of a vessel bursting because of seam failure or crack 
propagation from a flaw wherein only one "fragment" is generated, the vessel 
itself. This fragment, from a very slow reaction, can be propelled by re- 
leasing the contents of the vessel. At the other extreme is the conversion of 
the vessel and parts near it into a cloud of small fragments by an explosion 
of the contents of a vessel at a very rapid rate, similar to a TNT explosion 
(refs. 4 and 5). For most accidental vessel failures, the distribution of fragment 
masses and shapes undoubtedly lies between these two extremes. The modes 
of failure of the vessel may be dependent upon details of construction and the 
metallurgy of the vessel material. Some of the masses and shapes are dic- 
tated by the masses and shapes of attached or nearby appurtenances. In any 
event, assessment and prediction of these parameters undoubtedly is much 
more difficult than is true for the better understood phenomenon of shell 
casing fragmentation. 

Once the masses, shapes, and initial velocities of fragments from 
liquid propellant vessels have been determined in some manner, then the 
trajectories of these fragments and their losses in velocity due to air drag or 
perforation or penetration of various materials must be computed. This pro- 
blem is one of exterior ballistics. It differs from conventional exterior 
ballistic studies of trajectories of projectiles, bombs, or missiles in that 
the body in flight is invariably very irregular in shape and is usually tumbling 
violently. Exact trajectories cannot be determined then in the same sense 
that they can be for well- designed projectiles. Only approximate trajectories 
can be estimated, usually by assuming "equivalent spheres" oi other geo- 
metric shapes for which exterior ballistics data and techniques exist. But, in 
some fashion, one can predict the ranges and impact velocities for fragments 
which were initially projected in specified directions from the bursting vessel 
with specified initial velocities. An example of results of such analysis is 
given by Ahlers (ref. 6). 
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This problem is not complete until one can assess the effect o f frag- 
ments from the burst propellant vessels on various "targets". For a pioper 
assessment of hazards, one should consider a wide variety of targets, in- 
cluding human beings, various classes of buildings, vehicles, anc per laps 
even aircraft. Problems of this nature are exceedingly complex, not on y 
because of the inherent statistical nature of the characteristics of the im- 
pacting fragments but also because the terminal ballistic effects for large 
irregular objects impacting any of the targets described are not very well 
known. In most past studies of fragment damage from accidents, the investi- 
gators have been content to simply locate and approximate the size and mass 
of the fragments in impact areas and have ignored the important problem o. 
the terminal ballistic effect of these fragments. 

Related Work 

Extensive studies have been carried out over many years regarding 
the potential failure of nuclear reactor vessels from a variety of causes. The 
source of energy causing a reactor vessel failure can be the stored compres- 
sive energy in a liquid or gas within the containment vessel, the chemical 
energy release, or the uncontrolled nuclear energy release. The latter source 
is, of course, not present in the failure of liquid propellant vessels, but the 
first two sources are present. Although many of the studies of nuclear re- 
E actor vessels have concentrated on the design of the pressure vessel and the 

» attachments to it to obviate failure, many other studies also have been eon- 

cerned with shock and fragmentation effects in the event tnat failure does 
occur . The literature in this field is far too voluminous to cite other than to 
1 give a few examples which are indicative of the parallels between these studies 

I . and those reported here. The specific references given all relate to produc- 

; t; on of or containment of fragments caused by vessel failure. 

i 

* The first of these references cited is a review paper by Gwaltney ( ref . 

! 7) on missile generation and protection in one class of a nuclear power plant. 

I Various types of vessel failure are considered and reviewed and formulas are 

l j given for estimations of velocities with which fragments will be ejected after 

! > failure. Effects of impact of the fragments are discussed, and a number of 

! : empirical penetration formulas for metal missiles penetrating and perforating 

f i steel an J > oncrete are given. The paper also includes a simplified discussion 

; i cf poss-hip shock wave effects caused by the release of energy after vessel 

« rupture. The second paper is reference 8. This is the final report of a multi 

year experimental and analytical investigation by Stanford Research Institute 
of the problems of generation and containment of fragments generated by a 
runaway reactor. More details of various phases of the investigation are 
l given in additional 6-month progress reports predating reference 8. Many of 

I the aspects of the work reported in reference 8 are similar to the study re- 

I 1 ported here. Attempts were made to simulate the energy release rates in the 

■ event of reactor runaway by use of slow-detonating explosives and fuses. 
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These sources of energy release were used to apply pressure within model 
containment vessels and models of surrounding materials such as concrete. 
Failure of these model vessels was observed using high-speed cameras to 
determine velocities and initial trajectories of the fragments. In a parallel 
investigation, the Stanford Research institute staff conducted a series of ex- 
periments simulating impacts by long, slender missiles such as reactor con- 
trol rods. Penetration formulas for such rods striking steel plates were 
developed as part of the effort. 

A good general listing of the classes of problems considered in nuclear 
reactor containment studies can be found in reference 9, which reports papers 
on reactor safety given at the Second International Conference on Peaceful 
Uses of Atomic Energy. In particular, reference 10, one of the papers in the 
Proceedings of the conference, discusses various sources of energy release 
and gives approximate limits to the magnitudes which can be expected, con- 
siders ways of attenuating blast energy and of stopping fragments, and gives 
in general a good overall review of the range of problems one must consider 
in reactor containment studies. 

The explosive behavior of bombs, grenades, mines and warheads has 
always commended wide attention, and the most commonly used bombs are 
usually constructed from suitably corrugated steel casings, either fully or 
partially filled with explosives. Interest in the mechanics of fragmentation 
has largely been directed towards trying to predict the influence of casing 
material and wall thickness, the size of the explosive charge, and the type of 
explosive on the fragmentation velocity and the expanded radius of a casing at 
which fracture occurs. 

In a series of papers (refs. 4, 5, 11 through 14), there has evolved a 
simple approximate treatment for the acceleration of fragments by high ex- 
plosives. The basic assumption made was that the potential energy in the 
charge before detonation was equal to the kinetic energy of the charge and 
casing after detonation and expansion. It was also assumed that, after detona- 
tion, the gaseous detonation products were equally dense at all points and ex- 
panding uniformly. Formulas for fragment velocity at the radius for case 
breakup (essentially the maximum velocity) are presented in these references 
for various regular geometries of cased explosive charges. All are of the 
form 


U - f(H e , M/C) (1 

where J is velocity and is a function of heat of explosion H e , total casing 
mass M, and mass of the explosive charge C. 

If the results of accidents invoHag explosions of liquid propellant 
vessels are well documented, they can provide useful data to assist in the 


4 


assessment of this problem. Some have indeed proven useful sources for our 
study, as we will document in later sections of this report. Although acci- 
dent reports are useful in documenting the gross effects of vessel explosion 
determining the maximum ranges to which fragments are projected, and indi- 
cating shapes and masses of fragments, they are often of less value in assess- 
ing this problem than are controlled experiments. Because they are accxden a. 
usually there is no measure of blast yields, fragment trajectories, and oii.r 
data that would be useful in analysis of vessel failure and subsequent effects. 

Project PYRO involved many test explosions with liquid propellants. 

The purpose of Project PYRO (refs. 15 through 17), was "to develop a reliable 
philosophy for predicting the damage po+ential which may be experience rom 
the accidental explosions of liquid propellants during launch or test opera- 
tions of military missiles or space vehicles". Three combinations of propel- 
lants and oxidizers were chosen for test and evaluation, and at least seven 
agencies were involved in the program. The primary objective was to est 1 
mate blast yield and its effects. The effects of fragmentation were secondary 
in the study. But, Jeffers (ref. 18) analyzed a small number of the photo- 
graphic records to determine fragment velocity. As is apparent in later sec- 
tions of this report, the films from Project PYRO, when studied carefully, 
provide the primary source of data for initial velocities for liquid propellant 

explosions. 

There are a number of experimental studies and analyses of the ef- 
fects of bursting pressure vessels which fail under the action of internal 
energy sources other than liquid propellants. A number of these can provide 
useful information for the problem at hand. Some specific examples follow. 


Hunt, Walford, and Wood (ref. 19) have conducted an experimental 
study of the failure of a pressure vessel containing high temperature pres- 
surized water. In this study, the authors observed the failure of a vessel 
with high-speed cameras and also located a number of blasi transducers near- 
by to measure the resulting shock wave generated in the surrounding air. 

They also generated equations for calculation of velocities of the fragments 
resulting from burst of the vessel. In a somewhat similar study, Larson and 
Olson (ref. 20), measured the air blast effects from bursting pressure vessels 
containing high pressure gas. In this study, the authors also observed t :e 
flight of fragments from the bursting vessel and developed an empirical 
method for estimating fragment velocities based on an energy balance and 
knowledge of the strength of the shock wave generated by the bursting vessel. 
In a somewhat different category than the two previous studies are analyses 
and predictions of the effects of rupture of pressure vcisels containing hign 
pressure gases. An excellent example of such analyses is a compendium of 
gas autoclave engineering studies, edited by C. E. Muzzall (ref. 21). This 
compendium is an exhaustive study of the possible hazards associated with 
failure of a large, cylindrical vessel containing high pressure, high temper- 


ature argon. Estimates are made of both the blast and fragmentation hazaids 
in the event of failure of the vessel, of the effects of both fragments and blast 
on a test cell within which the vessel is operated, and recommendations fo - 
redesign of the test cell to withstand both blast and fragmentation effects. A 
second study of this same nature, but on a much more limited basis, is a re- 
port by Baker, et al. , (ref. 22), on the possible effects of failure of a high 
pressure helium vessel while under test in a NASA vibration and acoustic test 
facility. H 'e, blast and fragmentation effects were estimated in the event of 
failure of the vessel, blast loading and response of the walls of the test fa- 
cility were computed, as were possible penetration effects by fragments of the 
vessel. The report concluded with recommendations for modification of vest 
procedures to obviate the very real hazards in the event the helium pressure 
vessel failed. 


Purpose of Present Work 


The purpose of the work reported here is to assemble fragmentation 
data for bursting liquid propellant vessels, analyze these data, and develop or 
modify methods of prediction of fragmentation effects of such explosions;. An 
additional purpose is to enter all relevant reports, data, etc. , into a data bank 
at the NASA Aerospace Safety Research and Data Institute (ASRDI). 

Scope of Present Work 

The primary purpose of this study and analysis is the retrieval, assem- 
bly and recording of available data regarding fragments from exploded vessels 
that contained liquid propellants or substances that have similar chemical 
properties. These data cover both test and accidental explosions and include 
blast yield, fragment sizes, fragment trajectories, fragment velocities, and 
a description of damage caused by the fragments. 

The second part of the study and analysis shall be that of reducing the 
available data to a form most readily usable by aerospace engineers, for es- 
timation of the integrated shrapnel hazard to which neighboring structures 
will be subjected. This reduction of the da*a shall be in the form of equations 
that will relate the blast yield to the nature and quantities of fuel and oxidants 
and to various parameters describing the type of explosion, and in the form of 
equations that will relate the blast yield to distributions of fragment s.ze, 
initial fragment velocity, and initial direction of fragment motion. 

Significance 

It is believed that this report contains the first comprehensive 
study of fragmentation effects of bursting liquid propellant vessels. The re- 
sults of this work should provide a better understanding of such fragment 
characteristics as initial velocity, mass, shape, and range as thev relate to 


6 




I 


\ 


estimated blast yield of exploding liquid propellant tanks. These charac- 
teristics used with munition fragment equations could predict the termina 
or impact velocities of these fragments. Thus, it is possible to derive cri- 
teria that allow for the prediction of fragment hazards to people or the risk 
of damage to nearby facilities from the impact of these fragments. These 
criteriacould be used to arrive at safe distances between populated areas 
at launch or test stands for fragments and overpressure hazards that are 
caused by exploding liquid propellant tanks. Moreover, intra-line separation 
distances can be established that would decrease the risk of damage to near- 
by facilities or systems from fragment impact; or, with a prediction of im- 
pact energy of fragments, barriers could be properly designed to protect 
these facilities at intra-line separation distances. 


Statistical fitting to data on fragment range R versus measured 
terminal blast yield Y gave the following equation: 


zv 0.2775 

ft = 314. 74 Y 


where R is the mean range in feet, and Y is the terminal blast yield in 
percent TNT equivalent. Fitting to data on fragment weight W and mean 
presented area A gave the equation 


/v o / ? 0. 78 

ft = 9.864 (A/W^ 3 ) 


where W is fragment weight in pounds and A is mean fragment presente 
area in square inches. These two equations can be used for prediction sub- 
ject to restrictions and limits noted in the body of the report. Also included 
in Part IV are distribution functions for fragment initial velocities for var- 
ious types of accident, which can be used to predict distributions of fragment 
sizes and masses for other postulated accidents. 


Estimates of the distributions of the initial velocities for lour sped u. 
combinations of conjurations and propellants were derived. In n . 
cases, the initial velocity (U > in ft/sec. followed a lo« normal d.st. .button 
function. The distribution functions for the four cases are S ,ven uelow. whet, 
ft and* 5 are estimates of means and standard deviations of U, rospectlt ■- h 


1. Confined by Missile (CBM). LO^LH, propellant 

A A r 


f (U.; u, ") - (1/2- 475 3 C. ) exp [- ( 1 n U . - 6 . 464) 2 / 1 . 9503 )] 


Z. CBM, LO /RP-1 propellant 


f ( U ; 4, " ) = (1/1.58 24 IT) exp (In U . - 6.713) / 0.1° t l] 


m 


3 . Confined-by-Ground-Sarface (CBGS), LO^LH., propellant 


A A 

f ( M, a ) 


= (1/1. 9339 IL) exp [ - !inU.-6. 129) 2 / 1. 1904] 


4. CBGS, LO^/RP- 1 propellant 


f (U ; M, (J) = (1/1.6010 U.; exp I - ( n U. - 5 . 962) /0. 8 159] 

i i L i 
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I. RETRIEVAL OF FRAGMENTATION DATA FOR 
LIQUID PROPELLANT VESSELS 

The first task in this contract consisted of a series of contacts and 
visits with various government agencies and contractors to ascertain the ex- 
tent of data available on fragmentation from liquid propellant explosions, 
either accidental or from planned tests, and to obtain pertinent data and re- 
ports for entry into the data bank at the Aerospace Safety Research and Data 
Institute (ASRDI) of NASA. The contacts and visits were supplemented by a 
conventional literature search of the open literature and the Defense Docu- 
mentation Center (DDC). 

The work commenced with an initial visit to ASRDI, and temporary 
trarsfer to SwRI of pertinent documents already acquired by the ASRDI staff. 
Potential sources of data and individuals to contact at various agencies, m 
audition to those already known to SwRI staff, were also identified during this 
visit. We then made a series of telephone contacts to determine whether 
specific agencies or firms had applicable data and could send such data to is, 
and if a visit was desirable. All major NASA centers, several AEC labora- 
tories the military service ordnance laboratories. Air Force Eastern and 
Western Test Ranges, the Department of Defense Explosive Safety Board, and 
a number of commercial and other contractors were contacted during this 
initial telephone survey. More than thirty such contacts were made. A num- 
ber of the calls led to blind alleys, with no data available, or individuals who 
might have had data or known of it being no longer present. But, other calls 
unearthed potential sources of data and allowed appointments for visi' s to 
review these data. 

Following the initial phone contacts, several SwRI staff members 
visited those agencies which were potential sources of data. As much as 
possible, trips were combined to agencies in the same general geographical 
area. W'hen the visits yielded applicable or potentially applicable documents, 
reports or data, we tried to obtain them for permanent retention or loan, or 
tried to arrange for them to be transmitted to SwRI. The results of these 
visits are summarized in Table I. which lists the agencies visited, individuals 
contacted at each agency, and the type of applicable data found to be available 
to us at each agency. Several agencies had sufficient data or information to 
warrant a follow-up visit to further discuss the data or to attempt to obtain it 
for use on this contract. These agencies are indicated by an asterisk in 
Table I. Of particular importance to this contract is the library ot motion 
pictures of the Project PYRO tests available at Air Force Rocket Propulsion 
Laboratory. We obtained these films, and they form the primary data base 
for determination of initial fragment velocities from bursting propellant 
vessels. Individuals visited at the various agencies were for the most part 
very cooperative and helpful in obtaining or transmitting applicable data and 
documents. In at least one instance, however, we have not been able to obtain 


TABLE I - SUMMARY OF AGENCIES VISITED TO OBTAIN 
FRAGMENTATION DATA OR DOCUMENTS 


Agency 


Individuals 

Contacted 


Data Available 

Docu- Accident j 

ments Reports J Films Other 


Aerospace Safety Re- 
search and Data r 
tute, NASA- Lewis 


NASA- Kennedy’ 


I. 1. Pinkel 

C. D. Miller 
P. M. Or din 

D. Forney 

J. H. Deese 
A. H. Moore 
F. X. Hartman 
A. J. Carraway 


Air Force Eastern 
Test Range" 


L. J. Ullian 
CPT R. P. 

Welborn 
T. Fewell 
Maj. K. Cailer 


Dept, of Mech. Eng., 
Univ. of Florida 


Prof. E. A. 

Farber 
Prof. E. Watts 
Prof. J. Smith 


Air Force Rocket 
Propulsion Lab. 
Edwards AF Base 


J. G. Wancheck 
R. Thomas 


Aerospace Corp., 
El Segundo 


Air Force Space and 
Missile Systems Org. 

Gen. Elec. Co., 

Bay St. Louis, Miss. 


R. Wolfe 
J. Smith 
R. Vega 

CPT K. C. 
Tall man 

P. V. King 



This agency was v 


isited twice to obtain data identified during the irst visit 


0 



TABLE I. (Continued) 


Agency 


Individuals 

Contacted 


Dat d A va 1 1 ab l e 


Docu- 

Accident 


ment s 

Reports 

Films 

h 


t NASA- Mar shall 

[ 

Central Propulsion 

Info Agency (CPU), 

Johns Hopkins Univ., 

App 1 ied Physics Lab 


Dept, of Defense Ex- 
plosive Safety Board 
( DDESB) 


I URS Research Co. , 
j San Mateo. California 

i 

j U. S. Army Ballistic 
' Research Labs. 


W. A. Riehl 


E. A. Cad- 
wallader 


X 


X 


Dr. T. A. Zaker 


X 


R. Pe rkins 


C. Wilton 


X 


I D. Dunn 
| C. Kingery 
| A. J. Hoffman 


X 


j 

t 




copies of accident reports which would provide useful fragmentation data, and 
have not been able to include these data in our review and further analysis. 


As documents and data were received at SwRI as a result of our initial 
visit to ASRDI, our subsequent visits to other agencies, and our library and 
DDC literature searches, we reviewed each document, completed ASRDI 
form 102 A for the documents and forwarded these forms to ASRDI. A total of 
168 documents were reviewed in this manner, with various SwRI staff mem- 
bers completing the Forms 102A for documents which fell within their tech- 
nical specialties. 


We believe that we have discovered and reviewed most of the pertinent 
literature, data, and accident reports pertaining to fragmentation of liquid pro- 
pellant rockets and vessels. There is, however, one possible exception. 

There may be a body of fragmentation data in accident reports in the Air Force 
Inspection and Safety Center at Norton Air Force Base, California, which we 
could not review or obtain for legal reasons. These reports were reviewed by 
staff members of the Center, and we have been notified that they contain no 
data which could be used in our study. Because we were not allowed to re- 
view the reports ourselves, we have no way of comparing them with other 
accident reports which have provided useful data, nor do we know the criteria 
applied by Norton personnel in assessing the potential value of specific reports 
to this project. 
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II. 


DETERMINATION OF BLAST YIELD 
A. General 


A prerequisite to estimation of fragmentation effects xor liquid prope - 
lant explosions is the estimation of energ' released during the exp OS1 ° n ’ 
which is synonymous with explosive yield. Furthermore, to pr°perl> e-L n 
fragment Velocities of appurtenances which can be accelerated by the bias 
wave from propellant explosions, one must know the time hl8t °'- 8 a ° f c ; a f ”^ S 
physical parameters describing the blast wave as a function dist 
the explosion. We must therefore consider blast effects in some 
though this is a study of fragmentation. 

Accidents with liquid propellant rockets, both during static finn 0 on a 
test stand and during launch, have shown that liquid propellants can genera e 
violent explosions. These explosions "drive" air blast waves which can 
cause direct damage and can accelerate fragments or nearby objects. The 
ta" eh at the Air Force Eastern Test Range (ETR) have for a number of 

years been instrumented with air blast recorders to measure the 
sures generated during launch pad explosions, so some data are 
the intensities of the blast waves generated. Such measurements, an e 
common practice in safety circles of comparing explosive effects on the basi. 
of blast waves generated by TNT, have led to expression of blast yields of 
propellant explosions in ec.uiva.ent "pounds of TNT" (Al.hougn a direct c 0 , 1 - 
version of pounds of TNT to energy can easily be made - 1 »> m o. TNT equal. 

1 . 4 x 1 O* 5 ft -lb - this is seldom done. ) 

Liquid propellant explosions differ from TNT explosions in a number of 
wavs, so that the concept of "TNT equivalence" quoted in pounds of TNT is 
from exact. Some of the differences are described below. 

(1) The specific energies of liquid propellants, in stoichiometric 
mixtures, are significantly greater than fov TNT (specific 
energy is energy per unit mass). Table II, taken from ref. L 
gives specific energies for a number of hquid-propellant 
oxidize r mixtures, as ratios to TNT specific < nergy Note that 
all of the energy ratios in Table II are greater than 1, and 
range as high as 5. 3. 

(2) Although the potential explosive yield is very high for liquid 

or opel 1 ants . the ac.ual yield is much lower, because propellant 
and oxidizer are never intimately mixed in the proper propor- 
tions before ignition. 
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BLE II - DETONATION ENERGY EQUIVALENTS FOR 
SELECTED LIQUID PROPELLANTS (REF. 3) 


Specific 

Energy 




Oxygen/Fuel 

Ratio 

Density 

(gm/cm 2 ) 

(Relative 
to TNT) 

Fuel 

Oxidize r 

Aluminum Triethyl 
Oxygen 

A1 (C, H ) 
°2 ‘ 

1.5 

1.019 

2.61 

Fuel 

Oxidizer 

Aluminum Trimethyl 
Oxygen 

A1 (C H ) 
°Z 

2. 0 

1.012 

2.80 

Fuel 

Oxidizer 

Pentaborane 

Oxygen 

Nitrogen Tetroxide 
IRFNA* 

B of 9 

n ;° 4 

^1. 14 °3.24 

2.35 

3.35 
3.35 

1. 075 
1. 220 
i. 352 

2. 29 
3.36 
1.92 

Fuel 

Oxidizer 

Monomethylbydrazine 

Oxygen 

Nitrogen Tetroxide 
IRFNA 

Chlorine Trifluoride 

CH..NH NH 
°2 

N 2°4 

c ^,.H°3.Z4 

1.00 

2. 17 
2.45 

3. 13 

1. 007 
1. 253 
1. 460 
1.578 

?.. 1 1 

2. 38 

3. 14 
2. 3 7 

F uel 
Oxidizer 

UDMH 

Oxygen 

Nitrogen Tetroxide 
IRFNA 

« CH 3»2 N Z H Z 
°Z 

N ° 4 

*™1. 14 °3. Z4 

1.70 
2.55 
2. 85 

1. 009 
1. 246 
1.365 

2. 85 
2.33 

3. 16 

Fuel 

Oxidizer 

Hydrogen 

Oxygen 

Nitrogen Tetroxide 

N 2°4 

5. 00 
5. 25 

.968 
1. 2il 

5. 30 
3. 87 

Fuel 

Oxidizer 

Ammonia 

Oxygen 

Nitrogen Tetroxide 

NH 

° 2 

• : 2 °4 

2. 0 
6. 0 

.964 
1. 312 

3. 19 
1. 57 

F uel 

Oxidize r 

Hydrazine 

Oxygen 

Nitrogen Tetroxide 
IRFNA 

n 2 h 4 
% 

N 2°4 

H ^l. 14 °3. 24 

1. 00 
l t 30 
1. 50 

1. 070 
1. 245 
1. 342 

2. 78 
2.36 

3. 17 

Fuel 

50% N_ H -50% UDMH 

_ 2 4 

Oxygen 

Nitrogen Tetroxide 

1. 875 N 2 H 4 +(CH 3 
°2 

N 2 °4 

2. 0 

1. 043 
1. 262 

2. 70 
2. 40 


* Inhibited Red Fuming Nitric Acid 
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(3) Confinement of propellant and oxidizer, and subsequent effect 

on explosive yield, are very different for liquid propellants and 
TNT. Degree of confinement can seriously affect explosive 
yield of liquid propellants, but has only a secondary effect on 
detonation of TNT or any other solid explosive. 

(•) The geometry of the liquid propellant mixture at time of ignition 

can be quite different than that of the spherical or hemisphenca 
isometry of TNT usually used for generation of controlled blast 
waves. The sources of compiled data for blast waves from TNT 
or Pentolite such as references 2 and 3, invariably rely on 
measurements of blasts from spheres or hemispheres o 
plosive. The liquid propellant mixture can, however, be a 
shallow pool of large lateral extent at time of detonation. 

(5) The blast waves from liquid propellant explosions show differ- 

ent characteristics as a function of distance from the explosion 
than do waves from TNT explosions. This is undoubtedly 
simply a manifestation of some of the differences discuss ec 
previously, but it does change the "TNT equivalence' of a 
"liquid-propellant explosion with distance from the explosion. 
Fletcher (ref. 37) discusses these differences and shows them 
graphically in Figs. 1 and 2. These differences are very evi- 
dent in the results of the many blast experiments reported in 
Project PYRO (refs. 15-17). They have caused the coinage of 
the phrase "terminal yield", meaning the yield based on blast 
data taken at great enough distance from the explosion for the 
blast waves to be similar to those produced by TNT explosions. 
At closer distances, two different yields are usually reported, 
an overpressure yield based on equivalence of side-on peak 
overpressures, and an impulse yield based on equivalence of 
side-on positive impulses. 

There exist at present at least three methods for .'stimating yield from 
liquid propellant explosions, which do not necessarily give the same ’ 

tions. One method is based on Project PYRO results (refs. 15-1 0. and the 
other two are the "Seven Chart Approach" and the "Mathematical Model of 
Farber and Deese (ref. 23). We next discuss each method and some back- 
ground information. 
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Project PYRO and Related Experiments 


Project PYRO was a joint NASA / USA F project 
period 1965-1967 with the purpose of determining the 
acteristics of three liquid propellant combinations m 


conducted during the 
blast and thermal char- 
most common use in 




military missiles and space vehicles. It included 270 tests w*th .o.al heights 
military missile. F fo i 0 0 000 lb . Most of the tests were 

of propellants ranging from 200 lb to 100. UUU H n1 pn , Fdwards 

conducted at the Air Force Rocket Propulsion Laboratory (AIRPL) a, Edwards 
AFB, California. Prime contractor for much of the effort was L . V* 11 
Corp. , Burlingame, California. The project was supervised by a Stein 
Committee of representatives from several NASA centers, tit ir 
Eastern Test Range, and the Sandia Corp. 

The emphasis in Project PYRO was almost exclusively experimental ^ 
Tests were designed to simulate various types of accidents v. it. i cou c t. t 
mixing and ignition of the propellants. The primary blast instrumc n a^ion 
an array of blast pressure transducers whose outputs as funt ions 
recorded on magnetic tape. (Although fragmentation effects were int iae < 
the program objectives, high-speed motion picture cameras pnotographed m 
of the tests, and our data on fragment velocities are all obtained from most 
films ) The results of the program were reported in a three-volume 
report, with Vol. I (ref. 1 5 ) de scribing the program and giving oveial J" - 
Vol II (ref. 16) giving detailed test data, and Vol. Ill tret. ' g lxlI .h 
«:» IL- based orf the program re.uh., U «hc PYRO ellor. . - - 

n-pes of accidents were simulated. The first type constated of fatlur, of . 
interior bulkload separating fuel and oxidizer in a missile slag, . = 

termed Confinement by the Missile (CBM). The second type of «tH.« 
included impacts at various velocities of ,h missile on the ground, with all 
tankage ruptured, and subsequent ignition. This was termed Con it., mm 
the Ground Surface (COGS). The third type was High Velocity Impact « 1 U) 

after launch. 

Although Project PYRO generated much more daca on explosive yieh.* 
of Ucuid propellant explosions than all previous studies combined 

r Her oLorimental programs did give useful data and should be n-n -n 
here. 7 rthur D. Little, Inc., (ref. 24) conducted a series of blast tests sim- 
ulating spills and ignition on the ground of various combinations o pu e l*n, - 
in ‘■he Saturn vehicles. The tests were designed to produce the maximum 
possible blast yield for this type of accident, with thorough mixing o i u s -■ ml 
oxidizer and delays of ignition until such mixing was complete. - 

men;' investigators reported bl - :c wave characteristics ^entical to 
T\T explosions in the overpressure range of their measui « ■ ' ^ j 

and blast yields ranging from 0. 21 to 1. ‘>8 lb TNT /lb propellant They also 

estimated maximum potential yields (even including effect s o! ‘fierburning oi 
unburned fuel with oxygen in the air) of considerably less than in -ible II. 
Their predicted maxima are given in Table III. 
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TABLE III 


PREDICTED MAXIMUM 

BLAST YIELDS (REF. 24) 

Propellant 

Ratio of Blast Energies, 

lb TNT /lb Propellant 

RP-1 /LO z 

1.25 

lh 2 /l° 2 

3.70 

RP-1 /LO /LH 

2.75 


Another experimental effort prior to PYRO is reported by Pesante and 
Nishihazashi (ref. 25). These investigators measured the blast waves gener- 
ated by fuels and oxidizers which were violently mixed by explosively s a - 
tering dewars containing one component, while the dewars were ™ . 

a bath of the other component. Blast yields ranging from 0.-3 
lb propellant were obtained in these experiments. (These tests a so inc u 
attempts to measure velocities of objects placed near the blast wave, but no 
useful data were obtained.) 


The final set of tests prior to PYRO is apparently a series reported by 
Gayle, et al. (ref. 26). Fuels and oxidizers were mixed by several different 
methods (anticipating the two primary simulation methods in PYR ), an as 
wave properties measured as a function of distance. These investigators 
showed a much greater spread in blast yields for other methods of mixing than 
spill tests, with much smaller yields being observed in most tests for simu- 
lated bulkhead rupture, etc. Yields for LH 2 /L0 2 combinations were most 
affected by the change in methods of mixing, being no greater than 0.014 lo 
TNT / lb propellant for any of the tests . 


From the test results reported in references 15 and 24 through 26, a 
number of observations can be made regarding blast yields; from liquid propel- 
lant explosions . 


The yield is very dependent on the mode of mixing of fuel and 
oxidizer, i.e., on the type of accident which is simulated. 
Maximum yields are experienced when intimate mixing is 
accomplished before ignition. 
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i m dA, 


(2) Blast yield per unit mass of propellant decreases as t t p 
pellant mass increases. 

(3) The character of the blast wave as a function of dist “ c ' 
between propellant explosions and TNT explostons, as noted 
before There is some evidence that these differences are 
greatest for low percentage yield explosions. (They were not 
observed in tests of ref. 24, for example.) 

m On many of the LH ? /LO z tests (regardless of investigator s ) 

H) Spontaneous igniting occurred very early in the munng process, 

resulting in very low percentage yields. 

(5) Yield is very dependent on time of ignition, even ignoring the 
possibility of spontaneous ignition. 

(6) Yield is quite dependent on the particular fuel anci oxidizer 
mixed. 

(7) Variability in yields for supposedly identical tests was great 

' compared to variability in blast measurements of conventtonal 

explosives. 

The PYRO blast yield prediction methods given in referenc 

i^cuuTy^rAruirp^- 

S' 

CBM) confinement by the ground surface (CBGS), and high velocity impact 
(CBM), coniinemeiu s . , , of overpressure 

(H VI). Equivalent TNT yields are , ^YaTt data for TNT (ref. 2), with a cor- 

propellant explosions. Unfortunately . »e feel th« the P^'ction ^ 
given in reference 17 are oversimplified (for example, tney 
different correction factors for impulse 

tances whereas the data in references 15 and lb show a coni nu 

jh,. 

M^ef^ Km “*- 

fin of blast yield by early autoignition when large masses of propellants art 
mixed is ignored in the PYRO prediction schemes Considering th e care ful 
and well-documented experimental work reported tn references 15 and In. 
prediction methods of reference 17 arc quite d.sappo.nt.ng. 


C. Work of Farber and Deese 


More or less concurrent with the PYRO work, but continuing to the 
present, Farber at the University of Florida and Deese at NASA-Kennedy have 
conducted a combined theoretical and experimental program of the study of the 
physical and chemical processes involved in mixing, ignition, ~nd - plosion o 
liquid rocket propellants. The results of this work are reported in a number 
of papers, with the efforts through October 1968 being best summarized and 
reported in reference 27. Later work is reported in references 28 through 30. 


In the work reported in references 28-30, the problem o mixing, igni- 
tion and explosion of liquid rocket propellants is subdivided into a number of 
sub-problems, and each sub-problem studied more or less independently. The 
results of the subsequent analyses and experiments are then combined in sever 
al prediction schemes for explosive yield, of which the most detailed is the 
"Seven Chart Approach" of reference 23. 


The sub-problems into which the overall problem was divided by Farber 
and Deese were: 

(1) Determination of the potential maximum explosive yield obtain- 
able if the liquid propellants present were mix 2 d in an optimum 
manner (yield potential function) 

(2) Determination of amounts of propellant which would be mixed as 
a function of time after "spill" (mixing function) 


(3) Determination of most probable times of delay of ignition and 

detonation (delay and detonation times) 


Farber and Deese (ref. 23) also evolved a method of empiric ally fitting to 
experimental data a four- parameter probability function vlieh would predict 
the probability of explosive yield to various levels ol conf.dence, and have 
more recently evolved an hypothesis of a critical mass of nixing propellants 
for which autoignit on is certain. 

The yield potential function is calculated by Farber, ct al (ref. 27) on 
th. basis of chemical kinetics considering boiling and freeing of fuel-ox:.dizer 
components as a function of time after an assumed ir startaneoc s mixing. 

Meat values for the various chemical reactions which coi Id occur at various 
times, considering states and amounts of reactants present, are then calcu- 
lated. Some details of the manner in which these calcu’. itions are made are 
given in reference 27, and a typical result for an initial mixture oi L0 2 /LH 2 / 
RP-1 is shown ir. Figure 3. There are no experimental data which directly 
confirm such theoretical calculations of the yield function. 
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The mixing function, expressing the fraction of mass of fuel plus oxidi- 
zer actually mix( d as a function of time after spill, is the best established of 
Farber's sub-problems. No less than four experimental methods were used to 
establish this function (ref. 27). A typical mixing function is shown in 
Figure 4. This function peaks as mixing becomes more complete, and then 
decays because the mixed propellant evaporates. An optimum time for ignition 
(optimum in that it will produce maximum explosive yield) is therefore heavi y 
dependent on this function. 

The least well-determined of the sub- problems is the definition of ex- 
pected delay times for ignition. A statement from reference 28 apoears ap- 
propriate "The ignition time for prediction purposes, can be a controlled 
value, a known value based upon the characteristics of the propellants, a 
statistical value with confidence limits, or it can be a value determined by the 
critical mass method. . . . ". In other words, the ignition time is apparently 
anyone's best guess. It is clear from the results of liquid propellant mixing 
tests that autoignition always occurs for mixing of sufficiently large quantities 
of propellants Farber (refs. 27-30) has hypothesized that a source ot ignition 
whi ;h is always present upon mixing is electrostatic build-up of voltage and 
subsequent electrical discharge through a gas bubble, and that a critical mass 
exists for a given propellant mixture and set of initial conditions which pro- 
vides a short upper limit on ignition time, and therefore an upper limit on 
explosive yield. Recent experimental work by his group (ref. 29) is directed 
specifically toward measurement and verification of this hypothesis, and veri- 
fication is claimed (although not conclusively proven) by work reported in 
re :e re nee 30. 

The two prediction methods developed by Farber, et al, are termed the 
"c even Chart Approach" and the "Mathematical Model". Each will give an 
estimate of explosive yield y , expressed as a fraction of the total heat of 
combustion of a stoichiometric mixture of fuel and oxidizer. In the "Seven 
Chart Approach", a graph such as Figure 3 is normalized by dividing by the 
maximum heat value reached by the mixture, and then converted to a normal- 
i .ed yield potential y versus time plot such as Figure 5. The fraction 
nixed x as a function of time (Figure 4) is then multiplied by y p at approxi- 
mate times to give expected yield y (see Figure 6). Finally, some estimated 
gnition time, with suitable confidence limits, is superimposed on Figure o to 
/ivo the final estimate of yield y . (Note that only four charts arc discussed 
lere. The remaining three charts in the seven-chart method show intermedi- 
ate steps). 

The "Mathematical Model" method (ref. 28) consists of fitting to exper- 
imental data a relationship between normalized yield y and mixing function x 
of the form 
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I I GU Rr' 5. YIELD POTENTIAL AS TIME FUNCTION 





FIGURE ACTUAI. YIELD FOR RANDOM IGNITION AND DETONATION (R EF. 



( 2 ) 


y - 


b + c 


where b, c. and d are parameters to be adjusted. An auxiliary function is 

also introduced which inserts a fourth parameter, a . Using the physical y 
realistic limits of zero yield for zero spill, and a maximum for y w en 
x - 1 which is defined by the fraction of mixed propellants which represent 

stoichiometry (y< 1), fixes the parameters b and c. Farber origina yuse 

data from reference 24 to estimate the parameter d , and has not change * 
parameter since. So, essentially only the parameter a remains to be vane , 
and Farber claims that the parameter a represents a "scaling parameter 
related tototal mass of propellants. With parameters chosen from fitting data 
of reference 24, Farber claims (ref. 28) prediction of upper bounds on explo- 
sive yield which cover all available data through 1969 (see Figure (,. 

After review of the work of Farber, De.se and co-workers, we feel 
that their efforts have their strengths and weaknesses, just as does t e 
work. The greatest strengths are the excellent physical insight into the com- 
plex processes which occur during mixing and ignition of liquid propellants, 
and the division of the complex overall problem into sub-problems which can e 
studied separately. Farber and co-workers were also apparently the first to 
realize that explosive yields for large quantities of propellants were always 
limited by earlv autoignition. The primary weaknesses lie in lack of experi- 
mental verification of the physical processes, sometimes doubtful claims of 
usefulness or applicability of limited test techniques or data . and reiteration 
of the same material in succeeding reports. The two methods for prediction 
of blast yield are well described and understandable, but both give an estimate 
of explosive source energy without consideration of the nature of the blast 
waves generated by liquid propellant explosions which were evident in Project 
PYRO. It is also not clear how the critical mixing function such as r igure 4 
is obtained for various types of full-scale accidents or teste, or how it is 
scaled from laboratory experiments. 

D. Fstimation of Blast Wave Propert ,-s 

We present here methods of estimating blast yields and blast wave 
prop-rtict for liquid propellant explosions. based prin ari',- on PYRO result, 
and on tho work of Farber and Dense. Although our pred.et on method, retain 
many of the features of the previous work, they also differ somewhat where w ■ 
feel changes are appropriate. Furthermore, factors which appear to have t 


In particular, response times, 
phe nomena from thcr rr.ocouple 
seem doubtful. 


time resolution and identification of rhvsical 
grid measurements claimed ir. -eference Z: 


Zb 


EXPLOSIVE YIELD, 





secondary effect on blast yield, such as L/D ratio of tankage, are ignored. 
The concept of "TNT equivalency" is used only to estimate energy of a liquid 
propellant explosion, and not to predict detailed blast wave characteristics . 
Blast is strongly dependent on type of propellant, type of simulated accident, 
impact velocity, and ignition time, so these factors must be accounted for in 
estimating blast wave characteristics and yield. 

Throughout the PYRO work, blast yield is expressed as percent yield, 
based on an average of pressures and impulses measured at the farthest dis- 
tance from the source when compared to standard reference curves (ref. 
for TNT surface bursts (terminal yield). Hopkinson's blast scaling is used 
when comparing blast data for tests with the same propellants and failure con 
ditions, but different mass of propellant. So, the blast parameters P (peak 
side -on overpressure) and I/W 1 ' 3 (scaled impulse) are plotted as functions 
of R/W 1 / 3 , (scaled distance) after being normalized by the fractional yield. 
This procedure is equivalent to determining an effective weight of propellant 
for blast from: 


W 



(3) 


where W T is total weight of propellant, Y is terminal blast yield in percent, 
and W is effective weight of propellant. Because the data are normalized by 
comparing to TNT blast data, the effective blast energy E can be obtained by 
multiplying W by the specific detonation energy of TNT, 1.4x10 ft lb/lb m . 
We will use smoothed curves through the scaled PYRO blast data, and Equa- 
tion (3), to obtain blast wave properties for any particular combination of 
propellants and simulated accident. We will consider each propellant combina- 
tion separately. 

1 . H v per golic Propellant - The hypergolic propellant in widest use, 
and used in the PYRO tests, is a fuel of 50<7c N 2 H 4 -50% UDMH and an oxidizer 
ol Ni O, in a mass ratio of 1/?.. Hypergolic materials, by definition, ignite 
spontaneously on contact, so it is not possible to obtain appreciable mixing 
before ignition unless the fuel and oxidizer are thrown violently together. 
Ignition time is therefore not an important determinant of blast yield for hyper - 
golics, but impact velocity and degree of confinement after impact are impor- 
tant factors. Project PYRO results and resulting prediction methods which 
are given in references 15-17 concentrate on these factors and can be used 
directly to obtain estimates of blast yield. The only modification which we 
propose is to use smoothed curves from PYRO results for peak overpressures 
and impulses, rather than multiplying factors. 
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The procedure is then as follows: 

(1) Consider failure mode, or impact velocity and type of surface 
impacted . 

(2) Obtain terminal yield Y in % from Table III or Figure 8 (from 
reference 16). 

(3) Calculate W from Equation (2) knowing total weight of propel- 
lant and Y . 

(4) At distances R of interest, compute Hopkinson-scaled distance 
R/W 1/3 . 

(5) From appropriate smooth curves in Figures 9 and 10, r ® a<3 

overpressure P and scaled impulse I/'W ^ • Multiply scaled 

impulse by to obtain I . 


TABLE IV - ESTIMATE OF TERMINAL YIELD (REF. 16) 


FAILURE MODE TERMINAL YIELD RANGE (%) ESTIMATED UPPER LIMIT 


Diaphragm rupture (CBM) 
Spill (CBGS) 

Small explosive donor 
Large explosive donor 
Command de struct 
310-ft drop (CBGS) 


0.01 - 0.8 
0.02 - 0.3 
0.8 - 1.2 
3.4 - 3.7 

0.3 - 0.35 

~1 .5 


Note that the blast yields are very low (a fraction of one percent to a 
few percent' for all but high velocity impacts, which is not surprising in view 
of the small amount of mixing which is possible before ignition. Possible 
error in estimation of yield is also substantial, as can be seen from the 
ranges of yields in Table IV and data scatter in Figures 8, 9 and 10. 

2 Liquid O vypen-H vdrocarb on Propellant - The second propellant 
combination which we will consider uses Kerosene (RP-1) as a fuel, and liquid 
oxygen (LO z ) as the oxidizer in stoichiometric mass ratio of 1/2.25. Because 
this liquid propellant is not hypergolic, considerable mixing can occur in 
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various tvpes of real or simulated accidents, and time of ignition after onset of 
mWng is an important determinant of blast yield. Other important parameters 
have been shown by the PYRO and other results to be the mode of fa.lure or 
simulated accident, impact velocity, and propellant mass or *»elg • „ th _ 

important parameters appear to be geometry of tankage expressed 1 as a length 
to-diameter (L/D) ratio, propellant orientation and area of rupture of mte 
bulkhead for CBM case. The estimation methods which we gw 
based largely on the PYRO test results, but also include conclusions and/or 
physical reasoning of our own and of Farber and other investigators. 

For the case of milting and an explosion within the missile * a "^ago 
(CBM), time for ignition and mass of propellant are t e principa e 
of blast wave properties. The scaling of ignition time assumed for PYRO ’ 
not proven by the PYRO test results (see Appendix A), so we simp yf 
ffooth curve through PYRO results for blast yield Y as . .««« • 

F i„ ure 11 We also use Father's physical reasoning in p.otting this cur 
i e for zero time for mixing, yield must be zero, and for long enough time 
yield must decrease. A direct plot against ignition time is ^ed, mdepe 
of mass of propellant, because it fits the data as well as scaled time P l °‘= a " d 
also serves to indicate that scaling of ignition time has not yet been verified 

experimentally. Once blast yield Y has been determined from Figure l .or 
an assumed ignition time, effective weight of propellant W is then calculated 

from Equation (3) for known Y and W T , and blast “"LaXtl'T me 

are obtained from fits to PYRO data in Figures 12 and 13, m exactl, 

manner as for the hypergolic propellants. 

For simulated fall-back on the launch pad (CBGS), impact velocity as 
well as ignition time are important parameters in estimating yield 

PYRO prediction methods included fits to scaled time parameters and 
pact velocity to a fractional power close to one. As stated e ore, ime 
scaling is not proven by the data. Also, linear dependence on impact velocity 
L sampler than a fractional power close to one and fits the data just as welh 
A suitable fit of maximum yield Y^ to impact velocity, agreeing wi 
A of Figure 5-41 of reference 16, 


is : 


2. 08 


Y m 5 5 * + (F/S u i ' 




80 ft/sec 


(4) 


Blast data for this case 


where Y is expressed in %, and Uj is in ft/sec 

Lorn reference 16 aye normalized wii respec, «c Y m by the faclor # 

x. 100/Y and plotted versus ignition time in Figure 14. 

Leigh L’data can ,heu be used ,o predict percent yield using \ from 

Luafion ,4) to- impact eeloci.y U, . Again. "!"» V" which 

found and blast parameters determined from suitable fits to PYRO data, which 

are given here as Figures 15 and 16. 
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FIGURE 1 1. IGNITION TIME VS YIELD LO^RP- 1 CBM 




FIGURE 12. PRESSURE VS SCALED DISTANCE FOR 
LO /RP-i CBM CASE 
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SCALED POSITIVE IMPULSE I /W , psi-msec/lb 



FIGURE 13. SCALED POSITIVE IMPULSE VS SCALED DISTANCE 
FOR LO^/RP-l CBM CASE 
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FIGURE 14. IGNITION TIME VS NORMALIZED YIELD 

LO /RP-1 CBGS 
£* 


y 7 








In high velocity impacts of this propellant, the situation is scmew it 
simpler because there is little ignition delay and therefore only impact veloci- 
ty affects yield. Prediction methods from reference 16 can then be used o 
estimate yield Y (see Figure 17) and blast parameters obtained from Equa- 
tion (3) and Figures 15 and 16. 

3 . Liquid Oxygen-Liquid Hydrogen Propellant. - The final propellant 
combination is the entirely cryogenic combination of liquid hydrogen ( 2 

fuel and liquid oxygen <LC> 2 ) oxidizer in stoichiometric ratio by mass of /S. 

The rationale for predicting blast parameters for this propellant com- 
bination is identical to that for L0 2 /RP-1. For the CBM case Figure 18 
gives a plot of ignition time versus yield. After determining rom - 3 

plot and Equation (3), one can enter Figures 19 and 20 to obtain blast wave 

properties . 

For the CBGS case, a linear fit of maximum yield versus impact 
velocity gives 


Y =10% + 

m 


1 • ? 5 — u , 0<U< 80 ft/sec 

(ft/sec) I — I 


( 5 ) 


Using this equation to normalize ignition time versus yield, we obtain 
Figure 21. From the curve in Figure 21, we can find Y for a given ignition 
time, and then obtain W from W T and Equation (2). Finally, blast pressure 
and impulse are obtained from Figures 22 and 23. 

For high-velocity impact of this propellant, the blast yield is also 
dependent only on the impact velocity, and the prediction methods of refer- 
ence 16 can be used directly. The curve in Figure 24 gives the yield Y . and 
Equation (3) and Figures 22 and 23 will give predictions of blast wave proper- 

ties . 

4. Limit to Yield for Large Mass of Propellant - Any method of pre- 
dicting blast yield must, it seems clear, provide an upper limit on yield with 
increasing mass of propellant. All tests to date with large amounts of pro- 
pellant have shown autoignition sources which prevent mixing of even a large 
portion of the oropellant prior to ignition. The simplest way to incorporate a 
limit is to use the limit curve of Figure 7 generated by Farber. The curve 
labeled ' average value” should probably be used for our purposes. Explosive 
yield y in Figure 7 is only approximately related to Y in all of the preceding 
discussion, because Farber bases his yield on actual energy of combustion of 
the propellant rather than TNT equivalence. But, considering the multitude of 
other errors in these methods, we can igno.e the difference and simply assume 

that y = Y / 100. 
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FIGURE 17. TERMINAL YIELD VS IMPACT VELOCITY FOR 
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FIGURE 24. TERMINAL YIELD VS IMPACT VELOCITY FOR 

LO z /LH 2 (REF. 16) 


48 



III. DETERMINATION OF FRAGMENT 
VELOCITY DISTRIBUTION 

A. General 

In our review of accident and test data for liquid propellant explosions, 
we found no data on the velocities of fragments produced by these explosions, 
outside ofdata available in films of Project PYRO. Only the report by Jeffers 
(ref. 18) contained an analysis of fragment velocities, and this analysis was 
based on a study of films of three of the PYRO tests. 

The repository for all raw data, including f’lms, from the extensive 
series of PYRO tests is the Air Force Rocket Propulsion Laboratory at 
Edwards Air Force Base, California. In visits to this agency, SwRI staff 
members viewed all of the library of PYRO films and selected for data re- 
duction films of 94 tests in which individual fragments were visible and could 
be followed from frame to frame. These films were then loaned to the insti- 
tute for measurement of fragment velocities. For most of these experiments, 
several camera views of each test were available. Film speeds were accu 
rately known, with timing marks at known repetition rates impressed on tne 
edges of most films. The 94 tests for which we reduced fragment velocity 
data represented a spectrum of propellant types, scale of test, and type of 
simulated accident. In addition, reduced blast data and measured blast yields 
are known for each experiment and reported in references 15 and 16. 

Because the PYRO films constituted our raw data for determination of 
fragment velocity and because reduction of these data required more detailed 
knowledge of test details than were available in references 15 and 16, we en- 
gaged a consultant from URS Corporation, prime contractors for PYRO. He 
provided us with sufficient additional information, drawings of test arrange- 
ments, etc. , to enable us to obtain accurate estimates of such needed informa 
tion as distance scales in the field of view of each camera and timing mark 
frequencies on films. He also helped us resolve unexplained discrepancies in 
reported PYRO test results. 

Although the PYRO films provide all of our basic data on fragment 
velocities and can allow correlation with measured blast yields, they do not 
in general provide any data on fragment masses, shapes, or ranges. In only 
one PYRO test were these parameters measured, and this test is the only one 
for which any overall correlation of fragmentation effects conceivably can be 
made. 


B. Reduction of Film Data 

1 . Camera Locations - Locations of the camera during PVRO tests 
from which fragmentation data were obtained are shown in Figure 25. 
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Cameras were located along radial lines at azimuth angles 0 , 240 and 340 
and on a tower directly above the event.* A few tests were also photographed 
from airplanes and from the tops of nearby mountains, but these cameras 
were located too far away from the event to provide fragmentation data. 

During five tests cameras located at position B of Figure 25 were focused on 
the barricade located along an extension of the zero degree leg. Although 
fragment data were read from these films, they were not processed because 
the data were too few to be statistically significant. 


The greatest number of films and the largest amount of data came from 
cameras located at the azimuth angles of 0° and 240 . These camera positions, 
also called positions "A :1 and "B", were located 420 feet from the center of t e 
test pad. From these locations, cameras filmed the event with several di er- 
ent focal length lens and at indicated camera speeds ranging from 64 frames/ 
sec to 1,000 frames/sec. Although overhead shots were available for most 
tests, they provided little usable data. The field of view of the cameras was 
small relative to the size of the fireball and the fragments were not visible in 
the flames. Cameras situated at azimuth position 340 were located 1,050 eet 
from the center of the pad and were mounted on a tower 110 feet above the 
ground surface. Since the elevation at the base of the tower is 36 feet above 
that at the test site, the cameras were actually 146 feet above the test pad. 

Only a few tests were photographed from this location, however, and not all 
of the films provided fragmentation data. Some cameras had a field of view 
which was so large relative to the size of the fireball that fragments could not 
be seen. Other cameras with longer focal length lenses provided good film 
data. 


Generally, most data were obtained from those cameras which had a 
field of view slightly larger than the resulting fireball and which had the high- 
est film speed. Because several camera views were often available from the 
same location (particularly for locations A and B), the view or views providing 
the best data were selected for reduction. For example, other factors such 
as field of view and quality of the film being equal, the film which had the 
higher framing rate was chosen for data reduction. Even at the higher fram- 
ing rates, the film length was sufficient to observe the growth and start of 
decay of the fireball. 

Most films were accompanied by data which identified the test number, 
camera number, camera speed, frequency of the timing mark generator, 
focal length of the lens, distance from the camera to the event, and the azi- 
muth angle of the camera position. We found some of these data to be incon- 
sistent. For example, we were informed by Project PYRO test personnel 


* Obtained from drawings of field test arrangements for Project PYRO. 
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that cameras located 420 feet from the event were located at position A or 
position "B" (of Figure 25) on radial lines which are 120 aparh We foun , 
however, that these positions were either labeled as 0 and 270 , respective- 
ly or as 325° and 130°, respectively. Also, leg A was sometimes identified 
by an azimuth angle of 65° and the tower (located 1, 050 feet from the pad) by 
azimuth angles of 300° and 340° (340° is correct). Consequently, we relied 
upon the background in the film itself to identify the leg upon which the camera 
was located. Except for cameras located directly above the event, (and those 
at distant locations which did not contribute fragmentation data), the carneras 
were assumed to be located at either position "A", "B", or at an azimuth 
angle of 340° as shown in Figure 25. This is consistent with the instructions 

of Project PYRO personnel. 


Film data reduction was accomplished using a Vanguard Film Analyzer. 
Data obtained from the films included the frame number relative to the initia- 
tion of the explosion, the X and Y positions of the fragment referenced to the 
frame number, the spacing of the timing marks, the height to the top and 
bottom of the tank, the tank diameter, and, if in the field of view, the heig 
of the tower above the test pad. The height of the tower was used for com- 
puting the scale factor. Since each test was viewed from more than one direc- 
tion (from camera locations A and B) an attempt was made to identify the 
fragments in camera views from both locations. When available, such data 
(which we have labeled "Duplicate View Data") permitted a fairly accurate 
determination of the flight path of the fragments. If the fragments were iden- 
tified only in a single view, however, (this has been labeled "Single View 
Data") it was not possible to determine their actual flight path. As recorded 
on film, the fragment position is a projection on a plane normal to the lens 
axis of the camera. Without additional data, such as from an additional 
camera view, no correction to the data to account for a flight path other than 
normal to the lens axis of the camera could be made. Thus, the fragment 
data were processed a though the fragments were traveling in a plane which 
was normal to the len , axis of the camera. After reducing the fragment data 
from the films we estimated that nearly all of the fragments sighted had a 
trajectory which was within ±45 degrees of the normal to the lens axis of the 
camera. If this estimate is correct, computed velocities for single view data 
should be within f0% to -30% of the true velocity. 


Since cameras at positions A and B viewed the event from approxi- 
mately the same elevation as the test pad. fragment positions were recorded 
for a vertical plane normal to the azimuthal position of the camera. How- 
ever, the camera on the tower, at azimuth 340°, was elevated above the test 
pad. A correction for this elevation was made so that the fragment position 
would be calculated for a vertical plane as for camera positions A and B. 
This correction is discussed in the processing for single view data. 
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2. Processing of Single View Data - As viewed from the camera, the 
fireball and fragment would appear as in Figure 26 with the positive X axis to 
the right and lying along ground surface and positive Y axis vertical throug 1 
the center of the tank. Subscripts on the axes refer to the camera position 
A or B from which the data are obtained. The view shown in Figure 26 is 
typical of that obtained by a camera in position A. The angle between the 
: nter section of the X-Y planes of cameras A and B is then 120 (this angle is 
required for processing of the duplicate view data, but not for processing the 
single view data). Fragment positions are measured on the film with re- 
spect to the X and Y axes. Processed data are computed with respect to the 
X ' and Y axes; that is, with respect to a set of axes which pass through the 
center of the tank. 

The first items calculated from the film data are the film speed in 
frames per second (fps), the scale factor (SF), and the height of the event 
(h e ). Spacing of the timing marks (tm) are read from the film in frames per 
timing mark. When spacing of timing marks is multiplied by the frequency 
of the timing mark generator in timing marks per second, the framing rate 
in frames per second is obtained. To compute the scale factor a known 
length in the field of view is divided by its length as measured on the screen 
of the Vanguard Analyzer. This give? the scale factor directly if the known 
length is the same distance from the camera as the event; otherwise, an 
appropriate adjustment must be made. The height of the event (h e ) can be 
determined from engineering drawings if they are available, but it is also 
determined directly from the film data. It is computed as one-half the sum 
of the distances to the top and to the bottom of the tank. Even though an 
attempt was always made to align the X axis (Y = 0 position on the screen of 
the Vanguard Analyzer) with the top of the pad, this alignment was never 
exact. Consequently, whenever it is available, h e obtained from the film 
measurement is used. Knowing the film speed, the scale factor and the 
height of the event, the X-Y positions of the fragment can be calculated from 
Equations (7) and (8). 


x. j - X j (SF) 
1 i 


(7) 



(SF) 


cos 


*V 1 + ~ 


■fi (SF) 


dc 


sin 


>v_ 


( 8 ) 


Subscript i refers to the frame number and superscript j refers to the frag- 
ment number. Lower case x and y denote the values as read from the film, 
and capital letters denote calculated displacements. The denominator of 
Equation (8) is a correction factor which adjusts for the condition that the lens 
optical axis may not be perpendicular to the vertical axis along which Y is 
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to be calculated. For the camera elevated above the event as shown in 
Figure 27, 9 V is the positive angle between the lens axis of the camera 
the horizon and dc is the horizontal distance from the camera to the event. 

To obtain the vertical position of the fragment relative to the center of the 
tank, h p , the distance from the X-axis (Y = 0) to the center of the tank is 
subtracted from the calculated displacements. Equation (8) is deri J ed cr 
the condition that the lens axis passes through the center of the tank (the 
assumed origin of the fragment trajectories). The error induced by neg- 
lecting the offset of the lens axis relative to the center of the tank is about 1%. 

A correction similar to that in Equation (8) is needed when the frag- 
ment is not traveling in the X-Y plane, but is traveling either toward or away 
from the camera. The error incurred is indicated in Figure 27. h is exag- 
gerated because the camera is shown much closer to the event relative to e 
height of the fragment than actually was the case when photographing the 
fragments. This schematic does show, however, that the calculated Y dis- 
placement is a projection of the actual displacement on a vertical plane 
through the center of the event. If the fragment is coming toward the camera 
the Y displacement calculated will be slightly greater than the actual \ dis- 
placement; whereas, when the fragment is traveling away from the camera, 
the calculated displacements will be slightly smaller than the actual displace- 
ments. For single view data, there is no way to correct for this error e 
cause one cannot tell which direction the fragment is actually traveling. (A 
similar error occurs for the horizontal displacements. ) For the duplicate 
view data, however, it is possible to calculate the true vertical and horizontal 
displacements. This will be discussed under "Processing of Duplicate View 

Data. " 

Once the X and Y displacements of the fragments have been calculated 
for each frame number and once the times corresponding to these frame 
numbers have been determined by dividing the frame number (event starts in 
frame number one) by the framing rate, the fragment velocities can be calcu- 
lated. Velocities are calculated in two ways. An average velocity is com- 
puted by assuming that the fragment has traveled in a straight line from the 
center of the tank to the position where it is first sighted. This velocity is 
given by Equation (9) 


j V( x ii) + (Wi) 
T A 


(9) 


where the subscript II indicates the frame number in which the fragment is 
first observed. In addition to this so-called average velocity, the velocity of 
the fragment during the period of time in which it is observed is also com- 
puted. To determine this velocity, the distance the fragment travels along 
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VERTICAL DISPLACEMENTS 
CALCULATED FROM EQUATION (3) 



FIGURE 27. SCHEMATIC SHOWING ACTUAL AND C ALC ULATED 

FRAGMENT VERTICAL DISPLACEMENT 








its trajectory after first being sighted is calculated from Equations 
and (11) 


( 10 ) 



x Li ) 2 + R 


Y>. 


i-1 


( 10 ) 


s! = S? , + AS^' 


(ID 


"i-l 


where the distance S is set to 0 at the time of sighting. . 

ment U-, is then determined from the slope of a straight line fit to the co 

• • Si and T 


Velocity of the frag- 
^ht line fit to the co: 

puted distance -time data points and T i# where Tj is adjusted to 0 at 

the time of sighting as was The subscript i on the velocity is used to 

denote an instantaneous velocity. While it is true that this veloc.ty is actually 
an average over the time of sighting, the elapsed time is short compare 
the elapsed time between the explosion and the sighting o t e ragmen 
which U A was determined and the velocity of the fragment varies only slight- 
7y during**. time ot sighting. Also, the position of the 

spends to is computed for the center of range over which the fragm 

tracked (as discussed below) and Uj should be very close to the mstantaneo 

velocity at this position. 

In addition to the fragment velocities, the position of the fragment 
relative to the center of the event and its direction of travel »” e j 

Fragment coordinates are given as a radius R, height Z. and an azimuth 
angle Y. Since the computed velocity. Ui, is taken to be the instantaneous 

velocity of the fragment at the center of its range of travel while in -tew t e 
cylindrical coordinates are computed for this position also. The expressions 
for determining these coordinates are 


R J 


Z J 


IF 


= 1 

( X IF + X Il) 

(12) 


1 .J v. j \ 

(13) 

= 4 

^ y if + y ii) 

v j 


Xjp 

/nn\ , W 

(14) 


WfI 


Y j 

<0.0 set Y J = 360 + 

(15) 


where the subscript IF denotes the final frame in which the position of the 

fragment was recorded. Notice that R and Z are simply the center posi- 
tion of the : mge of the X and Y coordinates, respectively. The azimuth 
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position i 3 determined from the azimuth angle, * c , of the camera viewing 
rte fragment and the sign of the X displacement. This is based on the assump- 
tion that the fragment is traveling in the X-Y plane which is normal to the line 
of sight of the camera. (For fragments which are photographed andproperly 
identified from two different camera positions, the azimuth ang e o 
of flight can be determined more exactly. ) The last parameter compu 
the fragment is the elevation angle of its flight path above the horizon. This 
angle, 9, as shown in Figure 26, is determined from the slope of a straight 
line fit to the X- Y data for each fragment. 

3. Processing of Duplicate View Data - After the single view data for 
each fragment had been processed, data for those fragments which could be 
identified in more than one camera view, that is, in views from camera posi- 
tions A and B of Figure 25, were combined to more accurately determine a 
true trajectory of the fragments. Fragments were difficult to identify m more 
than one view, however, and thus the duplicate view data are very limited. 

In addition, we have often found when processing the data that what appeared 
to be the same fragment in two separate views, apparently was not. One cri- 
terion for determining whether or not the fragment is the same in two views 
are the Y displacements. They should be approximately the same when view- 
ed from the two camera positions. As noted in Figure 27, there wil e some 
differences depending on whether the fragment is traveling toward or away 
from the camera; however, these differences should be small because the 
cameras are located 420 feet from the event and it is unusual to follow a frag- 
ment for a distance of more than 60 feet above the ground. Thus, the cameras 
are viewing the fragment at an angle of less than 10 above the horizon. Al- 
though no hard and fast rule was established, when the Y displacements of the 
two fragments differed by more than about 15%, it was assumed that two dif- 
ferent fragments rather than the same fragment had been observed. 


In duplicate view data, the same fragment is observed from cameras 
situated 120° apart. Thus, errors introduced by the fragment traveling 
toward or away from the camera tend to cancel out when averaging the Y dis- 
placements computed from the two different camera views. Consequently, 
the vertical displacements, 6, are calculated as the average of the vertical 
displacements from the single view data. This is shown in Equation (16) 

*; ■ t («; • »“") 

where the superscript DUP is the number of the fragment which provides 
the additional camera view of fragment "j". 


Calculating the horizontal position of the fragment is more complicated. 
Figure 28 shows the geometry in plan view. The X displacements of the frag- 
ment, computed from the two separate camera views, are indicated by X^. 
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and Xo, These are projections of the fragment horizontal displacement upon 
the X- Y planes for each camera position. The true horizontal displacement of 
the fragment, Hj. relative to the center of the event is desired. Using 
simple geometrical relationships, two expressions for H F were obtained ana 
are given by Equations (17) and (18) 


H 


X A d ^A 


F ~ X A cos a + dc A sin a 


for X A ± 0 


H 


- X B dc B 


F = Xg cos(^0" + a) + dc B sin (60 + a) 


for 


B 


* 0 


(17) 


(18) 


As indicated. Equation (17) cannot be used for X A = 0 and Equation (18) 
cannot be used for Xo = 0. In practice. Equation (17) was used whenever the 
absolute value of X A was greater than X B , and Equation (18) was used when 
the reve. se was true. Subscript i and superscript j have been dropped 
from the equations for simplicity. The angle, a, is found y equa mg qua 
(17) to (18), and is given by Equation (19). 


, X a X r cos 60°+ X a X b dc B /dc A 
tan* * : _ . : : — , ^ o „ T“ v 


(19) 


Equitions (17), (18), and (19) hold for fragment positions in all quadrants. 

The same procedures were used in calculating the average velocity, 
U A , and the velocity of the fragment after acquisition, Uj, as were used 
for the single view data. Quantities U and H were simply substituted for 
X and Y. Also, with the exception of the azimuth angle, the cylindrical 
coorc.inates of the fragment and the elevation angle 8 were computed as for 
the single view data. The azimuth position of the fragment is simply 360 
minus the angle 1 . 


4. Computer Program for Data Processing - A computer program was 
written to process the data as it was received in punched card torm from the 
Vanguc rd Analyzer. Although the details of the program will not be presented 
here, it solves the ec -a >ns given in the preceding sections and prints the 
results in a usable form. The program also creates a tape file for subsequent 
statistical analysis of the results, 

A summary of camera data and number of fragments whose trajectories 
were measured is given in Appendix A. The tests are grouped by type, as 
defined in the PYRO reports (refs. 15 through 17). and are identified by the 
PYRO test numbers. Numbers of fragments whose velocities and trajectories 
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were measured range from a minimum of 2 per test to a maximum of 36 per 
test. The maximum number of fragments identified in duplicate views was 

3 per test. 

5 Summary of Fragment Velocity Measureme nts and Calculations - 
The fragment velocity data derived from the PYRO films are available on 
punched cards and in computer printouts. The key designator for all data is 
the PYRO test number; for each such test having readable fragment ve ocity 
data, all of the quantities defined by the equations in Section 4 are calculated 
and printed for each fragment traced. The printouts are quite voluminous 
and so are not included in this report. 

The data for each test were further reduced using standard statistical 
procedures and a " anned" computer program available to us through our CDC 
library of programs. We chose the six output variables U A , Up 9, R, 

and Hp for statistical treatment, and computed the mean, standard devia- 
tion, standard error of the mean, maximum, minimum, and range for each 
variable. Again, the computer printout sheets are too detailed for inclusion 
in the report, but are available for study. 

A summary of results is given in Table V. Each test is identified by 
the PYRO test number, and the tests are grouped by propellant and type of 
simulated accident. Test conditions given are total weight of propellant, 
tank L/D ratio, and impact or drop velocity. Ignition time and measured 
blast yield from reference 16 are listed.' Fragmentation data are given in 
the last four columns, indicating number of camera views from which data 
were taken, number of fragments observed, mean value of fragment velocity, 
U- , and standard deviation of this same parameter o u . These mean velocities 
and standard deviations were determined from the single view data only. The 
number of fragments which were identified in more than one camera view and 
processed as duplicate view data was too small for statistical analysis. Be- 
cause the data represent a wide spectrum of test conditions, propellant types, 
and propellant weights, they should allow correlation with methods of pre- 
diction of initial velocity such as those presented later in this section, and 
snould also lend themselves to various types of statistical analysis. A limi- 
tation to the data which may render statistical studies difficult is that rela 
tively few fragments could be observed in any one test. Some grouping be- 
tween tests with different propellants, different blast yields, etc. . may be 
possible. To aid in rational choice of ways to group the test results, we con- 
ducted several limited model studies. These are included in this report as 
Appendix C. 


* Blast yield in percent is defined in Section II. 
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c. Methods of Predicting Velocity 


1 # Explosions IVithin Missile Tanks 


a. 


A Deterministic Method 


(1) Choice of Par ameter s and Idealization of the Problem - 
In this section, we discuss a method for predicting maximum velocities for 
fragments generated by explosions within missile tankage (CBM case). A 
method for predicting ranges is discussed in Section IV. In order to predict 
maximum fragment velocities, the problem must be "idealized" in a number 
of ways, and the parameters must be chosen such that the theory may be 
applicable to a spectrum of missile explosions involving different liquid tuels 
and oxidizers, volume of fuel oxidizer mixtures, tank wall specific mass, and 
tank geometries. To idealize the problem we consider only a spherical 


volume V OQ of radius 


R in which we have a stoichiometric mixture of that 
° - ’ • — — Further- 


amount of fuel and oxidizer which is mixed at time of detonation, 
more, upon detonation, all of this volume of fuel and oxidizer is assumed to ^ 
be converted instantly to explosive products and energy, forming a "hot gas 
which can be characterized by some ratio of specific heats. *, initial pres- 
sure, P OQ , and initial sound velocity a OQ . All of the fuel or oxidizer ex- 
ternal to this spherical volume is considered to add mass to tne container 
wall only and takes no part in the chemical reaction of the explosive process. 
The container is spherically coincentric to the sphere in which the explosion 
takes place, and the non- reacting fuel (or oxidizer) moves initially with the 
fragments of this container wall as if they were an integral part of the frag- 
ments (i. e. , they just add mass to these fragments). Compression of the 
liquid fuel and oxidizer, and shock transmission effects in these liquids are 

ignored. 

All fragments are assumed of equal size and have circular projections. 
If we pick a parameter, n, for the number of such fragments into which the 
containing sphere fractures, we can later solve the problem for various n, 
thus effectively varying the fragment size. It is further assumed that the sum 
of the concave inner surface areas of all the fragments is equal to the surface 
area of the original sphere of explosion. 


Ai = 4 


R 2 /n 


( 20 ) 


where R is the initial radius of the sphere of explosion [R - r i* T j for t 0 ], 
and A; is the fragment surface area. Although these are gross idealiza- 
tions of the fragmentation picture, they make the problem amenable to solu- 
tion and allow some estimate to be made of how sensitive the solution is to 
fragment size and number. 
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If the stoichiometric mixture of fuel and oxidizer occupies a volume 
V and has a mass M Q , then, at the time t = 0 the instant after the ex- 
plosion, we assume the gaseous products of explosion occupy the same 
volume and have the same mass. The sound speed in this medium at this 
time is then a function of K, one of the independent parameters. One must 
pick k, M , V (0), and P & (0) = P QO which characterize the sphere of 
explosion at t = 0 (x, M Q , and P QO will depend on the type of fuel and oxi- 
dizer, and V depends on the amount of fuel and oxidizer that have mixe-/. 
One must then also pick an R, n, and Mj. which characterize the con- 
taining sphere (where R is the containing sphere's initial mternal radius, 
equal to the external radius of the sphere of explosion; n is the number of 
fragments into which the containing sphere fragments, and M t is the total 
mass of the containing sphere, i. e. , the mass of the non-reactants, as well 
as the actual fuel tank mass). With this information the maximum fragment 
velocity for any of the equal sized n fragments is calculated ' / the method 
to be described. This maximum velocity, U m , maybe used as the initial 
fragment velocity for the fragment range calculation (Section IV) m which the 
fragment is assumed to experience no acceleration due to the explosion after 
U m is attained, but it experiences deceleration due to drag forces in the 
medium through which it is traveling. 


(2) Equations for the Mathematical Model '' - To obtain 
the initial velocities for the fragments from CBM explosions, an extension 
was made of the techniques of D. E. Taylor and C. F . Price (ref, 31) and of 
G. L. Grodzonski and F. A. Kukanov (ref. 32) relating to the motion of frag- 
ments from bursting gas reservoirs in a vacuum. These techniques were 
generalized from the case of spherical vessels which broke up into two per- 
fectly hemispherical fragments to the case of a spherical vessel which frag- 
ments into n fragments, each of equal size and having a circular projection. 
For the case where the containing sphere is thin-walled, the equation of 
motion for the i-th fragment of mass and radial displacement r^(T) 

having a projected area F and experiencing a pressure Pj(t) at time T, is 
given by: 


d 2 ri (t) 

Mi - FPj(t) 

dt2 


( 21 ) 


If Ai is the area of the curved surface of the spherical segment (for the i-th 
fragment), then 


The equations given in this section require a number of mathematical opera- 
tions for their development. Many of these operations are omitted for 
brevity. All nomenclature is given in Appendix I. 
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( 22 ) 


Aj = 2tt Rh = 4tt R 2 /n 

where R = r^O) and h is the segment height. Let the segment radius equal 
d/2, then it can be shown with the use of Equation (22) that 



or 


F ~ Aj 


for n 2 >> n 


(23a) 


(23b) 


The equation of state assumed for the gaseous products of explosion is given by 


P q (t) V o< T > = C < T) R 1 T o< T > 

The rate of change in the mass of the confined gaseous explosion products in 
the sphere as the sphere begins to fragment and gas is lost to the externa 
sphere region is given by 


dc (T) 
dT 


-k P„, a „ H w 

'!-• '■I s 


(25) 


where p„ and a„ are the critical gas density and velocity as they escape 
through the cracked surface of the sphere, n is the crack perimeter, w is 
the crack width, and k is a discharge coefficient. The crack area for t e 

i-th fragment given by (II w)i is equal to the difference between the area sub 

tended by the initial solid angle, 0, at the distance r^T) for time t and 
the actual area of the fragment A i# 


(Hwlj = 0 r- 2 (T ) 


4ttR 2 f r _£±2l 

n \ R 2 


(26) 


Inherent in Equation (26) is the assumption that all fragments travel in a radial 
direction and that the symmetry of the fragment motion is such that the equa- 
tions of motion are only a function of the magnitude of the radius, i. e. , the 
motion of the i-th fragment in the radial direction describes the motion of 
any other fragment in the radial direction and shear forces are small com- 
pared to dynamic forces. Thus, we may set x i = r and, from Equations (25) 
and (26), we have 


dc (r) 
dT 


-k p 


••• (T - 0 


4tt R* 


(27) 
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The volume of the gaseous explosion products at any time T after 
T = 0 is given by 

V q (t) =( 4/ 3)nr 3 ( T ) l28) 

It has been assumed that the confined gas immediately adjacent to the 
fragments is accelerated to the velocity of tne fragments but that this con- 


stitutes a negligible fraction of the gas, the great bulk of which is unaccel- 
erated. Thus, from one -dimensional flow equations 

h/(k- 1) 


Pi(T) 

= P 0 < T > j^ 1 ' 

/ H- 1 

v 2(a 0 (T )) 

P... = 

p o' T )(„ +1 ) 

. 1/(M-1) 

1 

a.,. = 

( 2 ^ 
a o (T »VK + U 

I 1 " 


(29a) 


(29b) 


(29c) 


Equations (29) can be non-dimensionalized by setting 

r(T) = Xg(5), T=es, P 0 (t> - P 00 p *<5> 


(30) 


For a confined gas that behaves adiabatically we have: 


folU 

Pqo 


Pq(t) \ K /T q (t) \ h/(K ~ 1} 

Poo / V T oo / 



(31) 


where the double zero subscripts refer to the confined gas condition at T 0 
From Equations (20), (29). (30), and (31), one obtains: 


Mi „ 
8 


(g')‘ 


(P 


a) (h-1)/k J 


x/(K-l) 


( 32 ) 


where primes denote derivatives with respect to Mi is the mass of 

the i-th fragment, and M t is the total -shell" mass. For n fragments of 

equal size. 


M i _ l (33) 

Mt n 
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thus. 


p r, < g'> 2 1 

! * L (p*)'*- 1 ’'* J 


h/(x- 1) 


(34) 


where 


X = 


f; 

F P 


oo 


oo 


(A)- 


Q = 


Mta 


oo 


F P, 


oo 


(A) 


1/2 


Differentiating Equation (31) and nondimensionalizing, one obtains 

(h-1)/2x 


K 


■a g“ + ag 


„ 2 / 
■3k g g 


(35) 


(36) 


where 


3(n) = R' 


■ n 2/_L.^ 
n V H- 1 / 


P oo 2 


— 2 4 
Mt a oo 


n 


and 


a = 



(k+1)/2(k- 1) 



- 1/2 


Equations (35) and (36) must be solved simultaneously with the initial 
conditions 


r(o) = R 


g(o) = R/X 


= o 

dT 


d gM_ o 

dS ' 


(37) 


P«(o) = P. 


OO 


P*(o) = 1 


for values of § until the fragment acceleration is small with respect to the 
acceleration -t T = 0. dg/d? is a maximum as the acceleration goes to 
zero. dg/dS is the normalized initial velocity for the fragments to be used 
in the calculation of the fragment ranges. 
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( 3 ) Solution of the Equations by Numerical Techniques 
and Co mputer Program - The solution of Equa- 
tions (34) and (36) may be obtained numerically for the initial conditions of 
Equations (37), using the Runge-Kutta method. A program which does this 
(in FORTRAN IV), as well as the definition of its symbols, is given in 
Appendix C. This program requires as input the characteristics of the gas- 
eous explosion products at T = 0 (A0, the speed of sound; P0, the initial 
pressure; and CAP1, the ratio of specific heats), and the characteristics of 
the vessel (RR, the internal ra.dius of vessel and unburned fuel; TM, the mass 
of the vessel and unburned fuel; and FN, the number of fragments). Also, a 
guess must be made on the elapsed time between detonation and zero fragment 
acceleration, XMAX.* A time interval spacing, AH, must also be chosen 
for the calculation. 

The program outputs the normalized displacement, velocity and accel- 
eration (dynamic variables) of the fragment as a function of time; the normal- 
ized pressure within the vessel as a function of time; and the fin^l values of 
time, displacement, velocity, acceleration, and pressure at XMAX. A 
sample program run is given in Appendix C. 

Figure 29 shows the results for computing the maximum fragment 
velocity as a function of number of fragments when all other parameters are 
held constant. The figure indicates that the number of fragments does not 
affect maximum fragment velocity except at very low a (at which point the 
assumption of Equation (2 3b) does not hold). Figure 30 shows how maximum 
fragment velocity varies with K and the mass ratio M 0 /Mj., The ratio of 
explosive products mass to shell and non-reactant fuel mass, a constant 
initial pres sure (P QO ), number of fragments (n), sound velocity (a Q ) f aid 
radius of explosion at T = 0 (R) are assumed in this figure. 

(4) Comparison of Results from this Method with 
Other Sources - The results of this method are 
compared to the cases described by Taylor and Price (ref. 31) in which: 

P,V 

a = 9 -g- = 2.55, 0. 1436 

M t a 0 2 


A "guess" is made of XMAX because the computer has to work with some 
finite time. If XMAX is guessed too short, then g "( l) (the normalized 
acceleration) is not reduced close enough to zero to yield an acceptable 
estimate of g'd). If XMAX is guessed too long then the computer can not 
calculate the final values of g" and g'. In practice we watched g ' for 
each program iteration and accepted a g' which remained constant to three 
significant figures between iterations as the value of g for g ' nearly zero. 
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AGMENT VELOCITY AS A FUNCTION OF 







10.0 in 

10 2 FRAGMENTS 



\GMENT VELOCITY 





in Figure 31. Generally it can be seen that our results preset somewhat 
lower fragment velocities than they do. Some discrepancy is to be expected 
because our assumptions on geometry were not as precise as theirs. Thxs 
is especially true where we assumed that the surface area of a spherical 
volume could be divided into n equal circular areas; they assumed only ^two 
hemispherical fragments whose projected areas were well defined. None- 
theless. the agreement is relatively good, especially for the greater a. 

Table VI gives a comparison between fragment velocities measured 
in experimental work relating to hazards from bursting high pressure tanks 
(ref. 33) and the predicted velocities using our deterministic method. 1 e 
experimental values of Pittman (ref. 33) were obtained by pressurizing 
spherical metal tanks with N 2 until they burst. Fragment velocities were 
measured by use of a breakwire system. In this system, a time inters a 
counter was initiated by the tank rupture (start pulse) and was stopped w en 
the first fragment reached the breakwire l foot away. Thus, the measured 
values actually are the mean fragment velocity for the first foot of travel 
(a stroboscopic photographic technique was also used to measure fragmen 
velocities in some cases). Where the experimental values were not prec sel> 
determined in tests D and E, limits were assigned to the fragment velocity on 
the basis of the data obtained. Input data to our progi am was based on the 
tank geometry and burst pressure, described in tn« report, and the proper les 
of N 2 . In general, our values agreed well with the measured values. 

We believe on the basis of these comparisons with independent data 
that our method ha. proven to be sufficiently accurate in predicling fragment 
velocities to be useful. An effort to predict data reduced from PYRO films 
using this method is described later in this section. 


TABLE VI. COMPARISON OF PREDICTED FRAGMENT 
VELOCITIES WITH PITTMAN'S DATA (ref. 3 3) 


Test 

— 1 S - ; 

( psi) 

Rad. 

(in.) 

Mass 

(lb-sec^/in.) • 

Measured 
(ft/ sec) 

Predicted Ilf | 

(ft/ 8 C C 1 j 

c 

I 

8 x 10 3 

r 

L 

-o 

rw L 

i 

_i 

I.63x 10“ 2 

1 .2 x 1 0 3 

, 1 

1.201x10’ | 

i 

1 

D 

8x 10 3 

27.0 

4.45 X 10" 1 

< 1.3 x 10 3 

1.61 x 10 3 

| 

E 

8.1 3 x 10 3 : 27.0 

4.4 3 x 10" 

> 1.17 x 10 3 

1 

i l.bO x 10 3 

L - 


7 5 


VELOCITY, in 'sec x 10 


15.0 



DISPLACEMENT, inches 


FIGi: HE 31. COMPARISON OF TAYLOR- PRICE TO SwKl 
SOLUTIONS FOR ADIABATIC CASE 


» 

i f 
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2. Appurtenances Subjected to Propellant Blasts 

a. Gener ' - "He situation considered here is shown in Figure 
32. A propellant explosion occurs after some accident which ruptures the 
tankage and causes propellants and oxidizers to spill, mix and ignite on the 
launch pad. We wish to establish a method of predicting velocities o w ic 
nearby objects (which we will call appurtenances) are accelerated by the 
passage of the blast wave. These objects can be parts of the launch l tow * ' 
storage tanks, vehicles, and objects in or attached to the upper stages of 
launch vehicle itself. 



FIGURE 32. SCHEMATIC FOR ACCELERATION OF APPURTENANCES 

BY PROPELLANT I? LAST 


b. Blast from Propellant Explosions - T t is clear from the re- 
sults of Project PYRO and other studies that the characteristics of blast 
waves produced by liquid propellant explosions differ significantly from the 
characteristics of blaot wave, produced by TNT or other conventional solid 
explosives. These difference, are discussed elsewhere (see Section . U but 
some will be reiterated heVe because they affect estimation of velocities > 
appurtenan 


The first characteristic of blast from propellant explosions which 
differs from that of TNT is that the blast wave le relatively weaker at all dis- 
tance. because the propellant and oxldi.er are almost never Intimately mixed 
before being ignited. Thus, the full explosive potential of the propellant- 
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oxidizer computed on the basis of energy of combustion is never realized. 

This difference is expressed, when comparing with blast waves from T 
explosions, as a percentage or fraction of "TNT yield" where the per- 
centage must be multiplied by the total weight of propellant plus oxidizer 
and then expresses the weight of TNT which will produce the same blast pro- 
perties far enough from the explosion source. 

The second characteristic of the blast from propellant explosions which 
differs from the blast from TNT explosions is that measured parameters such 
as peak overpressure and positive phase impulse vary with range in a different 
manner. Compared to TNT blasts, equal energy propellant explosions gener- 
ate blast waves with relatively low overpressures and relatively hig im- 
pulses (longer durations) at close-in distances, with gradual change to near y 
identical pressures and impulses with increasing range. Below about 10 psi 
overpressure, blast characteristics from the two types of sources appear to 
be essentially identical. The term "terminal yield" has been introduced to 
designate blast yield calculated from measured blast parameters in the low- 
pressure regime. The basis for estimating the terminal yield is usually a 
source of compiled blast data for hemispherical TNT charges detonated on 
the ground (see ref. 2). 


A third facet of difference between TNT explosions and liquid propellant 
explosions can seriously affect our ability to predict fragment velocities of 
appurtenances. For TNT explosions, there are available both a quantity of 
measurements of time histories of dynamic pressures, and also proven com- 
puter programs which can predict time histories of this and other physical 
properties of blast waves. For liquid propellant explosion®, the only com- 
plete body of data available is for time histories of overpressure, from Pro- 
ject PYRO (refs. 15 through 17). Methods are presented in Section !I for esti- 
mation of certain blast wave properties for several propellant combinations tn 
test configurations which simulate accidental spill and subsequent ignition 
(identified in Section II as CBGS - Confined By the Ground Surface). The 
parameters which can be estimated, given type of propellant, total weight of 
propellant W t , impact velocity Uj, time of ignition t, and distance R from 
the center of the explosion to an appurtenance, are the side-on peak overpres- 
sure P and the side-on impulse I. As we see later, these values must be 
used together with other data to estimate the pressures exerted on a body as 
the blast wave diffract® around it. 


c. Interaction of Blast Waves with Appurtenances - To be able 
to predict velocities to which appurtenances are accelerated by propellant 
blasts, we must know something about the physics of interaction of blast waves 
with solid objects. Some of this knowledge is briefly reviewed here. Figure 
3) show, schematically, in three stages, the interaction of a blast wave with 
an irregular object. As the wave strikes the object, a portion Is reflected 
from the front face, and the remainder dlffracte around the object. In the 
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FIGURE 33. 


INTERACTION OF BLAST WAVE WITH 
IRREGULAR OBJECT. 


diffraction proc..., the incident wav. front clo... in b *^ th * 

r; After paa.ag. of the front, the bod, i. inom.r.ed in a «im.-var,.n 8 . f 
flow field. Maximum preesure on the front face during th g P 

loading is the stagnation pressure. 

We are interested in the net transverse pressure on the ^ object as a 

function of time. This loading, somewhat idealised, is shown in Figu . 

(Details on calculation are given in The Effect, of ^f/y^^^Vrom /ero 
At time of arrival t a . the net transverse pressure rise, linearly from ce 

to maximum of P r in time (Tj-t*) (for a Hat-faced object, this time is 

zero). Pressure then fall, linearly to drag pressure in time (T 2 - T,). and 

then decay, more slowly to sero in time (T 3 - T 2 ). The time history o( drag 

pressure is a modified exponential, with maximum given by 


Cn Q 


t 0 , 


( 3 8 > 


where C n is the steady- state drag coefficient for the object Q i« P«-ik 

dynamic pressure, and 0. and u t are peak density and particle velocity re- 
spectively for the blast wave. The characteristics of the diffraction phase of 
the loading can be determined easily if the peak side-on overpressure or 

the shoe!; v-iocity U are known, together with the ehape and some charac- 
teristic dlmensio.. D of the object. The peak amplitude of the drag phase. 

CpQ. can also be determined explicitly from P § or u # . But. the time 




Pressure 



FIGURE 34. TIME HISTORY OF NET TRANSVERSE PRESSURE ON 
OBJECT DURING PASSAGE OF A BLAST WAVE 


history of the ensuing drag loading. C D q<t>. i- quite difficult to predict 
accurately for propellant blasts. 


d. Metho d of Eatlmatlng N et T.amm Preature - The 
method w. praaent*"herT"utlllxea the fit. of Section tl to PYRO data for e.de-on 
blaet parameter!, an a.eumed time hiatory of drag pr...ur. known to be r.a- 
aonably accurate for TNT and nuclear blaeta. e.timate. of diffraction Jimt. 
baaed on .hock tub. experiment., drag coefficient, from 
reflected and etagnatlon blaat front propertlea baaed on equation, which 

well known in blast physics. 


Side-on overpressure is often expressed as a function of time by tne 
modified Friedlander equation (see Chapter 1 of ref. 1). 


p (t) 


P(1 - t/T) e 


-bt/T 


(39) 


where T is the duration of the positive phase of the blast wave. Integrating 
this equation gives the impulse 


T 

f PT f, (l-«* b ) 1 (40 

I = j P(t» dt = — L 1 ’ b J 

o 

The dlmen.lonlea. parameter b la called the time conatant, la a function of 
.hock atrength. and la reported In Chapter 6 of ref. I. It la plotted graph- 
teally in Figure 35 for a range of shock strengths P. where 
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BLAST WAVE TIME CONSTANT b VS. DIMENSIONLESS 
OVERPRESSURE T~ 







p = P/Po 


( 41 ) 


and p n is ambient air pressure. The peak reflected overpressure P r and 
peak dynamic pressure Q are unique functions of P for a given ambient 
pressure p Q . For shocks of intermediate to weak strengths, P < 3.5, these 

functions are (see ref. 1, Chapter 6): 


2P 



and 


Q 


5 P 2 
2 7 t P 


where 


p r = P r /p 0 . Q = Q/Po 


(42) 


(43) 


(44) 


For the time history of drag pressure, a good fit to experimental data for TNT 
is a slightly modified form of that employed by Glasstone (ref. 34), 


q (t) = Q (1 - t/T) 2 e' bt/T 


(45) 


The procedure for determining the transverse loading blast parameters in 
Figure 34 which are independent of object size and shape is then as follows: 


(1) Obtain P and I from curves in Section II 

(2) Calculate P 

(3) Read b from Figure 35 

(4) Solve Equation (40) for T, knowing P, I. and b 

(5) Substitute P in Equations (42) and (43) to obtain P r and Q 

(fci Obtain P r and Q from Equation (44) 

(7) Substitute in Equation (45) for q(t), realizing that T - 1 3 - t a 
in Figure 34. 

The remaining quantities needed to define tin, time history of trans- 
verse pressure are dependent on the size and shape of the object. They are 
only well defined for objects of regular shape, such as right circular cylinders, 
flat rectangular strips, etc. Methods for estimating (Tj - t a ) and (Tj» - Tj) 
are given by Glasstone (ref. 34) for several such objects, and will not be re- 
peated here. One does need to know, however, the shock front velocity U. 


A2 


This is a unique function of the shock strength P, and is given by (see 
Chapter 6 of ref. 1) 


u 2 



( 46 ) 


Drag coefficients C D are available from Hoerner (ref. 35) for a variety of 
bodies over a wide r^nge of flow velocities. Estimates for the subsonic flow 
range which applies over the shock strengths of interest to us are given in 
Table VII. Melding these quantities, dependent on the size an [ ^ ° f * 
body, to the previous ones which are derivable from side-on blast P 

perties permits an estimate of the entire time history of transverse pressure. 

at least for bodies of regular geometry. 

e. Method of Predicting Appurtena nce Velocity - Once the time 
history of net transverse pressure loading is knot™, the prediction of ( appy- 

tenanc. velocity can be made quite easily. The basic assumptions are that the 

appurtenance behaves as a rigid body, that none of the energy in * h ' b >“* 
is absorbed in breaking the appurtenance loose from its moorings 
ing elastically or plastically, and that gravity effects can be ignored d " * 
this acceleration phase of the motion. The equation of motion of the object 

then 


A p (t) = M x 


(47) 


where A is the area of the object presented to the blast front, pltl is the 
net transverse pressure accoiding to Figure 54, M is t e to a mass 

object, and x is displacement of the object (dot. denote derivatives with re- 

spect to time). The object is assumed to be at rest initially, so t a 


(o) = 0. 


x (o) = 0 


(48) 


Equation (471 can be integrated directly. With use of the initial condition. 
(48). this operation yields, for appurtenance velocity. 


xlt) 


< T 3-‘a> . 

£ J P'" dl " M 

0 


(49) 


where 1,. is total drag and diffraction impulse. The integrations in Equation 
(49) can be performed explicitly if the pressure time history is described by 
suitable mathematical functions, or performed graphically or numerically if 
p(tl cannot be easily written in function form. In either case. Equation (49) 
yields the desired result — a predicted velocity for the object. 
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Cube, face-on 


Cube, edge-on 


Long Rectangular Member, 
face -on 
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edge -on 


Narrow Strip, face-on 
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D. Correlation of Velocity Prediction Methods With Data 

Prediction of PYRO Data by a Deterministic Method for C BM_ - A 
deterministic method of predicting initial fragment velocities was described 
in Section III. C. It was found that initial fragment velocities could be ob- 
tained by this method when some of the geometric characteristics of the ex- 
ploding vessel were known and when certain properties of the explosive pro- 
ducts formed upon detonation within the vessel were known. The data from 
the Project PYRO tests were not really sufficient to properly quantify the 
required independent variables for this method, because exact quantities of 
propellant which mixed and exploded were not known. It was felt, however, 
that certain reasonable values could be assumed for the unknown variables 
(along w'ith the measured values for other variables), and a feel for the appli- 
cability of this method could be obtained. Accordingly, results for initial 
fragment velocities were obtained from the computer program given in 
Appendix C when values for the independent variables to be described were 
used in it. These results were compared to PYRO data for mean initial fiag- 
ment velocities obtained from analysis of films of PYRO tests. Tne data from 
the PYRO tests were limited to those cases which best approximated the con- 
straints of the deterministic method. The correlation between the data and 
the theoretical values was good enough for "reasonable” values of the indepen- 
dent variables, that it may be concluded that a deterministic method of frag- 
ment velocity prediction may be useful if sufficient knowledge is obtained of 
the explosive product parameters in future testing. 

a. PYRO Data - The test numbers for PYRO tests considered 
in this correlation 7re given in Table VIII. Since our deterministic method is 
based on a spherically symmetric containment vessel, it was reasonable to 
take the cylindrical test geometry of the PYRO tests which most closely ap- 
proximated spherical symmetry, i. e. , L/D ratio closest to 1.0. Thus, only 
tests with L/D =1.8 (the smallest L/D ratio) were considered. Furthermore, 
only confined-by-the-missile (CBM) cases were considered since confinement- 
by-the-ground-surface (CBGS) cases do not represent internal explosions. To 
ass re that well- confined explosions were being considered, only D Q / D t 0.45 
cases were allowed. Finally, to simplify the assumptions made about the 
reactants and explosive products, it was decided to take the LH 2 /LO 2 cases 
only. 


Data for tests which meet these criteria are plotted in Figures 16 and 
37 for the 200-lb and 1,000-lb tests, respectively. In these plots, the mea- 
sured initial velocity is plotted versus a measured yield. The measured ini- 
tial velocity was obtained from measurements taken from films of the PVRO 
tests (see Section II. A). The measured yield Y M . is a quantity obtained in 
the following manner. Measured blast overpressures obtained along blast 
lines about the center of explosion for the various tests produced the PYRO 
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"yield" values, Y. This yield value is the percentage of the total propellant 
weight, W t , such that 

( 50 ) 

Y • 

where W XNT is the weight of TNT which would produce an equivalent sent- s 
of blast overpressure measurements in air to those of a given test. 1 e 
weight of propellant which would produce the same blast effect as TNI is 


Wr 


W TNT 


H TNT 

h R 


( 51 ) 


where W R is the weight of the reactants (or propellant) estimated to be in- 
volved in the explosion from the blast measurements, and H X pji' an R 
are the heats of explosion per unit mass for TNT and the reactants, respec- 
tively. For an LH 2 /LO 2 explosion in which fuel and oxidizer are mixe 

stoichiometrically. Equation (51) becomes 


w W TNT (52) 

Wr 3.7 

Measured yield. Y M . is the ratio of the weight of the reactants estimated to 
be involved in the explosion from the external blast measurements to the 
total known weight of propellant. Thus, from Equations (50) and (52), 

v W R Y (53) 

y m = Wt -3.7 

Values for Y M and W R appear in Table VIII along with measured initial 
velocity. 

Since measurements during the PYRO tests were all made externa l to 
the confinement vessel, it is not possible to estimate from them how much of 
the true blast yield was absorbed as kinetic energy of flying fragments. Thus, 
neither the PYRO yield value, Y, or the measured yield value, Y M . can be 
used to determine the actual quantity of propellant involved in the explosion, 
and W R is really only the quantity estimated by ignoring the effects of the 
confinement of the detonation. This is an important point, as will be seen in 
th ensuing discussion, because our deterministic method requires knowledge 
of the actual quantity of propellant involved and the characteristics of the re- 
sulting explosive products, and to obtain this information one is required to 
measure parameters internal to the confinement vessel, which was not done 
in the course of PYRO testing. 
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b. Data Predicted by Deterministic Method - in order to pre- 
dict initial fragment velocities with our deterministic method, it is necessary 
to have values for the following parameters related to the confinement vessel, 
explosion products, and propellant ullage: 

(1) The number of fragments, n, the vessel breaks up into, 

(2) The peak overpressure, P 0 , attained in the volume occupied 
by the explosive products immediately after detonation; 

(3) The ratio of specific heats, x, which describes the region of 
the explosive products immediately after detonation; 

(4) The speed of sound, a 0 , in the region of the explosive products 
immediately after detonation; 

(5) The radius, R Q , of the region of explosive products imme- 
diately after detonation; 

(6) The mass, M, of all the propellant not involved in the explo- 
sion process plus the mass of the confinement vessel. 

Physically, the explosion process is thought of as one in which a 
volume of the fuel and oxidizer mix stoichiometrically and then form a gas in 
the course of explosion whose volume and physical characteristics are des- 
cribed by parameters 2 through 5 above at the instant immediately after the 
explosion has occurred. All the rest of the propellant is thourht of as inert 
and merely augmenting the mass of the confinement vessel. 

As was seen earlier in this section. Figure 29, our deterministic 
method predicts that initial fragment velocities are insensitive co the number 
of fragments, n, when n is on the order of 10 or more. It is relatively 
insensitive down to n = 2. We picked n = 102 as a reasonable number of 
fragments for the PYRO tests investigated here. Also in this same section, 
Figure 30 indicates that above x * 1.2, the initial fragment velocities are 
relatively independent of the value of x. The parameter x ranges between 
1. 1 and 1.4 for most gases and tends to be lower for the^same gas at high 
temperatures (x = 1.3 for oxygen and hydrogen at 2,000 C, for instance, 
while x = 1.4 for these gases at 15°C). We picked x = 1.2 for our explo- 
sive products "gas", although a choice of a larger x would have changed the 
results very little. 

Using the value for reactant weight, W R , described under the pre- 
vious heading, we were able to calculate the energy of explosion for each test 
as estimated from the external measurements. This was simply 

E = W R • H r < 54 ! 
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Based on the calculated density of a mixture of L.H2/LO2 mixed stoichio- 
metrically, we obtained a volume of reactants related to Wr for each test 
from 


V 



where 


V HPH + v OX PQX 
V H + v OX 


(-5) 


and 

8 v hPh = v OX p OX 

where p and V are density and volume, respectively, and the subscripts H 
and OX refer to LH 2 and LO 2 * respectively. It should be stressed at 
this point that E and Vr refer to some weight of reactants that was 

thought to be involved in the explosion process from external measurements. 
Obviously, this E is less than the actual energy of explosion, since it repre- 
sents only the shock energy imparted to the air about the confinement vessel 
and does not account for kinetic energy imparted to the fragments, energy 
lost as heat and light to the air, energy lost through fragment heating, etc. 
The calculated volume of the reactants associated with energy of explosion E 
is also lower than the actual volume for the above reasons. Nonetheless, 
assuming these values for E and Vr are reasonable approximations for the 
real situation, for the stoichiometric mixture of L.H 2 /LO 2 one may obtain 
the peak pressure, P Q , from: 


P oo V R 

p o 

/ P o\ 

1/x- 

1 

X- 1 

__ Poo 

0 

0 

h 

1 


where P OQ is the ambient pressure (taken as 14. 7 psi). This was done for 
our values of E, Vr, and K using the Newton- Raphs on iteration technique 
for root finding. The program used (in the FORTRAN IV language) appears 
in Appendix D. The results predicted peak pressures of nearly 10 4 psi within 
the explosive products. Since E and Vr are both proportional to Wr, 
this result is independent of the weight of the reactants, and the same value 
(P 0 = 10 4 psi) was used for predicting fragment velocities from all PYRO 
tests considered. 

Picking a value for velocity of sound, a Q was complicated by the 
fact that this is a quantity that cannot be inferred from external measurements 
at all. It is an intrinsic property of the explosive products, a function of 
thermodynamic parameters such as temperature, density, and pressure, as 
wel 1 as the microscopic state of the "gas" at the instant of explosion. Since 
we had assumed gas-like properties for the explosion products, however, and 
since sound velocity would certainly be In the range 10 8 to 10 8 in. /sec for 
gases under most conceivable conditions, we calculated fragment velocity. 
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Uf, as a function of a Q and Rq for fixed n, M, Pq. K an 

estimate of the sensitivity of U £ to a 0 and Ro. relatively. This plot, 
appearing in Figure 14. shows that over the order, of magnitude chosen for 
a rt , Ur is relatively insensitive to a 0 and relatively sensitive to^. 
Accordingly, we chose a 0 = 10 4 in./sec as a constant value for al l tests. 
This value is a reasonable one for gases whose sound velocities vary in the 

range 1.0 to 5. 0 X 10 4 in./sec for the most part under normal conditions. 

Furthermore, our results would not be as sensitive to error in t 8 .? arai * 1 
eter as they would be to error in R^ If one assumes that the density of t 
explosive products gas is equal to the density of the LH 2 /L0 2 mixture deto- 
nated as reactants (i. e. , see Equation (55)) which formed it (certainly “ **' 
treme and unlikely case), nonetheless a c , for a given R<, calculated from 


a. 




(57) 


gives a fragment velocity very close to the same one obtained using a Q - 10 
in./sec (as seen from Figure 38). 


The mass, M. required by our deterministic method is obtained 


from 


M = M t + M^ a - M R 

where subscripts t. ta. and R refer to the total propellant weight, the tank, 
and the reactants, respectively. Again, since M R must be the mass of the 
actual quantity of reactants involved in the explosion process, it could only be 
estimated from the external PYRO measurements, that other internal mea- 
surements would be required to ascertain this quantity for each test. From 
the yield values it is apparent, however, that M R was probably small rela- 
tive to Mt and M*.. We estimated M R based on twice the calculated 
value W R . This produced corrections in M on the order of less than 10%. 
Values obtained for M are given in Table VIII. 

As was shown in Figure 38, the fragment velocities obtained from our 
deterministic method are most sensitive to values of R© chosen to describe 
the radius of explosive products at the instant immediately after detonation. 
Since Wn is obviously lower than the actual weight of reactants involved in 

the explosion process, at least a lower limit could be obtained for R*, baaed 

on 

( 4 /3)nRo 3 = V R = W R /p (59) 

from Equation (55). We assigned an upper limit on Rq based on the maximum 
radius a sphere would have which had the same volume as the confinement 
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vessel tank whose dimensions were obtained for the 200 and 1,000 -lb 
tests. Thus, Rq could be larger than the radius of the tank since the 
was a cylinder, but was of the same order of magnitude 


Within these brackets, an R*, was chosen to give a reasonable value 
for initial fragment velocities for one PYRO test for each of the tank sizes. 

P was then varied as the cube root of the measured yield for the other 
te° ~ :nd a corresponding U f was calculated by the deterministic method. 
Values for Rq are given in Table VHL Figures 36 and 37 show the com- 
parison of the measured and calculated fragment velocities for the 200 -lb and 
1,000-lb tests. 


c. Conclusions - Initial fragment velocities obtained by our 
deterministic method predict measured fragment velocities from Project PYRO 
tests relatively well when reasonable values are chosen for the unmeasured 
parameters. It should be noted that the values for Rq in Table VIII are all 
within a reasonable range relative to tank geometry. Our results indicate 
that a thorough knowledge of the internal parameters on just one test for a 
given tank geometry could probably allow reasonable predictions for initial 
fragment velocity, using this method, for all tests of that geometry based on 
external measurements alone (i. e. , internal measurements would only have 
to be made once). Finally, the great difficulty in any deterministic approach 
foi predicting fragment velocities for liquid rocket explosions (CBM) in 
general, is in obtaining data on how much of the reactants are actually going 
to be involved in the explosion. It seems that fragment velocity will be rela- 
tively independent of the number of fragments and even of the thermodynamic 
properties of the explosion products gas (except the peak pressure generated) 
but it is very sensitive to the quantity of reactants involved in the explosion 
process. 

2. Correlation of Velocity for Fragments from CBG S Tests with 
Predictions - The method given earlier in this section can, in 
theory, be used to predict velocities to which parts of the missile and appur- 
tenances are accelerated by explosions during CBGS tests of PYRO. In prac- 
tice, the method can only be applied to a limited number of the CBGS tests for 
which we measured fragment velocities. To use the prediction method: 


(1) We must know blast yield Y and total weight of propellant W T . 

(2) We must know the distance R of the fragment or appurtenance 
from the center of the explosion. 

(3) We must know the geometry and mass of the fragment, 

(4) Blast wave properties at the initial location of the fragment must 
lie within the range of the PYRO date. 
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We obviously satisfy the first requirement, but we can only estimate the aecon 
two. Also, the parts of the tankage which we believe constitute most of the 
visible fragments are located so close to the source that the blast is much 
stronger than for any of the PYRO blast date. Accordingly, only a rough 
"spot check" can be made to determine if the prediction method seems rea- 

sonable. 

The test chosen for prediction was Test 293, a LO 2 /LH 2 CBGS test 
with 1,000 lb of propellant impacting at 44 ft/sec. The measured blast yield 
was 3.9%. Twenty- six fragments were observed, with a mean velocity of 
464 ft/sec ± 226 ft/sec (Table V). From descriptions of the test apparatus 
and method in reference 16, the dimensions, material, and skin thickness of 
the tankage are known, as well as the rest position after drop impact. We 
believe that the majority of the observed fragments were relatively large 
pieces of the tank skins, made of 0.060-inch thick aluminum alloy. The lower 
part of the tank was about 3 ft above the center of the explosion, but the 
middle of the tank was about 10 ft from the center of the explosion. So, 
approximate values for input parameters are as follows: 


Blast 


R = 

10 ft 

w T = 

1,000 lb 

Po = 

14. 7 psi 

a o = 

1,088 ft/ sec 

Fragment 


Shape 

- circular disc of 12-inch diameter 

Thickness - h = 0.060 in. 

Density 

- p = 2.59 X 10’ 4 lb sec 2 /in. 4 


From Equation (3), Section II, the equivalent weight of propellant for blast is 
W = W T 1.000 . 39 lb 


Scaled distance is 

R/W 1 ^ 3 = lQ/39 1 ^ 3 = 2.94 ft/lb 1 / 3 


♦in the PYRO CBGS tests, the tank assembly was arrested by stops at the 
bottom of its fall after striking cutter blades which ruptured the tank 
bottoms. 






From Figures 22 and 23 in Section II, the blast pressure and scaled impulse 


P r 80 psi; l/w l/3 = 30 psi-msec/lb 
I - 30 X 3.4 = 102 psi-msec 
Dimensionless overpressure is 

P = P/ Po = 80/14.7 = 5.45 

From Figure 35, the time constant b can now be determined. It is 
b = 10 

From Equations (40), (41), and (44), we can calculate dimensionless reflected 
pressure J? r , drag pressure Q, and shock velocity, U. Multiplying by 
p o or a Q , the corresponding dimensional quantities are: 

P r = 49C p-i 

Q = 72 psi 

U = 3. 15 X 10 4 in./cec 

Solving Equation (40) for duration T, we have 


T = lb/ P [l - ^~~^j = 102 x 1(^80 £l - ^ 1 ~ r 14, 


2 msec 


To complete the description of the transverse pressure loading, we must know 
the drag coefficient C D and two characteristic times for the diffraction phase 
of the loading. The drag coefficient for a flat disc is, from Table VII. 

Cj) = 1. 17. The rise-time (Tj - T a ) is zero,* and 


<T 2 -T,) (-£) Jj 

where D is diamet r o i the disc. 

< T 2- T 1» * (t)(j.15 X 10*) “ °- 2 


msec 


*See Glasstone, ref. 34. 
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I d = 139 psi-msec 

Finally, predicted velocity is 

A Id 139 x 10~ 3 

U i = M 14 = OH = 2.59 x 10 - 4 x 6 x 10- 2 X 12 


Ui 


- 745 ft/ sec 


1 


Comparing this to the measured mean value. 


Uj = 464 ± 226 ft/ sec 


D 



..i- 1; •< sSilSSST: 

high. The blast parameters, and con.equen pred^ 

E. Frequency Distribution of Initial Velocity 

Since the number of fragment, which could be traced to each •••£«** 
f i f r nm Project P YRO T/ere relatively small, a logical group g 

r£5£3T2 order* to determine the initial velocity frequency di.tribu- 

tions. 

Since the data on fragment distance versus yield (%) in Section IV 
showed good correlation, a desirable grouping appeared to be over ^ medium 
yield (%? range within a configuration and propellant type. The groupings 
shown in Table IX and are: 


( 1 ) 


CBM, LO z /LH 2 , Yield from 10 to 29%. group sample size 10H. 


£<*> 
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(2) CBM, L0 2 /RP-1, Yield from 9.8 to 20%, group sample size 131. 

(3) CBGS (Vertical), L0 2 /LH 2 , Yield from 12 to 22%, group 
sample size 15. 

(4) CBGS (Vertical), L.0 2 /RP-1, Yield from 10 to 30%, group 
cample size 72. 

The data for * were ordered, and 10th through 90th percentiles (in 10% 
steps) were determined. 

The percentiles were plotted on normal and log normal probability 
paper. The log normal plots were the best fits. These plots are shown in 
Figures 39 through 42. The parameters for the log normal distribution can 
be estimated, according to Hahn (ref. 36) as follows: 

The mean, u, can be estimated as^the logarithm (to the base e) of the 
50th percentile. The standard deviation, o , can be estimated as two-fifths 
of the difference between the logarithms of the 90th and 10th percentiles. 

The estimated means and standard deviations are shown in Table X. Good- 
ness of fit statistics were calculated using the "W" test as described by Hahn 
and Shapiro (ref. 36). The method and calculations are described in 
Appendix F and discussed in more detail in Section IV. The calculated values 
of the "W" statistics for the four log normal distributions are: 

(1) CBM L0 2 /LH 2 - .925 

(2) CBM L0 2 /RP-1 - .947 

(3) CBGS L0 2 /LH 2 - .973 

(4) CBGS L0 2 /RP-1 - .989 

Referring to Figure 48, Section IV, the approximate probabilities for 
obtaining values as low as the calculated values of "W", given that the stated 
distributions are the correct ones, are: 


(1) 

. 415 

(2) 

.610 

(3) 

.915 

(4) 

.990 


*U A i* defined in Equruon (5), and is taken to be an approximation of the 
initial velocity. 
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TABLE IX. GROUPING OF TESTS BY PROPELLANT 
AND CONFIGURATION 


Propellant Configuration 


PYRO 
Test No. 


Yield. Y 

m 


Weight, W 
(lb) 


lo 2 /lh 2 


CBM 


22.0 

29. 0 

13.0 

25.0 

20.0 

17.0 

27.0 

10.0 


200 

200 

200 

200 

200 

200 

1,000 

1,000 


lo 2 /rp-i 


CBM 


48 

87A 

192 

193 

209 

270A 

278 

282 


9.8 

16.0 

14.0 

20. 0 

10.0 
13.0 
13.0 
13.0 


200 

200 

1,000 

1,000 

1,000 

1,000 

25,000 

25,000 


lo 2 /lh 2 


CBGS 

(Vertical) 


106 A 
115 
152 
184 
197 
211 
230 
266 
288C 


22.0 

15.0 

14.0 

17.0 

19.0 

12.0 
21.0 

14.0 

13.0 


200 
200 
200 
200 
200 
1,000 
200 
1,000 
25, 000 


lo 2 /rp-i 


CBGS 

(Vertical 


29.0 

10.0 

13.0 

14.0 


200 
200 
1, 000 
1,000 
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TABLE X. PERCENTILES, MEANS AND STANDARD DEVIATIONS 
FOR GROUPED VELOCITY DATA (fps) 


Percentiles 

CBM 

lo 2 /lh 2 

CBM 

lo 2 /rp-i 

CBGS 

lo 2 /lh 2 

CBGS 

lo 2 /rp-i 

10 

105.8 

290.0 

170.0 

160. 1 

20 

228.6 

414. 3 

219.2 

262. 3 

30 

393.3 

536. 8 

289. 3 

316.8 

40 

541.8 

714.2 

345. 1 

364.9 

50 

642.0 

822.8 

459. 4 

388.3 

60 

821.4 

973.3 

594. 3 

431.2 

70 

915.0 

1096. 3 

728. 8 

498.5 

80 

1102.7 

1216.8 

897.7 

573.4 

90 

1249. 3 

1405.6 

1169.9 

790.5 

Estimated 




5.962 

Mean* 

6. 464 

6. 713 

6. 129 

Estimated 

Standard 

Deviation* 

.9875 

.6313 

.7715 

.6387 


* Log Normal Diatribution, to base e . To convert to fps, take 
anti-logarithm 







Since a probability a- low a. . 10 i. conaldered inaolficient evidence 
to reject the choeen distribution., the fit. are assumed adequate, 
vation of data for Figure 48 in Section IV is given in Appendix G. 

The normal and log normal density function equations are given in 
Equations (62) and (63) in Section IV. 



IV. DETERMINATION OF FRAGMENT SIZE AND RANGE 
A. Retrieval of Data from Accidents and Tests 

All of the accident and test data which were retrieved were reviewed 
for pertinent information on fragment size, distance, and distribution. 

The criteria for each fragment were that the range, weight and maxi- 
mum projected area be specified. The nature of the test or accident, blast 
yield, type of propellants, etc. , were also of interest. 

From approximately 168 reports and memos reviewed, only eight 
events, listed in Table XI yielded sufficient information to meet the abov';- 
listed criteria. There were many other reports which had partial information, 
such as distance listed for fragments over 5 pounds, with fragments under 5 
pounds listed in number per square yard. 

The eight events listed in Table XI can be classified in three groups as 
follows: 

(1) Events 1 and 2 were Saturn IV confined by missile (CBM), 
LO^/LiH^ explosions. 

(2) Events 3, 4, and 5 were spill tests using three tanks, on 120° 
radials with LO z /LH 2 /RP- 1, and mixing on the ground (CBGS). 

(3) Events 6, 7, and 8 were mixing tests using two tanks with 

LO /LH and pouring the contents of one tank into the other. 

2 2 

Event 1 was Project PYRO test number 62, a Saturn IV vehicle with 
91, 000 pounds of L0 2 /LH 2 propellant, with self ignition, and a 5% yield. 
Complete details of the test can be found in Reference 16. 

Event 2 was also a Saturn IV vehicle, with 101, 198 pounds of L0 2 /LH 2 
aboard. The yield was estimated to be 1. 1%. Complete details of the mal- 
function can be found in Reference 38. 

Events 3, 4, and 5 are described fully in Reference 24. The tents were 
conducted using 1/25 scale quantities of the propellants for the Saturn C-2 
vehicle (i. e. , 7880 lb of RP-lj 32, 928 lb of LO,; 3032 lb of LH 2 ). The indiv- 
idual parts of the propellant were placed in spill tanks placed on 120 radials 
from some central location. The spill tanks were tipped toward this central 
location so that all three parts combined. They were then detonated from the 
central location by a small charge of C-4. 
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TABLE XI - CHART OF EVENTS 


Event 

No. 

Ref. 

No. 

Test or 
Accident Type 

Propellant 

Type 

Total Propellant 
Weight (lb) 

Yield Y, 
Percent 

1 

16 

FYRO Test #62 
(SATURN IV) 

L °2 /LH i 

91,000 

5.0 

2 

38 

S-IV 

L °2 /LH 2 

101, 198 

1. 1 

3 

24 

J2 Spill Test 

L0 2 /LH 2 /RP-1 1,754 

23.0 

4 

24 

J3 Spill Test 

L0 2 /LH 2 /RP-1 1,754 

24.4 

5 

24 

J 1 Spill Test 

LO./LH /RP-1 1,754 

b “ 

62.6 

6 

39 

Mixing Test 

lo 2 /lh 2 

240 

86. 0 

7* 

— 

Mixing Test 

“V LH 2 

240 

70. 0 

8* 

— 

Mixing Test 

L ° 2 /LH 2 

240 

73. 0 


* Data from these tests were furnished by the NASA test director, Mr. J. H 
Deese. There is no formal reference. 










The method for event 6 is described in reference 39. Events 7 and 8 
were similar. These were autoignition mixing tests using 240 pounds of laO^ 
LH , pouring one tank into the other. Photographs of each fragment, along 
witfl fragment maps and fragment weights were supplied by the NASA Test 
Director, Mr. J. H. Deese. 


B. Data Reduction 

The data from each of the events were reduced by careful analysis of 
each fragment to assure that size (maximum projected area), weight, and dis- 
tance were specified. In some cases, it was possible to fill in missing items 
by estimating weight and/or size from information supplied by descriptions 
or photographs. Because of the paucity of fragment data, considerable effort 
was expended to extract as much fragment data as possible, without undue loss 
in accuracy of the parameters. Distances for each fragment were determined 
from fragment maps. 


For each event, data including event code, fragment number, distance 
(R) in feet, weight (W) in pounds, width in inches, length in inches, maximum 
projected area (A) in square inches, area divided by weight (A/W), and drag 
coefficient (CD) were entered on key punch sheets. The drag coefficient was 
estimated from photographs and descriptions, and subsequent comparison with 
Table VII. Cards were keypunched and used in the analyses which are de- 
scribed in the following sections. A listing of all fragment data by event is 
available if desired. 


C. Statistical Studies 

1. Computer Routines - Using the Biomedical Computer Programs 
(Ref. 40), with a CDC 6500 computer- the data from each event was subjected 
to the following routines: 

(1) Means and variances were calculated for R, W, A, A/W, 

CD, Log R, Log W, Log A and Log A/W. 

(2) Histograms were constructed for the parameters listed in (1) 
above. 

(3) Correlation plots were made for R versus W, R versus A, 

R versus A/W, R versus Log W, and R versus (A/W) CD. 

2. Preliminary Analysis - The output from the correlation routine 
was studied to determine if there were discernable patterns of correlation 
between parameters within an event. While there was some general pattern 
in some cases, as a whole, the scatter was so great as to discourage further 
inquiry along these lines. This result could be explained by not considering 


(because of a lack of knowledge) the flight angle and initial location of the 
fragments. 

The histograms were studied to relate the parameters R, W, A, and 
A/W to probability frequency distributions. Since the sample size varied 
from 31 to 1056. the histograms varied in information content in about the 
same ratio. However, the form of some of the histograms suggested that a 
normal probJbility density function would supply an adequate fit, and others 
offered the possibility of a by a log normal distribution. 

The normal distribution can be written as: 

f (x) = (1 h/Tri o) exp [- (x-g) 2 /2o 2 ] . - oo< x < <*> * 


- oo < U » c < oo , 

where 0 is the standard deviation and U is the mean. 

The log normal distribution can be written (ref. 36): 
f (y) s (Ihjzrra v)exp ^(iny - w ) 2 /2a 2 )] . 0 <y<ao, 


- oo > M < 00, 0>O, ' ' 

where a is the standard deviation, H is the geometric mean, and x is iny. 

The Weibull distribution was also considered, but later results showed 
that better fits were obtained by the normal and log normal distributions. 

Table XII presents the mean, standard deviation, standard error of the 
mean, sample size (number of fragments), maximum, minimum and range 
(maximum value minus the minimum value) for R, W, A, A/W, Log 1(J W, and 
Log 10 A. 

3. Probability Density Functions - The data for each parameter of 
interest within an event was sorted in ascending order and the value for the 
parameter for the 10th to the 90th percentiles in 10% steps was identified. 
Table XIII is a listing of these values. The use of order statistics tends to 
equalize the effects of varying sample size from event to event. 

A plot on normal and log normal probability paper was then made for 
each parameter for each event. Figure 43 is a plot of distance for event 3 
on normal probability paper. Figure 44 is a plot of the same data on og 
normal probability paper. Since the points on the normal probability p-per 
lie closer to a straight line than those on the log normal paper, the normal 
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TABLE XIII - PERCENTILES FOR PLOTTING PARAMETERS 
OF EVENTS 1 THROUGH 8 


Weight (W) 
(lb) 


Area (A) Area/Weight (A/W)| 


(in 2 ) 


(in /lb) 



120 . 0 
180 . 0 
216 . 0 
360 . 0 
576 . 0 
864 . 0 
1152.0 
2304 . 0 
4320 . 0 


26.0 
44 . 0 
55 . 5 

67.2 

76.0 

88.0 
100.0 
120. 0 
144 . 0 

1.2 
6 . 1 

19.3 
33.6 
43.2 
55.9 
72 . 0 
90. 0 

113 . 1 


3 . 3 
3.8 
4.4 
7 . 5 
14 . 0 
26.2 
27.9 
28 . 1 
28.3 






Percent 

«tance (R) 


<ft) 




925 

1035 

1120 

1220 

1730 

45.0 

82.0 
106 . 0 
127 . 0 
148 . 0 
169 . 0 
192. 0 
233 . 0 
281 . 0 



* t+jhJI 


15.0 

18.0 
2 l . 0 

48.0 

55.0 

65.0 

97.0 
110 . 0 
125.0 

. 005 
. 010 
. 018 
. 028 
. 044 
. 069 
. 108 
. 210 
.627 

. 025 
. 051 
. 090 
. 100 
. 100 
. 135 
.209 
.400 
1.833 

.022 
.039 
. 064 
. 100 
. 113 
, 193 
. 380 
. 721 
2 . 009 


72.0 
144 . 0 
192 . 0 
240 . 0 
288 0 
360 . 0 
432 . 0 
576 . 0 

1296. 0 

0.4 
0 . 5 

1.3 
2.0 
2.7 
4 . 1 

6.3 

12.0 

24.0 

1.0 

2.0 

3.0 

4.0 

6.0 
8.0 

14 . 0 

24.0 

60.0 

1.5 
2 . 5 
3.8 

5.5 

8.0 

14.0 

22.8 

40.0 

90.0 


Area/Weight (A/W) 
(in^/lb) _ 


2 . 1 
3 . 5 
3.6 

3 . 7 
3.8 

4 . 0 
7 . 2 

13 . 4 
27 . 0 


14 . 286 
27. 996 
54. 762 
70 .. 000 
80. 176 
88. 888 
100. 000 
119. 469 
160 . 894 

10. 000 
20 . 000 

30. 000 
45 . 557 
57. 600 
69. 767 

80. 000 
93 . 023 

127. 660 

10. 000 
28 . 000 
44 . 280 
61 . 303 
78. 571 
96 . 154 
103. 131 
125. 000 
166.667 
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FIGURE 43 . EVENT 3 PROBABILITY DISTRIBUTION (NORMAL) , DISTANCE 
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FIGURE 44. EVENT 3 PROBABILITY DISTRIBUTION (LOG NORMAL) , DISTANCE 









distribution appears to be abetter fit. Figures 45, 46 and 47 are examples of 
the other plots. From the plots over all 8 events, it appeared that the normal 
distributions adequately fitted the distance (R), and A/W and that log normal 
distributions best fitted the weight (W) and area (A). A complete summary of 
the fragment data is given in Table XIV, giving the estimated standard devia- 
tion (S), and mean (M) for the respective distributions for each parameter in 
each event. 

A "W" statistic for goodness of fit for each parameter was calculated 
using the methods outlined by Hahn and Shapiro (Ref. 36). 

The approximate probability of obtaining the calculated test statistic, 
given that the chosen distribution is correct, was then determined, and is 
shown in Table XV. Figure 48 is a plot of Table XV, and can be used to de- 
termine the approximate probability of obtaining a value as low as the calculated 
value of "W", given that the selected distribution is the correct one. The 
values in Table XVI were calculated using the method and formula outlined by 
Hahn and Shapiro (Ref. 36). 

Table XVI is a summary of the "W" test on normality for R, W, A, and 
A/W for the eight events. The method, and sample calculations for the W 
statistic and the probability of obtaining the calcul ted value of W are given 
in Appendices E and F. 

Referring to Table XVI, we see that there are 32 distributions, one 
each for R, W, A, and A/W for each of the eight events. The probability of 
obtaining the calculated value of "W M is greater than 50 percent for all except 
the A/W distributions for events 3, 4, and 5, indicating adequate fits for all 
except these three distributions, as it is customary to consider values exceed- 
ing 2 to 10% as adequate grounds for not rejecting the hypothesis that the data 
belong to the chosen distribution. It is interesting to note that each of the 
parameters i 3 distributed in the same family (i.e. , normal or log normal) 
across all eight events. That is, distance (R) has a normal distribution func- 
tion in each of the eight events, indicating a repeatable pattern. 

The estimate for the means and standard deviations for each of the dis- 
tributions is given in Table XIV. 

D. Methods of Prediction of Range Versus Fragment and 
Blast Yield Parameters 

1, Determination of Mean Range of Fragments Versus Blast Yield by 
Regression Analysis - The mean distance R versus the yield Y in percent and 
equivalent pounds of” TNT, W T for the 8 events was plotted on log-1 ’ 
paper and is shown in Figures 49 and 50, respectively. As can be seen, 
when events 6, 7 and 8 which are the mixing tests are excluded, there is more 
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scatter on the TNT chart than on the percent yield chart. A regression 
equation was derived to describe a linear fit to the points. Events 6 . 7 and 
8 were excluded on the basis that the propellants were mixed differently than 
in events 1 through 5, The equation is: 

R = 314.74 Y 0, 2775 

This equation should be limited to the range of weights and configurations of 
events 1 through 5. 

The line is drawn on Figure 49, and the predicted versus observed 
values of the mean distance for the 5 events can be read from the chart. On 
Figure 49, the maximum observed distance is also shown for each event. 

In Figure 51, the upper dashed line shows the estimated distance which 
should contain at least 95% of the fragments. Table XVII shows the upper 95% 
confidence limit (CL) on the estimate of the mean <M), the upper 90% confidence 
limit on the estimate of the standard deviation, and the various quantities 
necessary to calculate these confidence limits# 

The confidence limit on the mean was calculated using the following 
formula: 


CL * M +■ 


r- t ; 95, 
rn n 


n is the number of fragments and t n ; 95 is the value of the t distribution with 
«i degrees of freedom at the 95th percentile. 

The confidence interval for the standard deviation was calculated using 
the following formula: 


CL * 


[ EX 2 ! - (IX i) 2 / n 1 

L X (n- 1 ); 90 J 


Where Xi is the distance of the ith fragment, n is the number of fragments, 
and X 2 . j. 9Q is the value of a chi square distribution with n- 1 degrees of 
freedornat tfie 90th percentile. 

Then, using the new upper confidence level values of M and S, the 
range 1 l 95 in which 95 % of the fragments should fall was calculated as follows: 

R 95 8 M + S *n; 95 

The interval from the mean (M) to H 95 is indicated for each event on Figure 51 
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by a bar. 

v. j parallel to the regression line , and just touching 

A line was then drawn p , this line could be expected to 

2. Determinati o n of Mean Ran B* of F * Ll^ l ! Jlo M nn TT- 7 hp 
Equations - The range of a fragment produce y ; ^^ afion . g avai i ab le 
determined from ballistics equations if sufficient info Specifically, 

about the initial fragment velocity and fragment characteristic. pe 
it is necessary to know the following parameters. 


( 1 ) 

( 2 ) 

(3) 

(4) 

( 5 ) 

( 6 ) 


The fragment mass, M 

The initial fragment velocity, U f 

The initial elevation of fragment trajectory, 0 

The fragment cross-sectional area, A 

The drag coefficient for the fragment, C D 

The density of air, P . 


Having given these parameters, ttal tlme of il.ght of the 

fragment. This method »«•«»££ A missile exploding 

derived by J. J. Oslake, . density remains constant along 

on the launch pad with the assump gravitational forces 

the trajectory of the fragments produced. Only drag a g 

are assumed to act c the fragments where 


pa. (1/2) p U AC , Fg = -Mg 


(65) 


u . . raf thP drap force and the gravitational forces, 

are the vertical u" = vertical velocity of the fragment. > 'ine 

I^SioHf motion fo* vertical Icceleration may be integrated to obtain t 

as a function of U , which gives 


t = 4— tan '' (U 2 o‘'/ r7 ® ) 

r -ycg 20 

for the case in which U z = 0 < l * e * * t i 
In Eq. (66) U zq is the initial vertical 


( 66 ) 


is the time of rise for the fragment), 
fragment velocity and 
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(67) 


p A C 

id 


The equation of motion for the fall of the fragment may be similarly integrated 


to obtain 



i n 


2CZ M 

(2 e 


- 1 ) + 



2 CZ M 
e 




( 68 ) 


where 



The total time of flight of the fragment is 


(69) 


where t and t. are functions of initial conditions only. Finally, the equations 
of motion for racial acceleration of the fragment may be integrated to obtain 
the range R in terms of T t 


R* 7 in <1 + c u ro t > 


(7 1) 


The range may be determined from our 
function of M, A, C , P a £ r (^*1* (^7) )> 
given initial velocity, U^: 


initial parameters, since c is a 

U is the radial component of the 
Ro 


V 


Ro 


= U f cos 0 


V = U, sin 
Zo f 


(72) 


and thus is a function of U f and 0; and T is a function of initial conditions U f , 
0, c and g. 


A computer program in FORTRAN IV which computes fragment range 
as a function of the independent parameters mentioned above is given in 
Appendix H. This program was checked against the results of J. J. Oslake, 
et al (Ref. 41); the check appears in Appendix H. The symbols used in the 
program as well as the input and output parameters are also discussed there. 

This method of range prediction requires considerable information 
about the fragments emanating from the missile explosion. In general, the 
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initial fragment velocity, U , may be determined by one of the methods describ- 
ed in Section III. B. Fragment mass, cross-sectional area, and drag coeff 
must be determined empirically. Drag coefficients for fragments at various 
velocities are described by L. D. Heppner and J. E. Steedman (ref. 42). 
Coefficients for regular shaped projectiles (spheres and cylinders) are e- 

scnbedby E. Richards ( ref. 43). A correlation between range predicted y 

this method and that determined statistically for a "mean (or average) rag- 
ment from the PYRO test data, event 062, is given in Section IV. E. 


nor relation 


1. Standard Statistical Tests 

a. Standard Deviation Versus Mean Distance - Figure 52 is a 
plot of the standard deviation of the distance versus the mean distance R for 
each of the eight events. As one might expect, the standard deviation in- 
creases as mean distance increases. 


b. Results of Scaling A/W Parameter - F rom Appendix B , 
formula B 13 is: 


G 7 ’ = A f /W f 


2/3 


(73) 


where A f is the fragment cross-sectional area (in 2 ), and W f (lb) is the frag- 
ment weight. 


Figure 53 is a plot of the geometric means of A/W versus geometric 
mean distance R (ft) for the eight events. Figure 54 is a plot of the geometric 
means of A/W 2 ' 3 versus geometric mean distance. 


The points, excluding events 3, 4 and 5 have been fitted to a line with a 
regression analysis which resulted in the prediction equation. 


R « 9. 864 (A/W 2/2 ) °* 78 


(74) 


2 

where R is the mean range (ft), A is the area of the fragment (in ), and W 
is the weight of the fragment (lb). 


The exclusion of events 3, 4 and 5 can be justified on the basis that 
A/W for these three events were noticeably low in probability of belonging to 
the log normal family of distributions, while A/W for the other events showed 
a high probability of belonging to the log normal family. 

Comparing the points in Figure 52 with those in Figure 54, it is ob- 
vious that, in this case, scaling of the A/W parameter improved the correlation 
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( 1 ) 

( 2 ) 

(3) 

(4) 


The mean weight, W s 6.797 lb 

The initial mean velocity, U f - 741 ft/&ec 

O 

The mean elevation angle, 0 = 12 - 77 


The mean cross-sectional area, A - 4.0 ft 


( 5 ) 


The mean range, R = 447. 4 ft 


w r 

, were obtained from 

referent Z 43° 

ZSLXZX.'Z SS2 >im " pro ' 

jectiles, ( ref. 42). 
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It can be seen from these results that reasonably accurate range values may be 
obtained if an appropriate is chosen. 

This range prediction method could be more useful if the fragments 
were broken down into classes of narrow mass ranges whose mean velocity, 
elevation angle, etc. , were obtained from empirical considerations. Then 
a predicted range mapping could be made, as a function of yield, with the use 
of these ballistic equations. To do a complete empirical study, however, 
would probably require arena tests using models of the actual missile tanks. 


134 


V. EFFECTS OF FRAGMENTS 
A. Introduction 


It is desirable to have the capability to predict the probable damage 
levels to humans, structures, vehicles, etc. at various ranges from an ex- 
ploding missile. Damage can be produced due to the blast wave emanating 
from the rupture of the missile or the interaction of the blast wave with ob- 
jects surrounding the missile. 

Damage caused by blast waves has been extensively studied and is a 
result of the peak overpressure or impulse of the blast wave depending on the 
nature of the "target". Ref. 34 gives examples of this kind of data for damage 
to humans as a function of peak overpressure. Since peak overpressure is a 
function of range and yield, damage levels can be associated with these para- 
meters. 

In order to predict damage as a result of fragmentation, it is necessary 
to know the value of the terminal ballistic parameters of the fragments (i.e. , 
velocity, range, mass, cross-sectional area, etc.). Insofar as the value of 
these parameters can be determined from the characteristics of the missi e 
explosion, damage levels can be predicted as a function of yield. 

In Sections II and IV of this report, we have tried to investigate 
methods of determining fragment initial velocity and range by statistical and 
deterministic methods. By and large, however, the empirical data are not 
sufficient to cover the fragment damage problem on the basis of existing mis- 
sile explosion data. A great deal of work has been done, however, on the 
problem of fragment damage from bomb explosions. We feel that much of this 
work is applicable to the exploding missile problem, even though the initial 
fragment velocities for the exploding missile would be primarily subsonic, 
whereas those for exploding bombs are usually supersonic. 

A considerable amount of investigation has been made of what the ter- 
minal ballistics characteristics of a fragment must be in order to cause 
damage to various "targets". Equations for penetration of wood, steel, etc. , 
of various thicknesses (as well as human simulators) have been developed 
from empirical data which may be used to predict damage levels as a function 
of fragment mass, velocity, and area. The probability of damage to these 
targets at any range would be a function of the probability of arrival of a frag- 
ment meeting the minimum specifications for the damage equations. The 
probable fragment damage for a given range as a result of a missile explosion 
is thus predictable when these two parts of the problem are solved: 

(1) The probability of arrival at any given range of a fragment 
of specified mass and cross-sectional area, velocity, etc. 
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(2) The probability of penetration or damage to a specific 
"target" struck by a fragment of the specified mass, 
cross-sectional area, velocity, etc. 

B. The Probability of Arrival of a Fragment of Specified Characteristics 

Versus Range 

1. Fragment Character istics - It is necessary to characterize the 
fragments in various ways such that trajectory analysis will allow a map ping 
of fragments from a unit solid angle about the explosion to a unit area on the 
ground surface at some range, R, from the explosion center. Specifically, 
it is necessary to know the fragment mass distribution; cross-sectional area 
or projected area distribution, and the drag coefficient for the fragments. 

Many of these characterizations are obtained from empirically derived 
equations for bomb studies. One accepted equation in bomb studies relating 
fragment mass and area is: 


, >, . ' 

* - & c 
. ?■ V 




M = k A ' 

where M is a fragment mass, A is its cross-sectional area and k is some 
constant ( ref. 44). For a mild steel fragment of "flattened shape k - 1. 45, 
for instince(whe re A is taken as the mean projected area of the fragment) 
(ref. 45). The distribution of fragment masses may be given by: 

N (M) = N o exp (M/U (76) 

where N is the number of fragments with masses greater than M, and 4 
equals 1, 1/2, or 1/6 the average fragment mass depending on whether Y 
equals 1, 1/2 or 1/3, respectively ( ref. 46). Equation (76) actually thus 
defines three different mass distribution laws which are applicable for various 
types of explosions (for a thick walled shell in a three-dimensional breakup, 

4 equals the average fragment mass and Y s 1, for instance). 

Appropriate values for the drag coefficient must be chosen in order to 
calculate the aerodynamic coefficient in the drag force equation, c, [see 
Section III.B-1, Equation (67)1 . Some values for C D were given in Section 
IV-C. For bomb tests, the values Cp 3 .48, .60 for sub- and supersonic frag 
ment velocities seem to be common (Refs. 46 - 48) giving: 


F , * . 0031 m 
d 


-1/3 


subsonic 


- 1/3 2 

F . * . 0039 m U supersonic 
d 


for the drag forces (mass in engineering units and velocity in ft/ sec. ) 


Empirical equations such as those obtained Jrom bomb studies should 
be developed for missile explosions where the fragments are arger an av 
lower initial velocities. Until such equations are developed, some good Pre- 
dictions can probably be made using the equations relating to subsonic velocities 
for bomb fragments of relatively large mass. 

2. Tr ajectory Analysis - A trajectory analysis may be carried out on 
fragments of known characteristics as was demonstrated in Section lV-. 

These techniques result in a range calculation but not a terminal veloc. y cal 
culation. Solution of the equations of motion for a fragment traveling at the 
speed U along its trajectory: 


X + 3 V X + g sin Ci = 0 

• • • 

"V + 3 v Y + g cos a = o 


(78) 


where 0 is the aerodynamic_coefficient with the samedefinitionasc m 

Equation (65). Sections II-B, X and Y are local coordinates of the frag 
at time t measured in a coordinate system moving with the fragment^ and a 
is the elevation angle of the fragment relative to the trajectory in the '<• 
coordinate system) will yield the terminal velocity and range. These equat 
are solved by Feinstein and Nazooka (ref. 48) using a technique in which 
gravitational effects are a perturbation on the drag equations. 

jj _1 2 (79) 

X + 3 X =0 1 ' 

o 

For trajectories in which the fragment, had an initial elevation angle less than 
that required for maximum range (low register trajectory) 

(80) 


X = (X Q + X p ) cos a o - Y sin 


a 


where 


X = An (1 +u )/ 0 
o 

X p = - (g/2) t 2 sin a (1 + u/3)/ (1 + u) 
u = 3 V t 

O 

"y * (g/2) t 2 cos a[u (1 + u/2) - log (1 + u)] /u 
t is the time of flight obtained from the equation 
Y**X +X sin a + Y co* a ■ 0 

n n O O 


(81) 


(82) 
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and a = a is the initial elevation angle. Solutions for other trajectories are 
obtained numerically for time interval steps t in which displacements occur. 


AX = X cos a - Y sin a 
AY = X sin a + Y cos a 


(83) 


Terminal velocity can be obtained by straightforward differentiation of Equa- 
tions (78) and 179) with respect to time, t, but Feinstein omitted this step. 


3. Mapping From Explosion Center to the Ground Plane - Feinstein 
uses Eqs.(80) through (83) to map the number density of fragments from^ombs, 
n., in a mass internal i from those passing through an area element R 0 cos a Q 
A l 0 A on a hemisphere immediately surrounding the bomb, to those land- 

ing in an area on the ground plane at range R, R A 0A R. For nj density 
of fragments at the elemental area at the hemisphere and m density of frag- 
ments at the elemental, area at the ground plane 


R cos a 
o o o 

n i " n i R I dR/d a 

o 


(84) 


The trajectory equations and initial conditions on the fragments yield the value 

for R + dR/d a . n.°, R , and a are related by fragment arena data, 

o i o o 

In a simpler view, n.° could be obtained from Eq. (77) and the explo- 
sion could be assumed to have spherical symmetry. As in Feinstein' s calcu- 
lation, Eq. (75) would be assumed to hold. Values for C D would be chosen 
appropriate to the type of fragments expected from missile explosions and the 
number of fragments in a mass interval classification as a function of range 
would be calculated using Equation (84) and trajectory equations. Finally, the 
probability of impact on a target of area A^ of a fragment in mass internal i 
may be obtained from Poisson statistical equations, and the probability of 
damage for a given "target" would be assessed on the basis of vulnerability 
criteria and the flux of fragments as a function of range meeting those criteria. 


Feinstein has actually solved the problem of hazards to humans as a 
result of several types of bomb explosions using his technique. His results 
(for the assumption that fragments of 58 ft-lb or greater of energy are hazard- 
ous to humans) are given in the form of probability of injury contours. We 
suggest that this technique or a simplified version of it coupled with more com 
plete data on missile explosions could be used to produce similar results for 
the missile explosion problem. 
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C. Vulnerability Criteria 


1. Structures - When a fragment emanating from an explosion arrives^ 
at the "target", the degree of damage it will produce is generally a function o 
its terminal velocity, mass, cross-sectional area, and the characteristics of 
the target. The relationships between these quantities and some measure of 
damage level for .tructurea are empirical equation, which are unrelated to the 
aource of the fragment. Thus equations of thia type which have been exten- 
sivelv developed tor prediction of damage from bomb iragment3 would be 
applicable to the exploding missile problem. 

An accepted form for vulnerability equations is: 


U = k 


a P 

e A 

y- 

M 


(85) 


where U is the fragment velocity. A is the presentation area of the frag- 
ment, M is the fragment mass, and e is the thickness of the target mater- 
ial. This equation of the "de Marre type" predicts the terminal v e loclt Y of a 
fragment required to penetrate the distance e into a structure where the con- 
stat a , P , Y , and k are a function of the material from which the struc- 
ture and fragment are made. These constants are obtained empirically for 
various materials, and example values for the constants can be seen in 

Table XIX. 


Some fragments from missile explosions are sufficiently large that 
penetration of a structure is not likely to be the primary damage mode. For 
surh fragments the kinetic energy of the fragment may be of most interest 
in relating to probable structural damage. In some cases, the mass of the 
fragment would be sufficient to crush a structure were it to land on top of it. 
Clearly the damage potential of fragments from missile explosions is related 
to the types of structures which are in danger of damage and the materia s 
and techniques used in their construction. To classify and investigate all 
such structures is beyond the scope of this report. 


2. Personnel - The simplest criterion for fragment damage to per- ^ 
sonnel is the so-called "kinetic energy criterion" which is attribed to Cranz 
"Lehrbuch der Ballistik". This criterion states that a fragment having 58 ft-lb 
of kinetic energy is capable of producing a human casualty. This criterion 
was widely used during World War II to predict casualty levels and was used 
by Feinstein ( ref. 44) in 1972 to obtain probability of casualty contours from 
the equations described in Section V-B. Like most empirical values its 
justification is that it seems to give correct results more often than not. 
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Other empirical values for the characteristics of fragments which 
will produce an incapacitating wound have been obtained from 


Mj U 


a 

1 


a 

= m 2 u 2 


( 86 ) 


where M is the mass of one of two fragments, i, and U i is its velocity. 

This equation attests that having determined the mass and velocity of one 
fragment capable of incapacitating personnel all other fragments satisfying 
Eq. (86) also would be potentially incapacitating. The value for a is various 
ly given as 2. 5, 3. 0, etc. 


A "standard" assessment for fragments capable of producing casualties 
which has been widely used, is the wood penetration test. In this test, pine 
wood panels of thicknesses ranging from 1/2 to 1 inch are present in the range 
of fragments from an explosion. Where the panels are perforated, it is 
assumed a casualty would have occurred, otherwise not. 


At this time, it wculd seem that acceptance of an empirical "standard 
such as the kinetic energy criterion of 58 ft-lb would be th- best idea in inves- 
tigating potential fragment damage to personnel as a result of a missile exp o- 
sion. Acceptance of this value would allow the use of bomb fragment data 
which was based on that value. Ultimately, it would seem reasonable to try 
and obtain data on the probability of injury as a function of kinetic energy 
of the fragments, and extend existing bomb fragment data to the case of more 
massive fragments of the type which may be encountered in missile explosions 


VI. DISCUSSION OF RESULTS 

A major part of the work consisted of efforts to retrieve data relevant 
to the basic problem of fragmentation of liquid propellant vessels by literature 
search and visits to various government and private agencies. We believe 
that we were successful in assembling the majority of such data which are sti 
available in the literature and in the files and archives of NASA and Department 
of Defense agencies, and in entering the data in the ASRDI data bank. The 
primary sources of data were the Project PYRO experiments (refs. 15-17). 
recent experiments conducted at NASA-Kennedy Space Center (ref. 39), an 
several accident reports (refs. 24 and 38). These sources, however, yielded 
only partial information with one exception. They gave either data on initial 
fragment velocities together with measured blast yield, or data on final loca- 
tions, masses, and shapes of fragments. The only exception was PYRO Test 
#62, which gave both types of data. 

The following observations are appropriate in order to place the results 
obtained from fragment velocities measured from PYRO films into the proper 
perspective. 

(1) Velocities were measured for many different combinations of 
tank size (and total propellant weight), tank height to diameter 
ratio, propellant types and test type, that is, confined by 
missile or confined by ground surface. 

(2) Velocities were measured at different positions in space de- 
pending on where they were first sighted. The instantaneous 
velocity of the fragment immediately a'Ler sighting was deter- 
mined, and an average velocity of the fragment was determined 
by assuming that it traveled in a straight line from the center 
of the tank to a point at which the fragment was first sighted. 

No attempt was made to estimate the initial velocity of the frag- 
ment at the tank wall from the velocities calculated at a later 

time. 

(3) Velocities determined from data obtained from a single camera 
view are only approximate because the fragments often had 
components of velocities directed toward or away from the 
cameras which could not be detected. The necessary assump- 
tion that such fragments travelled normal to the camera lens 
axis introduced possible errors of +0% to -30%, i.e. , velocities 
were generally underestimated. Likewise the coordinates of 
the fragment in terms of the radial position, height, and azimuth 
angle could not be determined accurately. 


142 


(4) For fragments which could be observed from two different 
camera locations (spaced 120° apart), their true trajectory 
and thus their true velocity could be obtained. Unfortunately, 
fragments were difficult to identify in more than one camera 
view and so true velocities were obtained for only a few frag- 
ments. 

(5) It was rot possible to relate on a one-to-one basis fragment 
velocity and heading with fragment size, shape or final range. 

In spite of these limitations to the fragment velocity data, we feel that 
the data we obtained were very useful, and reasonable correlations with pre- 
dicted values of velocity have been obtained. Nevertheless, we feel that the 
velocity data warrant additional study directed perhaps toward predicting the 
initial fragment velocities and those we had measured. We also feel that con- 
trolled experiments should be performed which would permit correlation be- 
tween measured fragment velocities and the resulting missile masses. 

Techniques could be employed to allow initial velocities to be correlated with 
the final fragment size, weight, and position. An overhead camera located 
above the event with a wide field of view should facilitate velocity measure- 
ments and tracking during such a test. Tag markings would also aci l a e 
fr igm^nt identification. 

The fragment initial velocities, as estimated by the calculated para- 
meter U, (see Section III) for the four pooled groups (CBM-LC> 2 /LH 2 , CBM- 
LO /RP-1, CBGS-LO./LH,, CBGS-LO ? /RP- 1) were fitted by log normal 
distributions. The means and standard deviations of these distributions are 
given in Table VIII. 

The goodness of fit test, using the "W" statistic (ref. 36) showed pro- 
babilities from .415 to .99 of obtaining the calculated values of W, giving 
us no cause to reject assumed distributions. Thus, the log normal distribu- 
tions appear to be an adequate fit for the fragment initial velocity distributions. 

These probability density functions can be used to represent fragment 
initial velocity densities for the subject propellant combinations within the two 
configurations , CBM and CBT S. 

Complete details of all data sources and analyses are given in the main 
body of the report. 

In Section II, we reviewed duta pertinent to blast yield for propellant 
explosions, and also presented methods for estimation of blast yield and the 
free-field blast wave properties for various types of explosions and propellants. 
These methods relied heavily on the extensive series of test results from 
Project PYRO (refs. 15-17), but also Incorporated much of the physical 
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insight apparent in the work of Farber and co-workers (refs. 27-30). The 
methods allow estimation of blast wave peak overpressure and impulse, 
given propellant type, total propellant weight, type of accident, mixing time, 
and distance from the center of the explosion. Predictions are limited, how- 
ever, to the range of scaled distances of the PYRO tests. A major difference 
in our prediction methods and those of reference 17 is in choice and scaling of 
ignition time. We believe that the scaling for this time assumed in references 
15 and 17 is not verified by the test results, and discuss this point at some 
length in Appendix A, together with a dimensional analysis of ignition time 
scaling. 

Density functions for fragment distance R, weight W, area A, and 
area/weight A/W, were estimated for each of eight events. A complete 
listing of the form and estimated parameters is given in Table XIV. Goodness 
of fit tests using the ,: W" statistic supported the selected density functions, 
in all but 3 cases out of the 32. These three cases were A/W distributions for 
events 3, 4, and 5. 


The form of the density functions was constant for each parameter. 
That is, for all eight events, R and A/W followed normal density functions, 
and W and A followed log normal distributions. 


Distance R versus percent yield Y for events 1 through 5 shows a very 
good correlation in Figure 48. A prediction equation was derived, and is: 


A 

R 


= 314. 74 Y 


0. 2775 


(63) 


Confidence intervals were calculated for the means and variances for R for 
each of the eight events, and are shown in Table XVIII. In addition the 95th 
percentile was calculated for R for each of the first five events, and is shown 
on Figure 5 1. 

Using Equation C13 from Appendix C, dimensional analysis for fragment 
trajectories, a new parameter was calculated: 

g; = A f /w f 2/ ’ (73) 


Plotting the geometric mean distance R of events 1, 2, 6, 7 and 8 ver- 
sus G' (Figure 54 demonstrated a good correlation between these two variables), 
the scaling of W by taking the 2/3 root brought the distance R into a good 
linear correlation (on log-log paper) with G^ . A second prediction equation, 


R ■ 9. 864 (A/W 


2/3^0.78 


(74) 


results in a correlation of 0. 97. Thus, within the limits of th • accuracy 
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allowed by the fragment data available from the»e A 5 events, the prediction 
equation could be used to predict mean dietance R for a fragment of given 
maximum projected area A, and weight W* 

Data sources and analysis are given in the body of the report. The 
damage potential from fragments as a function of range for exploding mi^les, 
can be as important a consideration as the damage potential due to the blast 
wave. To obtain useful data such as probability of damage as a function of 
range for various missile types and "targets” of interest, 8 ^*?; ter ! on 
damage to the particular target must be accepted and the probability ^arrival 
at any given range of a fragment meeting that criterion must be calculat . 

The former can perhaps be accepted from bomb fragment damage studies, 
although fragments from liquid propellant explosions are larger and s ow . 

The latter must be obtained from empirical equations developed to define the 
distribution of fragment characteristics at the instant the missile explosion 
occurs and from suitable solutions of the trajectory equations. Using these 
techniques, a "mapping" of the fragments from missile to target and P r °ba b - 
ility of damage to a specific target as a function of range may be calcul ^ ed ‘ 

We suggest an expanded data base relating to distribution of fragment charac- 
teristics from missile explosions. The supplementation of b omb * ' 

data to cover larger and more massive fragments will probably be neces y 
before accurate damage probability maps can be made. 


- I 
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VII. CONCLUSIONS 


This report presents the results of an extensive study of fragmentation 
effects of liquid propellant rockets or vessels which explode after accidental 
mixing of the propellants. The work was entirely investigative and analytic, 
and included no experiments. 

Some specific conclusions are as follows: 

(1) Accident reports and results of tests simulating accidental 
explosions provided a significant source of data on fragmenta- 
tion effects of exploding liquid pxopellant vessels. 

(2) Data on initial fragment velocities, masses and shapes of 
fragments, and ranges of fragments showed wide statistical 
variations. When the data were fitted to statistical distribution 
functions, however, good statistical fits were obtained for these 
parameters or combinations of these parameters. In some 
instances, the combinations of parameters were chosen by 
dimensional analysis. 

(3) Methods were developed or adopted for prediction of blast 
wave properties, initial velocities of fragments, and fragment 
range. These methods were compared with data and appeared 
to give reasonably good correlation. 

(4) Some approximate methods for estimating damage to various 
"targets" by impacting fragments are presented. These could 
not, however, be correlated with data obtained during the data 
retrieval part of our work because no data were available 

on terminal (or impact) velocities of fragments from bursting 
propellant vessels. 

(5) The analyses and empirical fits to data included in this report 
do allow prediction of blast yield and subsequent fragmentation 
patterns and velocities for the common propellant combinations 
over a wide spectrum of type of accident, weight of propellant, 
and time of ignition. 

We believe that the work reported here constitutes the first relatively 
comprehensive study of fragmentation effects from exploding liquid propellant 
vessels. As noted above, predictions can be made of some of these effects 
using results from this report. But, there are limitations imposed by limita- 
tions in the fragmentation data - which has, after all, been retrieved from 
sources in which study of fragmentation effects was secondary or even entire- 
ly incidental. There seems to be little doubt that one could better test the 
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prediction methods presented here if he were able to design and conduct exper- 
iments with the specific purpose of observing and measuring fragmentation 

effects. 
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VIII. recommendations 

The results of this study can, we believe, form the basis for develop- 
ment of relatively simple methods of assessing hazards to people and damage 
to facilities from the impact of fragment from liquid propellant explosions. 

But, the study does not in itself generate such methods because its primary 
aim was the collection and analysis of fragmentation data from past tests and 
accidents. We have pointed out previously that none of the data collected 
came from experiments designed to obtain initial or terminal fragment effects. 
There is also an almost complete lack of terminal ballistics effects uata or 
methods for assessing hazards of the relatively large, odd- shaped, low- vel- 
ocity fragments which predominate in liquid propellant explosions. The 
authors therefore feel that, although simplified methods for overall estimation 
of fragment hazards can be generated, some additional experimental work 
would be very desirable to validate these methods. 

We recommend the following studies: 

(1) sing the initial velocity prediction methods developed here, 

the statistical fits to data on fragment mass and shape, and 
exercising the trajectory equations of Section IV, expected 
variations in impact velocities for a spectrum of fragments 
could be computed for a number of hypothetical explosions. 
These impact conditions could then be overlaid on expected 
densities of human observers or bystanders and on nearby 
structures to estimate fragment hits. Using the very much 
oversimplified assumption that a fragment hit on a person 
is a kill, or that a fragment hit on a structure causes some 
specified dollar damage, each "scenario" could then be 
assessed for fragment kills and property damage. In this 
way, a simplified damage assessment model could be gen- 
erated, based on the work in this report. It is recommended 
that this be done, but that the user be strongly cautioned that 
the model could be considerably refined and improved, if 
better data were available in certain areas. 

(2) A careful set of small-scale experiments be designed to obtain 
more accurate initial velocity data for tests simulating several 
types of accidents. Vessels or nearby objects should be 
carefully designed to produce fragments of known geometry 
and size, and their motions followed with at least two high- 
speed cameras. Such data should serve to Improve the data 
presented in this report, with much more "two-view" data 
yielding accurate trajectories. 
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(3) As a part of the initial velocity experiments, accurate 
missile maps should be made to determine fragment impact 
locations. The pre -formed fragments should be carefully 
marked so that impact locations of many of the fragments 
could be correlated with initial velocities and launch angles. 

(4) Studies of the terminal ballistics effects of impacts of relative- 
ly large, slow fragments of irregular shape would be very 
desirable for animals (to make estimates for humans), and 
for a variety of structures. Cost might be prohibitive. Some 
analytic studies could well provide partial answers, however. 

The order in which the above recommendations are listed is also our 
suggested order of priority. The first study is purely analytical, but yields 
an approximate hazards assessment method. The other three are primarily 
experimental, and serve to generate data which should refine the method. 
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APPENDIX A 

MODEL ANALYSIS FOR MIXING OF LIQUID ROCKET PROPELLANTS 


In the Project FYRO studies, a basic assumption was that, for any par- 
ticular type of simulated accident, the time of ignition to produce maximum 
blast yield could be scaled as 


t 



(Al) 


where T is scaled time, t is time of ignition delay after missile failure, and 
W is total weight of propellant in the missile. It is not clear that this is the 
proper scaling, and no model analysis is presented in the PYRO final repor s 
to justify such scaling. There are instead statements that the experimental 
data appear to agree with this scaling, but the inherent scatter in blast yields 
makes this conclusion doubtful. We thought that a model analysis snouldbe 
conducted to ascertain, if possible, how ignition times for maximum yield 
should scale. 


From both the Project PYRO work and Dr. Farber's work, it seems 
clear that intimate mixing of fuel and oxidizer is needed to maximize the blast 
yield for a given type of accident. The time for ignition must be great enough 
to allow as much mixing as possible, but not so great that all of the most 
volatile liquid can have evaporated. The physics of the Processes which occur 
on mixing of LH, and LO z , RP- 1 and LO,, and LH 2 /RP- l/LO z have been 
studied most exhaustively by Farber and his associates. They have, however, 
considered the dynamics of impact or other accident leading to the mixing in 
only a cursory manner. To conduct a model analysis, we should be able to 
list the physical parameters affecting this problem, and their dimensions, 
which is the prerequisite to conduct of the analysis. 


A number of the physical processes which occur are a function of the 
particular propellant combinations, and conditions just prior to an accident, 
rather than the details of the accident. These processes should be governed 
bv parameters wnich will be common to all possible types of accidents, and will 
thus be considered fi’-t. Those propellants and oxidizers which are cryogenic 
will be essentially at their boiling temperatures, while fuel such as RP-1 will 
be at or near the temperature of the ambient air. Once the propellants start 
to mix, violent boiling of the colder liquids will occur, while the warmer 
ones will be chilled and perhaps frozen. The processes will involve conductive 
heat transfer, convective heat transfer, and eventually radiative transfer to 
the ambient atmosphere. Gravity is perhaps important in convective processes 
and in rate of upward migration of bubbles formed during boiling. Latent heats 
of fusion and boiling are obviously important, as are temperature differences 
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and gradients. Parameters affecting these processes, and their dimensions 
a M, L, T, 0 (mass, length, time, temperature) system are listed m Ta e 


in 

1 . 


To complete the list of parameters for model analyses, we must con- 
sider specific types of accidents, and also add other parameters known ^to 
affect the blast yield. Consider first the case identified by PYRO as CBM 
(confinement-by-the-missile). This case is shown schematically in Figure A . 




FIGURE Al. 

SCHEMATIC OF CBM CASE 


A rupture is assumed to occur in the common bulkhead, idealized as a 
circular opening of diameter D Q . Oxidizer then spills through this opening 
under the effects of gravity, and mixes with the fuel. Geometry of the tankage, 
ruptured area, distance for oxidizer to fall, static head to force oxidizer 
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through the opening, etc., all seem important. These can be characterized 
by the tank dimensions shown, ullage volumes, initial masses of fluids, 
gravity, and a generic length l indicating location within the tankage. 
(Velocities of impact of fluids are important, but these are fixed once the 
other parameters described above are determined.) The specific heat of 
combustion (or explosion) for propellants mixed in stoiciometric ratios, to- 
gether with fraction mixed at time of ignition, should determine total energy 
driving the blast wave. Adding these parameters to Table Al, we have in 
Table A2 the total list of 20 parameters for this type of accident. These 
parameters will yield 20 - 4 = 16 dimensionless groups, when the methods of 
dimensional analysis are applied. One such set is given in Table A3. 

The model law in Table A3 can be used to express any one of the 
dimensionless groups (tt terms) as a function of the remaining ones. Because 
we are interested in time scaling, which is contained in tt^j, w e would prob- 
ably express rrjj as a function of through n ig, and n i2 tlirou 8 l1 
The law can alsobe used to fix interrelations between scale factors for physical 
quantities. In its present form. Table A3 is too general to provide much 
guidance. It states that all of the dimensionless groups must be maintained 
constant between model and prototype tests for accurate scaling. Let us con- 
sider, however, the practical limitations of testing and the effects of these 
limitations. 

Two of the physical parameters in Table A2 are constants of nature^ for 
testing on earth, and cannot be altered. Scale factors for these quantities , 
the Stefan-Boltzmann constant a and earth's gravity g, are therefore unity, 
i. e. , 



Also, we wish to employ the same propellants under the same initial tempera' 
tures and atmospheric conditions (or, at least, this is what was done in 
Project PYRO). This renders a number of other scale factors unity, namely: 


X 







p J i 




1 


(A3) 


These limitations cause several tt terms to be identically satisfied, namely, 
tt j and tt 12 . Furthermore, by making the mod si geometrically similar to 
the 1 prototype in all respects, a number of other tt terms will be satisfied. 


♦ A scale factor X is defined as the ratio of a quantity in the model to the 
corresponding quantity in the prototype. 


TABLE A1 - PARAMETERS AFFECTING HEAT TRANSFER OF 
MIXING LIQUID PROPELLANTS 


Parameter 


Symbol 


Dimensions 


Ambient Air Temperature 

$ 

Initial Temps, of Liquids 


Temp, in Mixture 


Heats of Fusion 


Heats of Boiling 


lV 2 

L V 2 


Masse: 


Thermal Conductivities 


Kinematic Viscosities 


MLT -3 0‘ 1 
L 2 T~ 1 


Volumetric Heat Capacities 


‘"Vi 


Time 


Gravity 

Convective Heat Transfer 
Coefficient 

Stefan-Boltzmann Constant 


-3 „-l 

MT 0 


Mr 3 o' 4 


Subscript i denotes a number of similar properties 







Parameter 

Symbol 

Dimensions 

Ambient Air Temperature 

9 

a 

9 

Initial Temps, of Liquids 

0. 

1 

9 

Temp, in Mixture 

6 

m 

9 

Heats of Fusion 

H fi 

L V 2 

Heats of Boiling 

«bi 

l 2 t" 2 

Heat of Explosion 

H 

e 

l 2 t" 2 

Masses 

M. 

l 

M 

Thermal Conductivities 

K. 

l 

MLT" 3 9 1 

Kinematic Viscosities 

V 

i 

l 2 t _1 

Volumetric Heat Capacities 

<pc P , i 

ML" 1 T" 2 9~ 1 

Time 

t 

T 

Gravity 

g 

lt" 2 

Convective Heat Transfer 

h 

MT" 3 9" 1 

Coefficiency 



Stefan- Boltzmann Constant 

a 

-3 -4 

MT 9 

Tank Diameter 

D 

L 

Tank Length 

L 

L 

Opening Diameter 

D 

o 

L 

Generic Length 

l 

L 

Ullage Volumes 

V . 
ui 

X. 5 

Total Mass of Propellant 

M t 

M 




f, 

< 


TABLE A3 - DIMENSIONLESS GROUPS FOR PROPELLANT 
MIXING FOR CBM ACCIDENT 
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These are tt^, tt^, tt^ and ^jq. Employing these restrictions, the re- 

maining n terms require the following interrelations between scale factors. 


TT 

2 



1 


TT , 

3 

V TT 5 

TT 

— - X = X 1 

11 

t L 

TT 

II 

t 

13 

h L 

TT 

— .X 3/2 = 

14 

L 

TT 

1 

:v 

ii 

15 

L 

TT 

— x 1/2 = 

16 

L 


1 = X , 


1/2 


(A4, 


A quick examination will show that these interrelations can only be satisfied if 
all scale factors, including the length scale factor X^, are unity . That is, 
no sub-scale model is possible which correctly scales all of the phenomena 
which were initially assumed to be important! This impasse is not unusual 
in attempting to model complex phenomena. One must now consider those 
effects which will hopefully be of minor importance, and let them deliberately 
go out of scale. 

For the problem we are considering, radiation heat transfer to the 
outside atmosphere can perhaps be safely neglected, because the mixture is 
confined within the tank up to the time of interest. So, let us ignore term 
. For the mixing fluids, conduction and convection are probably the 
primary modes of heat transfer, so we wish to retain Nusselt No. , be- 

cause this term is a ratio of these two effects. Inertia effects are undoubtedly 
important and should be retained. Gravity is triply important because it 
affects mixing impact conditions, convection, and bubble migration. Whether 
viscosity is important is not certain. We will assume that it is not, and ignore 
tt 14 . Heats of fusion, boiling and explosion appear to be important, so we 
should retain TT 3> v 4 and n _. This reduced form of the law can be written: 
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geom. similarity, 


ti*t 



\m 3/2 

V M t g 


l i / 


(A5) 


But, we still have problems! The second term ?n the bracket in Eq. (A5) re- 
quires that ^ = 1, which negates sub-scale testing. If we relax this require- 
ment, the last term in the bracket makes 

X 1/2 = X (A6) 

L M t 


But, = if identical materials are used, and (A6) is only satisfied if 

Xj^ = 1 , *Xo obtain any model law which c oes allow small-scale testing, we 
must throw out all effects except inertia, "he at conduction, and heat convection. 
This very much reduced law is 



(A 7) 


It says that, if geometric similarity is maintained, a model will have similar 
temperature distributions during the mixing process at similar locations and 
similar scaled times, provided the film coefficient h scales inversely as the 
length scale factor. Recalling that ^ = 1 and X L = X Mfc ' 3 , the dimension- 

less time parameter (t g ' 2 / L ) requires that time scales as 


7=<t/M t 1/6 ) < A8 > 

'«’here T is not dimensionless, but uniquely determines the dimensionless time 
parameter mentioned before. 

This conclusion regarding time scaling is, of course, dependent on the 
series of assumptions used to r ji -ct the basic me. .el law. One critical 
assumption was that gravity effects were important and must not be allowed 
to go out of scale. Although the physics of the mixing process seem to be 
strongly dependent on gravity, let us examine the consequences of letting 
gravity go out of scale. We will make the same assumptions as before regard- 
ing use of the same fluids at the same initial temperatures, and assume geome 
trie similarity. Let us also modify the n terms somewhat by combining and 
substitution. We can, for example, square TT U and multiply it by "through 




tt to give new terms which can be substituted. In a similar manner, g can be 
eliminated from all terms but "jj. The resulting modified terms are: 


"4 = «bl tZ/hZ 

TT ' = H t 2 /L 2 
5 c 

' = V. t /L 2 
14 1 

" ' « OQ L /k. t 6 

l!> a 1 

TT, ' =k. 0 t 3 /M L 
16 1 a t 


(A9) 


These terms dictate the following interrelations between scale factors 



(A10) 


As before, only the assumption that all scale factors are unity will satisfy all 
of these terms. But, heats of fusion, boiling and explosion will be properly 
scaled by choosing X = X So, by letting gravity, viscosity, radiation, and 
ratio of conduction to inertia go out of rcale, we can generate another lestrict 
ed scaling law. 

(H. t 2 /L 2 ) = f 3 ^geom. similarity, (0 m /0 a ), j^-j < Al 

Time scaling related to total mass of propellants scales in this law as 


T- <t/M t 1/3 ) 

This can be seen to be the scaling for time to produce maximum yield which 
was assumed in Project PYRO. 
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Which time scaling is correct, Eq. (A8) or Eq. (AlZ)? If all phenomena 
are to be properly scaled, neither is strictly correct. If gravity effects are 
truly important, Eq. (A8) is more nearly correct. If scaling of heats of fusion, 
boiling, etc. is more important than scaling gravity effects, then Eq. (A 2) 
is more nearly correct. To test either hypothesis, sufficient test data must 
exist to fit data scaled in either manner over a relatively large range of total 
propellant weights for any given combination of fual and oxidizer. Even then, 
it may be difficult to determine which scaling law to use because the propellant 
mass is raised to a small fractional power (1/3 or 1/6) in either case, and the 
dependence of the scaled time on propellant mass is therefore quite weak. 

We reviewed the PYRO test results to ascertain whether those data sub- 
stantiated correlation of a specific law for scaling of ignition time with blast 
yield. 


The PYRO tests indicated that, for hypergolic propellants, ignition 
time is unimportant because ignition occurs immediately on contact of the 
fuel and oxidizer. Blast yield did depend on impact velocity for high velocity 
impact tests of hypergolics. 

For non- hypergolic propellants, one would expect the blast yield to be 
a function of type of simulated accident, type of propellant, impact velocity, 
etc. , in addition to time of ignition. Blast yield is expressed in the PYRO 
reports as a percent of an equivalent weight of TNT, based on the measured 
terminal yield from each experiment. This procedure essentially normalizes 
the results for all tests with respect to mass of propellant so that data on 
yields versus ignition times or scaled ignition times can be easily compared. 
We will make such comparisons for the propellant combinations L0 2 /RP-1 
and LO_/LH_, and for the CBM and CBGS type of test. 

In Figure A2, we see blast yields Y plotted as a function of ignition 
time t for all LCL/RP-1 CBM tests*. Different symbols are used fo- differ- 
ent masses of propellant. It is immediately evident that the data scatter is 
large, but the data do indicate an increase in yield up to 150-200 msec, and 
then a substantial decrease for longer ignition times. The data are plotted as 
a function of time, rather than time scaled by division by either or 

M 1/6. Re plotting the data with time scaled in e ilher of these manners 


* For physical reasons, finite blast yield at zero ignition time is paradoxical. 
If propellants are not given time to mix, no explosion should be possible. 
Data points for zero ignition time probably represent a small, but non-zero, 
igr.ition time. 
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FIGURE AZ. YIELD VS. IGNITION 



produces no less data scatter than is already present. Figure A3 show s a 
similar plot for LO /LH CBM tests. Again, the scatter is considera e. 
but a trend to increasing 1 yield with increasing time is e ' w ltl ? , n ° apP * 
ent decrease at long times. One isolated data point with a high yield 
indicate peaking at near 200 msec, but this is inconclusive. Again, replott g 
on scaled time based doesn't improve the scatter. 

Data for all LO,/RP- 1 CBGS tests are plotted in Figure A4. It is evi- 
dent that impact velocity for this propellant and type of test does have a sig- 
nificant effect on blast yield, and that yield increases for a m«imum ; and 
then decreases. The same trends are evident in the data for L0 2 ' LH 2 C 
tests plotted in Figure A5. Some type of normalization of yield versus ™pac. 
velocity would undoubtedly reduce the scatter evident in these two fl 8 u ^ * 
the spread would still be great. As for the CBM tests, there is no evidence 
that scaling ignition time would decrease the data scatter. 

We draw different conclusions from this study of the PYRO blast data 
than do the authors of that study. They concluded that, (or all bu^ the LO ? / 
RP-1 CBM tests, explosive yields scaled as a function t - (t/M^ • 

conclude, on the other hand, that the data do not verify any P^ rtlcular 3Cal ^ g 
for ignition time. The blast yield is certainly a function ot ignition time an 
that function can perhaps be estimated from the PYRO data for specific pro- 
pellants and types of test. Direct plots of yield versus ignition time appear to 
^ as valid a s^ plots versus (t/M t 1/3 ), (t/M t ‘ /6 ), or perhaps any other scaling 
involving a small fractional power of M t . 





FIGURE A3 . YIELD VS. IGNITION TIME LO /LH CBM TESTS 
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APPENDIX C 


MODEL ANALYSES FOR FRAGMENT VELOCITIES, RANGE, ETC 
FOR BURSTING LIQUID PROPELLANT VESSELS 


In an attempt to facilitate comparison* of the data which have been un 
earthed during this study, let us conduct several limited model analyses re- 
lating to various aspects of the problem. To avoid too complex model laws, 
we will divide the problems into several subproblems. 


Initial Velocities for CBM Case 


Consider first the bursting and acceleration of vessel fragments by an 
internal explosion (CBM case). A list of physical parameters appropriate to 
this case and part of the problem are given in Table Cl, with dimensions in a 
Force-Length-Time (FLT) system. The twenty-one parameters in Table Cl 
can be combined into eighteen dimensionless groups (Pi terms) by 'he methods 
of dimensional analysis. One such grouping is given in Table C2. based on 
mass of propellant, heat of explosion i f propellant, and a characteristic engt 1 
as repeating variables. Table C2 includes descriptions of the pi terms. 


The dependent variables in Table C2 are essentially M jy ' J4’ ^ \t> 


md 

n Terms n«^Srough n are mass or density ratios, and terms through 
n , require geMetric similarity. Term n l2 dictates scaling of energy re- 
lease, and is equivalent to fractional blast yield y. Term can be consid- 
ered as a ratio of vessel burst pressure p v to the pressure which would be 
generated by contained reaction of all of the mas* of propellant. Finally. ^ 18 
requires identical equation of state for the reaction productions of the explod- 
ing propellant. 


The model law in Table C2 can be expressed in functional form as 


[(U f /ny z >. (pb 3 /M t ll ( ) , (my 2 /:.), nj 

3 ^ Density ratios, Geom. s rr . , (E/M^ll^). 

(p L 3 /M H ) , k 1 (Cl) 

If we confine our comparisons to tests with geometrically slmll ir tanks made 
of the same materials and filled with the same propellant mixtures then the 
density ratio terms and terms describing geometric similarity will all be the 




TABLE Cl - PHYSICAL PARAMETERS FOR VESSEL BURST 
BY INTERNAL EXPLOSION 


Dimensions 


ft 2 l'* 

ftV 

ft 2 l*' 


. 2_-2 
L T 


ft 2 l * 4 

ft 2 l ' 4 



Parameter 


Total mm of vessel 

Total mm of propellant 

Reactive mm of propellant 

Energy releaee daring ex- 
ploeion (blast yield) 

Heat of explosion of pro- 
pellant mixture 

Tank diameter 
Tank length 

Diameter of internal bulk- 
head rupture area 

Thickness of vessel material 

Density of propellant 

Density of vessel material 

Number of fragments 

Mass of individual fragments 

Mean presented area of indi- 
vidual fragments 

Velocity of individual frag 
ment 

Internal volume of tank 

Ullage volume of tank 

Ratio of specific heats 'o* 
explosion products 

Pressure within tank during 
explosion 

Time 

Buret pressure of vessel 





TABLE Cl - DIMENSIONLESS GROUPS FOR VESSEL BURST 
BY INTERNAL EXPLOSION 


Number 


Description 


M / M 
v t 

M / M, 
r t 

M f /M t 

p L 3 /M t 
P t 

p L 3 /M t 
v t 


Ma.«s or density ratios 


D /L 

o 

h /L 
v 


Geometric similsrity 


V /L J 
u 

V /L 3 

V 

E/M t H # 

U ,/H 1 '* 

f e 


pi. /MH 
r t e 

p L 3 /M II 

V t < 

.h''V 

e * 


Energy scaling 
Velocity scaling 


Pressure scaling 


Time scaling 


Number of fragments 


Ratio of specific hi ats 





V 


same independent of scale. X^ = X Mt * and k will not change. Also, the burst 
pressures p v for geometrically similar vessels of the same material are iden- 
tical (X p v = 1). Because we are using the same propellant. X He = 1. Under 
these restricted conditions, Equation Cl reduces to 

[u f , p. <t/u, n] = f (y) (c2) 

In other words. Equation (C2) says that initial fragment velocity, pressure rise 
within tank prior to rupture, time scaled proportional to length, and number of 
fragments generated should all be functions of fractional yield. 

Many other interpretations can be made of this law. If we do not make 
the simplifying assumptions above, then we rr.ay be able to compare s >me test 
results for different propellants. Terms and n J3 can be combined to form: 

"I 2 * u t M , /E (c3) 

which :tn replace in Equation (Cl). 

This term defines an equivalent dimensional group 



U 2 W t /Y 


(C 4 ) 


where W t is total weight of propellant and Y is blast yield in percent. We 
have recorded most of our fragmentation data in terms of these quantities, and 
so can plot Gi, , or its square root, versus Y (equivalent to Hj,) or against 
n. 5 , provided we use a consistent set of units. It is quite likely that k is only 
weakly dependent on type of propellant, no perhaps we can Ignore it. The 
model law does point oat that (D/L) and (D/D 0 ) are potentially Important geo- 
metric ratios, and thai mean fragment area may scale as the square of the 
length scale. 


Initial Velocities for CBGS and HV1 Cases 

In this case, the propellants spill, ml* on the ground, and then ignite 
or the missile Impacts the ground at high velocity. Solid objects and nea rhy 

~*>The symbol X denotes a scale factor, and a subscript letter denotes the 
physical parameter being scaled. 

**No,'e that these products are not now rTondimenelonal. because we have 
delated dimensional terms which do not change. 
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appurtenances are then accelerated by the resultant blast wave to some maxi- 
mum velocity ("initial" velocity). The pertinent physical parameters are some 
what different than for the CBM case. Geometry of the vessel is no longer im- 
portant, but gravity and impact velocity must be considered. A list is given in 
Table C3. Again using as repeating variables H e and L, we obtain a group 
of fifteen pi terms from the eighteen physical parameters. These are given in 
Table C4. 


Many of the terms in Table C4 can be seen to be identical to those in 
Table C2. These include the mass ratios ttj through tt,, area ratio n 4 , volume 
ratio n_, scaling for fragment velocity rt<j, energy scaling ttjq, pressure scal- 
ing tt j£, time scaling ttj 4 , and number of fragments ^ 15 . The terms which are 
new are n^, tt^, tt^, and nj., and some terms in Table C2 do not appear. 
This model law can be expressed as: 

[( U fl / H * /2 ), (pL 3 /M t H 4 ), (qL J /M t H # ), ftH 3 / 2 /U, (A {[ /L 2 ) , n] 

1 / 2 

■ f £ Density ratios, Geom. slm. , (E/M t H # ), C^, ), 

(gL/H # )] <C5) 


The term n n , if strictly adhered to during tests of the same propellant, would 
prevent comparison of testa conducted at different geometric scales, because 
X g = 1, \ H# = 1, and therefore X^ * 1. 

If we deliberately let this term go out of scale, which is equivalent to 
ignoring giavlty effects, then almost the same model law as for the CBM case 
results, with the addition of a requirement of scaling of Impact velocity from 
”q. If we assume that we cannot Ignore effect of gravity on such physical pro- 
cesses as spreading and intermingling of the fuel and oxldlser after spill, ini- 
tial trajectory of fragments, etc. , then Equation (C5) requires a changs in heat 
of explosion H, in order to change length scale. It is suggested that we simply 
try plotting fragment velocities from CBGS data versus blast yield for the same 
propellant. We know from the PYRO work that blast yield is a function of im- 
pact velocity Vj , so term may oe already accommodated by meisurlng 
blast yield. An alternative might be to plot the term 


", * v (l //7T (CM 

car sue blast yield for the same propel’ant. Correlation here would indicat*) 
that gravity effects were indeed impor .ant. 
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Wf 

M 

t 


*i * i 
FT L 

^ 2 , * l 
FT L 


lV 2 


Reactive mtn of propellant 


Mass of the solid body or 
appurtenance 

Energy release during explosion 
Heat of explosion of propellant 
Volume of propellant 
Impact velocity 
Acceleration of gravity 
Length ratios 

Characteristic length 

Number of fragments 

Mean presented area of I th solid 
body or appurtenance 

Velocity of same 
Drag coefficient of same 
Overpressure in blast wave 
Dynamic pressure in blast wave 
Time 
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TABLE C4 - DIMENSIONLESS GROUPS FOR FRAGMENT 
ACCELERATION AFTER SPILL 


Number 


Description 


M /M 4 
v t 


' /M. 


M ri /K: ( 

11 i 


A ( ./L‘ 


Mass ratios 


Geometry 


Drag coefficient 


V /L 3 
P 

1 e 

U„/H 1/2 
fi e 


E/M H 

t e 


Volume ratio 


Velocity scaling 


Energy scaling 


gL/H 


pL7MH 
t 6 

q L"*/M H 

t © 

.H ,rt /L 


Gravity scaling 


Pressure scaling 


Time scaling 


Number of fragments 




Trajectories of Fragments in Free Flight 

After the fragments have been accelerated to their maximum velocities, 
then the problem is one of exterior ballistics for each fragment. Most of the 
parameters governing the resulting trajectories differ from those in the previ- 
ous problems. They are listed in Table C5, and a possible set of pi terms is 
given in Table C 6 , 

The first four pi terms in Table C 6 are initial conditions, ^5 and n 6 
are aerodynamic coefficients, and ti 4 through n 10 specify geometric similarity. 
Term n q can also be considered as the dependent variable. Term tt l x specifies 
gravity scaling, which is essential in trajectory problems and cannot be ignored. 
Terms rr^ and are wind velocity and air density scaling, respectively. 1 ne 
model law can be written as: 

(r { /l) =f[(u f t/L), 9 f , Y ( . C Df , C Lf , (A f /L 2 ) . fx if /L) . i., 

(gt 2 /L), (U w t/U, Co a L 3 /M ( ) ] (C7) 

From our physical knowledge of exterior ballistics and the problem of frag- 
ment scatter, we can considerably reduce this function space. Range may be 
dependent on azimuth angle Y and scaled wind velocity (U w t/L), but this depen- 
dent is weak for high velocity and "chunky" fragments. Furthermore, our data 
from missile maps average ranges over all azimuths, so this dependence is 
not considered. The "lifting" characteristics of the fragments are represented 
by C L and by initial spin <*> f t. Again, these characteristics represent random 
and uncontrolled quantities which we cannot assess, so we again ignore them. 

A reduced form for Equation (C7) is then 

(R ( /I.) =i[(u f t/l.), 9 ( . C D , (a ( /L 2 ) , (x lf /L), lj. fet 2 /L), 


(0 a L 5 /M ( )] (C( 

Of these remaining parameters, the first two * nd the fifth are scaled initial 
conditions. The d*ag coefficient is a function \Aj/L ) and i- x and is there- 
fore superfluous. So, for fragments of similar geometry and s me scaled 
initial conditions, the law iurther 1 educes to 


(R ( /') > f[ (U ( t/1.) , (gt 2 /L), (o a L J /M 1 )_ 


(C 9) 
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TABLE C5 - PHYSICAL PARAMETERS FOR 
FRAGMENT TRAJECTORIES 


Symbol 

Dimensions 

Parameter 

U f 

lt" 1 

Initial linear velocity 

(JJ 

f 

t" 1 

Initial angular velocity 

C D f 

— 

Drag coefficient 

% 

— 

Lift coefficient 

A t 

L 2 

Mean presented area 

M t 

ft 2 l _1 

Mass of fragment 

9 f 

— 

Initial elevation angle 

T f 

— 

Initial azimuth angle 

X. . 
i i 

L 

Initial coordinates of 
fragment 

R f 

L 

Range of fragment 

g 

lt' 2 

Acceleration of gravity 

D 

a 

ft 2 l' 4 

Density of air 

t. 

1 

— 

Length ratios 

U 

w 

LT -1 

Wind velocity 

t 

T 

Time 

L 

L 

A characteristic length 


183 




TABLE C6 - DIMENSIONLESS GROUPS FOR 
FRAGMENT TRAJECTORIES 





requires that 


Now, because X^ = 1, 


x f - h. 


(C 10) 


The first term in Equation (C9) then requires that 


, , 1/2 _ , 

A u, l - x l' 

f 


or 


, _ ,1/2 (C 1 1) 

X U f * X L * 

i.e. . initial velocities should scale as the square root of the length scale for 
the same scaled range. The third term is automatically sat sfied because 

Xp a = 1 and X* = X Mf . 


We can combine tt^ and rr^ to form 
2/3, w 2/3 


/ 

'7 


TT„ = A p 


f "a 


/ M 


(C 1 2 ) 


Because 



= 1, 


a dimensional form of this term is 



A /W 


2/3 

f 


(C 1 3 ) 


where is weight of fragment. 

This may provide a rational grouping for plotting versus scaled range. 
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COMPUTER PROGRAM ENTITLED /W2 / IN FORTRAN IV 
Function: This program computes initial fragment veloc ity. 

Given the following input data: 

A) Characteristics of some explosive fuel-oxidizer gas mixture at moment 
t - 0 of detonation. 

(C API) Ratio of specific heats, K 

(A<J) Speed of sound in medium in in. / sec 

(P0) Initial pressure in psi 

B) Characteristics of containing vessel 

(rr) Internal radius of vessel + unburned fuel in inches 

(TM) Mass of vessel + unburned fuel in lb/ sec 2 /in. 

(FN) No. of fragments 

C) Calculator y requirements 

(AH) Time interval of each calculation in seconds 

(XMAX) Maximum time to last calculation in seconds 

Variables: The definition and units of the vari ies in this program are give- 
in the following table. 

TABLE (D-l) 


Program 

Variable 

Variable 

Definition 

Units 

FF 

F 

projected fragment area 

2 

in. 

C API 

K 

specific heat ratio for explosive 
mixture 

none 

AO 

a 

o 

speed of sound in explosive pr^xiucts 
mixture 

in/ sec 
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Program 

Variable 


Variable 


Definition 


Units 


PO 

P 

initial pressure after explosion 

psi 


oc 



FN 

n 

number of fragments 

none 

RR 

R 

radius of explosive products mixture 

in. 

TM 

M t 

mass of shell + unexploded fuel 

*> 

lb-in / sec 

FK 


coefficient 

none 

AH 


time interval 

sec 

XMAX 


maximum time 

sec 

FN 1 


if ^ 1 displays T-NORM (>» G 1 , G 

none 

FN2 


if < 1 displays normal pressure + time 

no le 

FN3 


if < 1 calculates maximum range 

none 

G 1 

U f 

distance to initial velocity 

in. 

G2 


initial fragment velocity 

ft/ sec 

G3 


initial fragment acceleration 

in/ sec 

G4 


final explosive product mixture pressure psi 

T 1 


time to initial velocity 

sec 

JJ 


counting variables 

none 

PI 


the constant n 

none 

CAP2 


the quantity (1- K )/ K 

none 

C AP3 


the quantity - 1/ H 

none 

C AP4 


the quantity (3 H -1)/2 K 

none 

XX 


displacement normalization coefficient 



Units? 


Program 

Variable 

THETA 

Al 
B 1 
C(J 

X 

Y(2) 

Y(3) 

Y( 1) 

NA 

F( 1), F(2), 
F(3) 

TT 

PS 


Variable 


a 

3 


Definition 


time normalization coefficient 
(see Eq. -35) 

the coefficient a (Eq. - 36) 

the coefficient ^(Eq. -36) 

normalized initial fragment displace- 
ment from center of sphere 

normalized time 

normalized velocity 

normalized pressure 

normalized fragment displacement 

number of differential equations to 
be solved 

differential equations solved (see 
Eqs. -34 a nd -36) 

normalized time 
normalized pressure 


sec 

none 

none 

none 

none 

none 

none 

none 

none 


none 

none 


The subprogram entitled (RUNGE) is described in the following: 
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FILE NAME: 
SUBROUTINE NAME: 
PURPOSE: 


METHOD: 


ARGUMENTS: 


■r fw’ 



FORTRAN IV 
RUNGE * KUTTA 


FIS 

RUNGE 


Thi. subroutine employs the Fourth Order Run f ff ^ n a ti ^‘^. 
to solve N simultaneous first-order ordinary dlff ' rent ‘*\ 
tions by calculating successive values of Y according to 
formula: 


Y fl' V i < K ! * 2 *V 2, V K 4> 


where * f(X., Yj) 


K 2 . f(x k ♦ §. y 4 * ^ > 

K 3 " R *i * 3!' 

K, « KXj ♦ h. y i * hK3 ) 


The subroutine is called by the calling program five times in 
order to approximate successive Y(I)'s: the first time to initialize, 
the second time to calculate KjU). the third time to calculate 
KVI) the fourth time to calculate K3U’ and the fifth time to ca 
cu’ate Kj(I) In addition, each time the subroutine is called, it 
calculates a new Y(I> and a new X(i) which are returned to tne 
calling program where the functions (first-order differential 
equations) are evaluated with the new X(I> and Y(I). These values 
of the function are then returned to the subroutine where thev arc 
used as K » (1 ). K 2 <1>. K 3 (I). ° r K 4«> and appropriately accumu- 
lated to obtain Y (1) in the 5 calls to the subroutine. 


The subroutine T.UNGE uses nine arguments: W, 7, F. X. H, M, 

SAVEY, Pin. K 

1. The first argument. N, represents the number of simulta- 
neous first-order ordinary differential equations to be solve . 

2 The second argument, Y. is the array name which the calling 
program uses to transmit the initial Y(I) values for each 
differential equation. Upon completion of the 5 calls to 
RUNGE, Y(I) vrlll contain the new approximated values for 
the Y l+1 (I)'a. 
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I. Tbt third •rgum*«t. F. is thv array which contains the cur- 
rant values at the differential equations calc lated by the 
main program, l.a. , F(J) contains the value of the J tn first- 
ordtr differential equation. 


4 The fourth argument, X, represents the independent variable 
which should be initialized in the main program before calling 
RUNGE. RUNGE increments X by the stepsise H. 

5. The fifth argument, H, represents the step size for X. 

6 # The sixth argument, M, indicates which of the five passes 
of the subroutine is to be executed. The main program must 
initialise this argument as 1. RUNGE then successively 
increments the variable by 1 up to 5. 

7. The sownth argument, SAVEY. i. used within RUNGE and 
must be dimensioned in the calling program to be of size N. 

8. The eighth argument, PHI. is also used internally by RUNGE 
but must be dimensioned in the calling program to be of size 
N. 


8 . 


The ninth argument. K. ia manipulated within RUNGE. K 
should be teiUO right after the call to RUNGE. in the calling 


program. 


When K*l, control should transfer to a set of code in the 
calling program which calculates new values for the rirst- 
order differential equations. F{1), with the current values 
of X and Y( I). Then RUNGE should be called again. 


When K*2, the approximation for Y(I) is completed. Values 
for the Yi+idVs are stored in Y(D at this time, and normal 
flow of the calling program should resume. 


LIMITATIONS 

AND 

COMMENTS: 


1. The calling program must dimension SAVEY and PHI. 

2. The calling program must set M*1 before calling HtJNGE. 

3. The calling program must set up the N first-order differen- 
tial equation values in an array F to be passed through to 
RUNGE when the subroutine returns with K*l. 

4. The calling program must set up separate ar rays if all X and 
Y values for the set of differential equations are to be saved, 
perhaps for plotting purposes. 
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PROGRAM /W2/ LISTING 


DIMENSION F(3),Y<3)a’M<3),W-'(3>aPS<50)aTT(S0> 

300 FORMAT <?/,20H READ IN TRAJ. ANGLE) 

3001 FORMAT (E10.3) nrutlTY) 

301 FORMAT <?/, 35H READ IN DRAG COEF. AND AIR DENSITY) 

3011 FORMAT (2F. 10.3) . 

302 FORMAT (2/,25K READ IN FRAGMENT MASS C0E1- . ) 

303 FORMAT <9/,25H CONDITIONS ON TRAJECTOR ) /CORIC FT.) 

30 A FORMAT </.21H AMBIENT AIR DENSITY- , El 0 . 3, I 7H LBS. /CUBIC F 

305 FORMAT </,16H AIR - DRAG COEF .* , E 1 0 . 3 ) 

306 FORMAT ( / 2 1 rl FRAGMENT MASS C 3 EF . = > E 1 0 . 3 ) 

307 FORMAT (/,1CH TRAJECTORY ANGLE-, h 1 0 . 3, Rh DEGREE j ) 

307 1 FORMAT (3/,9H CHECK C*a E l 0.3) 

308 FORMAT (3/,l5H MAXIMUM RANGE* , E I 0 . 3 , 9:1 F T • > PRESSURE ) 

309 FOR 1m r (2/, ASH READ IN <APPA, SOUND .-»?Zc.o, INITIAL PRE^S 

3H ‘c^Imh’iIEAU 1« 1'. IF FP.iiG’IENTS > RADIUS 0F EXPL.. 

MASS OF .ilVLL-FUEL, DISCHARGE C ,C -F . ) 

313 FORMA" ivl^i^READ IN Til- INTERVAL, MaEIMUM TIME) 

315 FORMAT Ivl^GAS CH/RACTK34 STICS, /p7H XAPPA.,E|0.3. ? 

I'H SOUND S °E E i )* , E I 0 . 5 , 7 H IU/SLC./aIOH J ' E1 ® ! 3 ' PSI ' 

2/ . 23 H mSW. C:l.iRACT-:»I5TICS. /. >M ,F S- , c 1 0 . 3 , 2 

SH MASS-,. 1 0.1. I LI .02.- >r.O. . i./I.I. . o 1 

316 FORMAT ( 2 / a 3 .3 .'i DISPLAY DY,),i ... 'hul. • /_; ' ' 

317 FORMAT (2/,33H DISPLAY MORM. PREoa.? fc. *1 NO A) 

318 FORMAT < 2 / a 3 3 , i 1AKF. Hai'JGE CALCULATION? Yc.S-1 ^” SEC • 

319 FORMAT (3/,l M FIWaL U.iLUES, / , SH r IMr. * , c. 1 0 . 3, 9rl SEC,/, 

10H DISTANCE* a r. 10. 3, AH IHSa/aIOH WELMC I TY= a E 1 0 . 3a 7;» FT /SEoa / a 
Lh ACCEL'S XT I * 10. 3, 10.1 I ,J/Si-SEC, /, 1 OH P.tta j'.- a »->•-' 

o20 ? FTRl\T mU C.i V3ACT .PISTILS >F MO" MU -IF FHaG .HITS ( 1 1 1. 'AL UED ) 

,2/a 1 0 i i r-.i Ia’a jMa 1 M": A I "-W 2lM * A / '.A ■ t >" A / ) 

321 FT HAT ( /.''•■ PRESSOR.. ' TJ.-UAL. i VE n ) a / / -4 a 7 i >-• -M, 

322 FORMA^ ( 1 OH INITIAL CO, ID I T I 0,15, / , 6 H < ( O ) * a E 1 3 . 9 a 7u ..( 0 ) = ,<-. 1 0 . 4 

• ( 0 ) * , E 1 O • •* , LI P- , i 1RM* > i’. I >.9,2/) 

JJ = 0 

'..•RITE ( 1 a 3 if) 

READ (1,310) LAP l » mOa 2 '■ 

WRITE (1,311) 

MEAD (0,312) FNaRRaTMaFK 
V:iITE <• 1,313) 

RF.OD (0,319) AH,XMAX 
WRITE (1,316) 

MEAD (0,3001) FU1 
WHITE ( 1 a 3 1 7 ) 

READ (0,3001) FN2 
WrUTfe (1,313) 

READ <0, 300 1 ) FN3 

WRITE (1,315) CAPIaAv),? i, (T.TMaFN 
Pi- J.191 5926535 

FF«A»0*P1*< RH**2'«00 )*<<I/FN)-CI /FN**2 . 0 ) ) 

XX-1W*<A0**2.O>*<2.O/<CAP1-I >>/<Fr*PO> 

THE TA* TM -AM * < < a . 0 / < C AP I - 1 . o > > *•0 . 5 > / < FF «■ P w > 

CAP2»( 1 .OO-CAPI )/CAPI 
CAP3»- 1 .0/1 AP2 


'CAP4* ( 3 .0*CAP1 - < 1 .-»> )/(2.0*CAPl ) 

A1 = 3 • 0 0 *FK*CAP 1 *< (2.0/ (CAP I ♦ ! • ))>**< ( Ch 


A I ♦ l . 0 > / ( 


0 * ( CaP 1 - 1 .0) > ) )*( <2. 0/ 




0 ) ) 

C£=RR/XX 

x=o.o 


y ( i )=cei 

Y<2)*0.0 
Y(3)*l .0 

WRITE ( 1,322) X,Y<1 ),Y(2),Y(3) 


NA=*3 

F<2)-F^Y(3>*(<I.00-(Y<2)**?.Q)*<VC3>**CAP2>>**CAP3: g*caP 1 *< Y( 2 )*Y 

F(3)*((Y(1 >**-3.0)*CY<3)**CA?4)*<Al*Bl-Al*( ■ ( 1 )**2.0) ) )-3.0 CA 

3) / Y ( 1 ) ) 

IF (FNl-1.0) 200,200,30 

200 '•:rite o,320) 

30 C \LL U'NGE (NA,X,Y,F, AH, KA,’U , W2 ) 

IF (XA- 1 ) AO, 50, AO 

F°2)»F ): )*Y(3>*((l.00-(Y(2)**2.0)*(Y(3)**CAP2))**CAP3) 

F('J)« ( (Y(I >**-3.0)*CY<3)**CAPA>*(AI*BI-Al <(Y(1 >**2 .0 ) > )- 3. 3*CAP1 < Y <-> 

<3)/Y( 1 > ) 

,33 T0 30 

AO IF (rN I - l .0 >AS, A5, 201 
45 WRITE (1,312) X,YC1 ).Y(2),F(2) 

201 CONTINUE 


JJ a JJ + 1 
TT ( JJ )*X 
P5(JJ)“Y(3) 

IF (X-XMAX) 41,10,10 
141 CONTINUE 
G0 T0 30 
10 CONTINUE 

IF (FN2-1.0) 202,202,203 

202 WRITE (1,321) 

WRITE (1,3011) (TT ( I ),PS( I ), 1 = I , JJ ) 

203 CONTINUE 
T1 =THE f \*X 
G1»XX*Y( 1 )-XX*C0 
G2» XX /THETA* (Y( 2 ) ) 

G2*G2/ 12.) 

G3« (XX/ (T,IETm)**2.0)*F(2) 

G4<*P0 *Y ( 3 ) 

WHITE < 1 ,319)T1 ,GI,G2,G3,G4 
IF (FN3-I.0) 204,204,205 

204 CONTINUE 

205 CONTINUE 
END 
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A SAMPLE RUN OF /W2/ 


p. 'N 

> IM <AP?A, 3 J’J -U SPEED, INITIAL PAESSJAE 
1 .4, 1 352 > . , 6000. 


?\r A: > 
CV T . 

1 J . . 


J # -IF ri^A ;.‘lE:;Toi KaDIUS EXP--.* IASS 3 r ^ar.U-F'Jtti 

• j • ‘4 ^4 5 ) 1 • ) 


•V.v II TME IJTE.-.'AL, 1 AX 1'1’J 1 'ME 
r 3 . < r * * *5 » *i • li t* — ft 2 


..I ;.v »v ’.iYmai i«. va:..? Y'S*1 :»?o 


•;I >:,.\Y “JIRM. PIES..? YES- 1 - 

1 


U.r. -I * CALC’JLATI* i? YES* I 2 


IAS C.ia..aCTc.*IS7ICS 

KAPPA* .kltf'H S0UNO SPEED 



.603 >0A PSI 

VESSEL 0 HA.'LACT Ei - 1ST I CS 

:iADI'»3 = 

270E+32 IMS 

. -IAS 3= 

; 1 . ■ p frm. 

I'lEM T 3 " . 1 

10S+M3 

IN I T I il. C ‘ 

'LIT I 'MS 


:< < • j > - . 0 >J 

, iC * . J 1 '(0) 

= . Ao 1 7 - 

C.IA'.ACTE-U 

STIC. !JF N0TIOM Jr ' 

r - ‘tip. i 


Cl 1 

• r * >0E- 0 2 

. /. -3C.-0I 

.39 ; ;E*0 ) 

. IOOE-Oi 

•617E-01 

.543E+00 

. 1 ^ -0 1 

.546? -01 

.600Ef0‘> 


. S76E-01 

.6 >3E+0 ) 

. 2 r ”> ' ) 

* 50} 0-0 1 

.6 32E+v)0 

. 3;# - 

63/E-OI 

• 636E+03 

. 3 C 0F- i 

.671E-01 

.637E^00 

.4D0F-01 

. 703E-0 1 

.638E+0J 

./I50E-0I* 

.735E-01 

. 6 38E+00 

. SOOE-O 1 

.767E-OI 

. 638E.00 


. 13520E'.. IJ/SEC 


• 4m5e+jo s • — sec • Su • / i t i , 


)i 5*0)* . lOOOE+Ol P-M.JR.1* 


( •' v.<MAl. I ZED ) 


. i6 7 K+'iW 
. I 661 ♦ , 'r! 
. 6821. +0 I 
. P o 2 t. ♦ 0 1 

.11 ".♦01 
. /|0 ' V ♦oo 

. 1 24E+00 
.'^751.-01 
.297E-02 
. 144S-04 


I, I SC i*. A ’ioE 


. I OOOf.*'! 1 



PRESSURE (NORMAL 1 2 ED ) 


T*N3RM 

P-NORM 

. 500E-02 

.992E+09 

. 100E-0I 

.70Sr.t-93 

. is6s-oi 

. 5 1 6 E 1 0 0 

.'’90E-0 » 

. 3 S7E+00 

. 250Z-3 1 

. 258E + 99 

. 3C30E-0 1 

. 1 3 1 L+O.. 

. 250E.-9 1 

.. 127_i-00 

. '400E-Q1 

.907E-01 

.A50E-0I 

.*522-01 

.500F-01 

. I74Z-01 


FINAL VALVES 
TIME* . 927E-03 SEC 
DISTANCE* . I 69Ei02 IN'S 
VELOCITY* . 1 6 l Ei-04 FT/SEC 
ACCELERATION* ,235Ei02 IH/30-SEC 
PRESSURE* . 379E+03 PS I 


i 


l 
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COMPUTER PROGRAM ENTITLED /ROOT/ IN FORTRAN IV 

Function: This program computes the root to the equation 

r ~ / « \ 1/ K T 


v r p Kr l 

oo o_ I 2. 

fTT" p ’ \ p 

oo X oo / 


(see Section II. C) for the following input data: 


(EN) 

(VO) 


Energy of explosion of reactants ft-lb 


Volume of reactants 


Variables: The definition and units of the variables in this program are 
given in the following table. ^ 


Program 

Variable Variable 


A, B, C 


Definition 

specific heats ratio 
ambient pressure 
coefficients of the polynomial 

■ - - - >'K C =0 


A P - B P„ 
o o 


Energy of explosion of reactants 
volume of reactants 
initial guess solution 
maximum number of iterations 
desired root 

^alue of the equation at X 

value of the equation derivative at X 


1 1 ■ - 

The subprogram entitled /ROOTN/ is described in the following: 
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Newton’s Root Finding 


This subroutine finds a root of an arbitrary differentiable equation F(X)— 0 using 
the Newton- Raphson iteration method. 


CALLING SEQUENCE: CALL ROOTN(X ( F,FP.XS,E,M,lFL) 


Input: 


XS = initial estimate of root 
E = error tolerance 

M = maximum number of iterations allowed 


Output: 


X = value of root 

F = function value at X 

FP = derivative of function calue at X 
I F L = error flag 

0 - if normal 

1 - if no convergence in M iterations 

2 - if derivative equals zero 


METHOD: 


Given a function f(x), find a root of f(x)=0 using an initial estimate x s . 
The iteration algorithm used is Newton- Raphson: 


where 


f '(Xj) is the derivative of f(x) evaluated at xj. 


The procedure has converged if: 
f(x j+1 ) < 100E 

and if either: 


x i+ 1 -x i 

1 i < E when I Xj + 1 I > 1 

x i + 1 


or: 

I Xj 4 j - Xj i < E when I Xj + y I < 1 
where E is the user supplied error toerlance. 
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COMMENT 


Refer r»ni e 


* 


i 

i 


The user suppliesa SUBROUTINE subprogram XNFUN(X.F,FP. 
which computes f(x) and f *<x). 


Hildebrand. F B. \nu<uUjcuan to Nu merical Analysts McGraw Hill 
New York . 1966 PP 447 450 


PROGRAM /ROOT/ 


C0NMON /ABC/ A, B, C, D 


1 Jr" 3RMAT 
R format 
3 FORMAT 
A V ORMAT 

5 FORMAT 

6 FIRM AT 

7 F 3 AM AT 


(2/,15H READ IN ENERGY) 

(El 2.5) 

(2/,15H READ IN VOLUME) 

(2/, AH EN=, E 1 2. 5 , AX, AH V3=,E12.5) 
(2/, 1 1 H READ IN X3) 

(2/, 10H READ IN E) 

(2/, 1 OR READ IN M) 


6 F 1RMAT (12) ' 

9 F 1RMAT (.2/ j 3H X=, El 2. 5, AX, 3H F= , E 12 . 5, AX, AH FP 
IN F3RMAT (2/,llH CALC. A3 ?) 

11 FORMAT (AH YES=1,aX, 5K N3-2) 

! FORMAT (2/, AH A3 = , El 2. 5 ) 


XX = 1 .2000 
P= 1 A . 7 
A= 1 .0/P 

0= ( 1 ,n/D)**( 1 ,0/XK) 

nr i n: (i,i) 

REAM (0,2) EM 
NRITE (1,3) 

READ (0,2) HO 


V u r 

( i 


EM, vr? 


c= ( i 

*( ’< 

7- i 

. 0 ) ) / ( P * 

V3* 1 AA . ) 

^ ( i 

.0/ 

v ;:o 

*'? 


' ' • I Tr. 

( 1 

, 5 ) 



T \ > 

(■', 

) 

•/ <■* 
. . -J 


'HI T7. 

( 1 

, h ) 



*7/v ) 

(0, 

2) 

Tr 


*• a r - 

( 1 

, 7 ) 



r- m i 

(■), 

/ ) 

’1 


:.v/. 


(M,F,FP, 

XS, E, 1, I EL) 

- n r:-; 

C ! 

, •) ) 

X , F , F P 


' • 1 1 T ! 

< 1 

, i > 

) 


I r ■ 

C I 

, 1 1 

) 

- 


:• i i ( i, o OH 
IF (00- 1 ) V), ’0, ‘‘3 

.;•) ■ ' ! T I.JV : 

' ' - ( ( M -* ( ( * l . 7 ; .•> • ; + 0 3 ) / ( ?. A . • > ) ) * * 0 . 3 
- ' (IT'. ( 1 , 1 ) .i ■■ 


: < ; ' 'T I ■ ; V", 


E12.5) 


v , • r ;r i "'••••;! ( ;,nf>) 

/ i i / / i , o , o , n 
• l . • ' )■» 

f= \* <- c " (•:-• >( l . o/ :•;•())) -c 

f' 3 - . t - ( r. '(•;*(( l ..)/.( ( >- 1 • ) ) ) ) 
F.T IRM 



REPRODUOBIMTY OF THE ORIGINAL PAGE IS POOR. 


SAMPLE RUN OF /ROOT/ 


READ IN ENERGY 
3.35E05 


READ IN y0LUME 
7.82E-02 


EM= . 33500E+G6 


V0 = .73200E40! 


i 

! 

j READ IN XS 

! 1.0 EOS 


j READ IN E 

I .1 

i 

i 

IF A D IM ‘ I 

I) 1 ' 


. ?-)^l + F = •1256 1E-03 FP= • Z *8:W, 
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CALCULATION OF "W GOODNESS OF FIT STATISTIC FROM REFERENCE 36 

Table of Coefficients (a^+i) used in W test for normality for n = 9 
0.5888, 0.3244, 0. 1976, 0.0947 


Step 1 Rearrange observations for ordered Sample X ^ . . X r from 

percentile values in Table V 


Step 2 Compute 


s 2 = s (Xi-X> 2 = z xi‘- (~) 
i-l i- 1 

where X is data mean. 

Step 3 For when n is odd, set k = (n- i)/2 then compute 


n 


.L 


2 X i 


= A n <V X 1 ) + A n- 1 (X n- 1 ' X 2 ) + * * * + A n-k+l (X n-k+l 


(X_ 


= ^ A n -1+1 (X n-l+l ‘ V 

i-l 


_ ir _ 

Step 4 Compute the test statistic W — 


Step 5 Compare the calculated value of W with the percentiles of the 

distribution of the test statistic. From Table X of reference 36, 

n 1% 2% 5% 10% 50% 

9 0.764 0.791 0.829 0. 859 0.935 

This table gives the minimum values of W that we would obtain 
with 1, 2, 5, 10 and 50% probability as a function of n, if the 
data actually came from a normal distribution. Thus, small 
values of W indicate non- normality. 

1. Calculation of W Statistic for Fragment Parameters R, W, A, A/W 

Event 1 Distance 


200 


Normal Distribution 


Event 


Event 


Event 



1,910, 052 


(3920) 2 

9 


1,910,052 - 1,707,377.778 = 202,674.222 


b = 0. 5888 (722-222) + 0.3244 (608-287) + 0. 1976 (521-331) + 
.0947 (449-372) 


b = 443. 1683 

b 2 = 196398. 1421 

b^_ _ 196,398. 1421 

W “ 2 202,674. 222 

s 


= .969 


Log Normal 
S 2 = 326. 895411 


(54. 14941195) 2 = 326 . 8,5411 - 
9 


325. 7954239 = 1. 0999871 

b = 1. 045378098 

b 2 = 1= 092815368 
k 2 

W = — r — = . 993 

s 2 

2 Distance (Normal) 

S 2 = 185, 405. 5556 

b - 426. 695 

3 Distance (Normal) 

S 2 = 2, 330, 136. 889 
b = 1481. 4121 

4 Distance 

S 2 * 628, 382.889 
b = 785. 0769 


b 2 = 182068.623 
W = . 982 

b 2 = 2, 194, 581. 81 
W = . 942 

b 2 * 616, 345. 7389 
W * . 981 


01 


Event 5 Distance 


s 2 = 1,748,338.889 

b 2 = 

b = 1300. 255 

w = 

Event 6 Distance 

S 2 = 44172 

b 2 = 

b = 208.9122 

W = 

Event 7 Distance (Normal) 

S 2 * 47, 296 

b 2 = 

b = 214.0407 

w = 

Event 8 Distance (Normal) 

S 2 = 38,674. 8889 

b 2 ^ 

b = 193.7642 

W = 

Event 1 Weight (Ln Normal) 

S 2 = 4. 02243396 

b 2 

b = 1. 959105363 

w : 

Event 2 Weight (Ln Normal) 

S 2 = 9. 55743464 

b 2 

b » 3. 087710481 

w 

Event 3 Weight (Ln Normal) 

3 2 = 8. 308249 

b 2 

b » 2. 808424493 

w 

Event 4 Weight (Ln Normal) 

S 2 = 5.4833992 

b 2 


1,690,633.065 

.967 

' 43644.30731 
.988 

= 45813.42126 

* .969 

= 37, 544. 5652 
: .971 

= 3.838093823 
= .954 

= 9. 533956014 

* .998 

= 7.887248133 
= .949 

= 5. 422546817 


b = 2. 328636257 


W * . 989 



Event 5 Weight (Ln Normal) 

S 2 = 5. 2896359 
b = 2. 207715255 
Event 6 Weight (Ln Normal) 

S 2 = 18.40217547 
b = 4. 27189715 
Event 7 Weight (Ln Normal) 

S 2 = 12. 28668662 
b = 3. 391847545 
Event 8 Weight (Ln Normal) 

S 2 = 16. 56249558 
b = 4. 018598852 
Event 1 Area/Weight (Normal) 
S 2 = 11, 152.03556 
b = 104. 89576 

Event 2 Area/Weight (Normal) 
S 2 = 11639.72 
b = 105.62921 

Event 3 Area/Weight (Normal) 
S 2 = 1012. 18 
b = 28. 43535 

Event 4 Area /Weight (Normal) 
S 2 = 1064. 262222 
b * 29. 01741 


b 2 = 4. 874006647 
W = . 921 

b 2 = 18. 24910526 
W = . 992 

b 2 = 11. 50462977 
W = . 936 

b 2 = 16. 14913673 
W = .975 

b 2 = 11, 003. 12047 
W = • 987 


b 2 = 11157. 53001 
W = . 959 

b 2 = 8C8. 5691296 
W = . 799 

b 2 = 842. 010083 1 
W = .791 


* 


203 


Event 5 Are a /Weight (Normal) 
S 2 = 515. 8288889 
b = 18.61245 

Event 6 Area/Weight (Normal) 
S 2 = 16338.98741 
b = 126.724354 

Event 7 Area/Weight (No. . ) 
S 2 = 11,373.60751 
b = 105. 1395562 
Event 8 Area/Weight (Normal) 
S 2 = 19, 558.42122 
b = 138.6416769 
Event 1 Area (Ln Normal) 

S 2 = 5. 4012421 
b = 2.317066189 
Event 2 Area (Ln Normal) 

S 2 = 29.313013 
b = 5.322848019 
Event 3 Area (Ln Normal) 

S 2 = 11. 8267759 
b = 3.380014318 
Event 4 Area (Ln Normal) 

S 2 = 11.6117262 
b = 3. 350700597 


b 2 = 346. 423295 
W = .672 

b 2 = 16059.0619 
W = .983 

b 2 = 11054.32628 
W 85 .972 

b 2 = 19221.51457 
W = . 983 

b 2 - 5.368795724 
W = . 994 

b 2 * 28.33271103 
W ■ . 966 

b 2 = 11.42449679 
W « . 966 

b 2 ■ 11.22719449 
W - . 967 
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Event 5 Area (Ln Normal) 


S 2 = 5. 5540099 
b = 2. 350202139 
Event 6 Area (Ln Normal) 

S 2 = 15.45026354 
b = 3. 821539598 
Event 7 Area (Ln Normal) 

S 2 = 13.00914442 
b = 3. 586886773 
Event 8 Area (Ln Normal) 

S 2 = 14.3944438 
b " 3. 752708615 

2. Calculation of W Statistic for Initial 

CBM L0 2 -LH 2 

S 2 = 4. 1582624 
b = 7.961656617 
CBM LC> 2 -RP1 

S 2 = 2. 217937 
b ■ ). 449243666 
CBGS L0 2 - LH 2 

S 2 = 3. 4304869 


b 2 = 5.523450094 
W = . 994 

b 2 = 14.6041649 
W = .945 

b 2 = 12.86575672 
W = .989 

b 2 = 14.08282195 
W = . 978 

Velocity Distributions 

b 2 = 3. 848096683 
W = . 925 

b 2 ** 2. 100307203 
W = . 947 

b 2 = 3. 338371974 


b “ 1. 827121225 


W = .973 



CBGS LO z -RPl 



707618 


b = 1. 299335448 


b 2 = 1. 688272606 
W = .989 
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CALCULATION OF APPROXIMATE PROBABILITY FOR OBTAINING THE 

CALCULATED VALUE OF "W" 


For n = 9, from Table XI of reference 36, a = -2.968, n = 1.400, e - 0.3900 


Approximate probability of obtaining the calculated value of W, assuming a 
normally distributed variable can be obtained by finding: 




z= Y + 

Z = -2. 968 + 1. 400 In 


/w - 0 . 3900\ 

\ I 


For W = . 942 
Z = -2. 968 + 1. 400 In 




968 + 1. 400 In 9. 517241379 

= -2. 968 + 1. 400 (2. 253 105036) = -2. 968 + 3. 15434705 = . 186347 
Pr (Z < 0. 186) = . 57 
For W = .994 
Z = 3. 4885 


Pr = . 999 
For W = .969 
Z = 1. 13 


Pr - . 87 1 


For W = .954 
Z = .54 
P = .7054 

The following values and table result from using the above methods. 



4 

> 


W = .921 
W = . 9 " 
W = . ^ 

W = . 

W = . 1 
w - 
w = . 

W - . 9V* 


Z * -.3 
Z = . 186 
Z = . 268 
Z ■ 1. 13 
Z = . 540 
Z = 1.63 
Z = 2.628 
Z = 3.4885 


P (Z^ -.3) = .382 
P(Z<. 186) = .574 
P (Z< . 268) = .643 
P (Z< 1. 13) * .871 
P (Z< . 540) = .705 
P (Z< 1.63) = . 948 
P (Zf 2.62) = . 995 
P (Z< 3.488) = 9. 99 


Summary of the Above Results 


w 

% 

764 

.010 

791 

. 020 

829 

.050 

859 

. 100 

921 

.382 

935 

. 500 

945 

.643 

969 

. 871 

978 

.948 

988 

.99 
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APPENDIX H 


COMPUTER PROGRAM ENTI'-i-ED /TEMP/ IN FORTRAN IV 

Function: This program computes the range of a fragment from the equations 
described in Section IV.D for the following input data. 


SI 

Initial trajectory angle of fragment 

radians 

S2 

Fragment drag coefficient 


S3 

Ambient air density 

lb -ft 3 

S4 

Fragment mass coefficient (ratio of fragment 
mass to tank mass) 


TM 

Tank mass 

lb-sec 2 /in. 

FF 

Projected cross-sectional area of fragment 

in 2 

G2 

Fragment initial velocity 

ft/ sec 


Variables: The definition and units of the variables in this program are given 


in the following table. 



Program 

Variable 

Variable 

Definition 

Units 

Si 

0 

Initial trajectory angle 

radians 

S2 

C D 

Drag coefficient 


S3 

P. 

air 

Air density 

lb -ft 3 

S4 


Mass coefficient 


TM 

M 

Fragment mass 

lb- sec^/in. 

FF 

A 

Fragment cross-sectional area 

. 2 
in 

G2 

u f 

Initial velocity 

ft/ sec 

PI 

TT 

The constant n 


GR 

g 

Gravitational constant 

/ 2 
in/ sec 
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Program 

Variable 

Variable 

Definition 

Units 

S5 


The quantity 2A / 

. 2 
in 

S6 

c 

Coefficient defined in Section 

. -1 



IV. D. 

in 

S7 

V Ro 

Initial radial velocity 

ft/ sec 

S8 

V 

zo 

Initial vertical velocity 

1/2 

ft/ sec 

S9 


The quantity V^ o (c/g) 

in. 

S10 


-1 / i \ 1/2 

The quantity tan ^ z0 ' c ' 8' 

1/2 

The quantity 1/ (eg) 

1 / 2 

The quantity cos tan" * V^ (c/g) 


Sll 


sec 

S 12 



S13 


The quantity 2.0 LOG (1. 0/S12) 


S14 


The quantity 2. 0 - 1. 0 


S15 

T 

Time of flight of fragment 

sec 

S16 

R 

Fragment range 

ft 

TR 

*R 

Time of rise of fragment 

sec 

ZM 

Z m 

Maximum height reached by fragment 

in. 

S 17 


The quantity 2. 0 e^ c Z m - 1.0 

sec 

TF 

*f 

Time of fall of fragment 

sec 

T 

T 

Time of flight of fragment 

sec 


Sample runs: Substitution of the following data in the program yielded the range 
values R (appearing in Figure (J|- 1) as X's). Thus, results of this program 
are in accord with Oslake, et al. 
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RANGE OF FRAGMENT feet 


FIGURE 



HI FRAGMENT RANGE VERSUS ELEVATION ANGLE 

FOR CONSTANT JL = 100 LB /FT 2 , AND 

CONSTANT u, C A 




TABLE H- 1 - DATA FOR PROGRAM CHECK 


51 =43.3° 02 

52 - 1. 2 

53 = 7. 48 x 10' 2 lb/ft 5 

54 = 1. 0 

FF = 4. 52 x 10 3 in 2 
TM s 6. 2 lb-sec / In. 


: 100 ft/sec 

R = 2. 89 x 

10 2 ft 

500 ft/ sec 

3. 19 x 

10 3 ft 

1000 ft/ sec 

5. 50 x 

10 3 ft 

5000 ft/ sec 

1.09 x 

10 4 ft 


A sample run for the case in which a "mean" fragment from test 062 of PYRO 
was considered ICq = .750) is given in the following. 

READ IN TRAJ. ANGLE 
1 2.77 


READ IN DRAG C0EF. AND AIR DENSITY 
.75, 7. ARE-02 


READ IN FRAGMENT MASS C0EF • 
1 .0 
1 .0 


READ IN N0 . 0F FRAGMENTS 


FRAG. MASS 
1 .7/SF-02 


ERAG. PR0J • AREA 
5 .77E02 


INITIAL VEL. 

7 .4IE02 

AMRTENT AIR DENSITY- .748C-01 LBS. /CUBIC ET . 



AIR DRAG rfoFF 


7S0F+O0 


FRAGMENT MASS G0EF.= . 100F+01 
TRAJECTORY ANGLE 3 . 188F+08 DEGREES 


check r- .H7SF-0R 
782 • AT 1 5AO/> 
l 63 .7R90A93 
2.9SM 1 SA6R 
1 .2AA59RR02 
. S«3S^ F+ 01 


MAXIMUM RANGF= .353E+03 FT. 
. 52J84E+01 
* STOP* 


.MAT 

' n \ \T 


<2/,23H REAL) IM TRAJ . AUGL£> 

<.vl35* > >£.«> IS DIUG CW. AMD A1S 8EKSITY) 


! V 

t \T 


( vi 3 

c i i 

c *4 / .. 


3 ) 

IN F 3 AG -1 1*1 '4 T 1 A 3 1 

I : «-«diti i:-js t.j truJectvry) 

( / , 2 1 ;i aM J IEUT AIR DENSITY* / *L 1 0 i 3 
C/,16H AIR DRAG CMEF . = > E 1 0 . 3 ) 


C3EF • > 


1 711 


i A 

,..T 


( / 3 1 

«; 

yr a 

In'. : .-. IT MASS 

( /, 1 

- 

! [ 


TTRY A3 

(3/ 


> 0 

CHECK 

C* / E 1 ' 

( 3 / , 

1 

3;{ 

MAXIMUM Ra'L' 

(V, 

j 

r > 3 

READ 

in :n . 

( Vi 

1 

7 H 

FRAG . 

PRO J . 

( / i 

1 

1 3 

“ RaG . 

1 A S 3 ) 

( /> 

1 

3 . i 

r 1 1 ? i 

AL ’ VL 


13.3) 

= , L 1 ' 3 . 3 , Si ; OF GRF, if. 3 ) 


i F FRAGMENT 3 ) 
.VEA) 




i i : 1 V , ,‘r) 


( 

! , 3 DA ) 


) 

,1131) 31 


> 

I 


i 

. 1 * 3 1 ) / ( • 

! .:J*PI ) 

( 

1 , V) 1 ) 


0 

, ; > i i ) .>• 

>,33 


( 1 , 33 ? ) 
(0,33)1) 3 0 

(1,311) 

' , , 3 301 ) FO 



*/****/• 


READ (0,3001) TM 

write (1/20D 

READ (0/3001) rr 
WRITE (1/202) 

READ (0/3001 ) G2 
WRITE (1/30A) S3 
1 WRITE (1/305) S2 

\ WRITE (1/306) SA 

| WRITE (i/307') SI 

t S5 = (2.0*(FF/PI) ) „ ■ 

\ S6=(0.5*(( (S3/GR)*S5*S2)/<SA*(TM/FN >>>>/( 12*0 **3.0) 

l ’.TRITE ( 1/ 3071 ) S6 

' S7=G2*C3S< <2.0*PI >*(Sl/360.0) ) 

* S8= G2*S IU ( ( 2 . 0*P I > *( S 1 / 3,60.0 ) ) 

i S9= ( S8* ( (S6/GR>**0.5) ) * ( 12.0) 

1 S 1 0 = AT AW ( S 9 ) 

£ SI 1 = ( 1 .0/ ( S6*GR) ) **0. 5 

| 512=COS(S10) 

1 513=2. 0*AL3G( 1 .0/512) 

• \ 3U*(2.0*EXP(S1 3) )- 1 .0 N 

315 = 51 1 *(S 1 0+ ( 0 . 5 *AL0G ( S 1 4*(((S1A**2.0)-1.0)**0.5>)>) 
’.-RITE (1/300 1) SI 5 

SI 6= ( 1 . 0/S6 ) *AL03 ( 1 .0 + (S6*S7*Sl'5*12.0)) 

. S 1 A = S 1 6 / 1 2 . 0 

WRITE (1/303) 516 
tr = s l 1 * Maw ( S') ) 

’ != ( 1 . 0/SA) *AL.1G( 1 . 0/ CO S ( TR/ S 1 1 ) ) 

51 7= (2. D*£XR(2.0*S6*ZM) )-l .0 

TF= 0.5*31 1 *AL0G(317+( ( ( SI 7**2 .0 )- 1 .0 ) **0 . 5 ) ) 

7=T? + TR 

'".I T 5 (1,333 1) T 


appendix I 


Symb ol 
a, a 

o 

LIST OF SYMBOLS 

Definition 
Sound speed 

* 

Units 

LT 

** 

a,b, c, d 

Parameters in a math model 


b 

Time constant of blast wave 

— 

d, D, D 

o 

Diameters 

L 

dc A' dC B 

Distances to cameras 

L 

g 

Acceleration of gravity, nondimens ional 
displacement 

LT , - 

h 

Convective heat transfer coefficient, segment 
height 

MT"V 

h 

e 

Height of explosion center 

L 

k 

Discharge coefficient, constant in fragment 
drag equation 

, ML 

j e 

Generic length 

L 

n, N, N q 

P > P 
*o oo 

q 

Numbers of fragments 
Ambient pressure 
Drag pressure 

FL" 2 

FL' 2 

r 

t.T,T i ,' r 

Displacement of a fragment 
Times 

L 

T 

u , u g , U , U j, 

U,. u„, u„ 

Velocities 

lt" 1 


>!<>',* 


" Units u..d are fore. (F). lengtML ). m... <M>, time <T> and temperature (0) 
rather than dimensions in any particular physical system. 

A dash in the units column indicates a dimensionk as quantity. 



Symbol 

w 



F 

G., G.' 
1 1 


H e' H R’ H TNT 

H f, 

k 

H b. 

L 

^d 

K. 

1 

L 

M, M , M_ 
i T 

p, P r ,Q 

R 


s. 

I 


T 


Definition 

Crack width 

Mixing function 

Coordinates of fragments 

Yield Fraction 

Areas 

Mass of explosive charge or confined ga3 

Drag coefficient 

Energy of reactants' 

Projected area of a fragment 

Dimensional parameters fixed by nondimen 
sional groups 

Heats of explosion 

Heats of fusion 

Heats of boiling 

Blast wave impulses 

Thermal conductivities 

Tank length 

Masses 

Peak pressures 
Distance from accident 
Universal gas constant 
Displacements 
Temperature 


Units* 


L 


L 


L 


2 


M 


FL 

7 

L“ 


various 


r 2 _-2 

L T 


T‘ 


L 2 T" 2 
FTL -2 

mlt" 3 c- 1 


L 


M 


FL 


-2 


L 

l 2 t -2 0“ 1 


L 


0 


2 16 


Definition 


v 


Kinematic viscosities 


Time ratio 


Crack length 
Densities 

Volumetric heat capacities 
Stefan -Boltzmann constant 
Standard deviation of velocity 
Solid angle 


L 

ML -3 
ML" 1 
MT 3 0 



Azimuth angles 
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